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SUMMARY 


Some  of  the  techniques  in  current  use  for  measuring  dynamic  stability 
derivatives  in  wind  tunnels  are  described,  with  emphasis  given  to  the 
Important  features  of  balance  system  design,  data  reduction  methods, 
Instrumentation  and  typical  balance  systems.  The  use  of  gas  bearings 
for  dynamic  stability  and  roll  damping  balances  is  treated  and  a  thre > 
degree-of- freedom  balance  system  employing  a  spherical  gas  bearing  is 
described. 


RESUME 


Qielques-unes  des  techniques  courantes  pour  la  mesure  des  d^rlvees 
de  la  stabilite  dynamique  dans  les  tunnels  ae'rodynamiques  sont  de'crivees 
en  soul ign ant  les  particularity  importantes  du  dessin  des  systemes  de 
compensation,  des  methodes  pour  la  r&iuction  de  donne'es,  1’  instrumenta¬ 
tion  et  des  systemes  de  condensation  typiques.  On  trait  1* emploi  des 
paliers  k  gaz  afln  d'obtenir  de  stablllte"  dynamlque  et  des  compensateurs 
pour  1'  amort issement,  de  roulis  ainsi  qu’ un  syst&me  de  condensation  en 
trols  degree  de  llberte  avec  1’ utilisation  d’ un  palier  a  gaz. 
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SECTION  1 
INTRODUCTION 


Analysis  of  aircraft  and  missile  aerodynamic  performance  and  stability  requires 
knowledge  of  the  forces  and  moments  which  act  under  conditions  of  steady  and  unsteady 
flight.  One  of  the  first  studies  of  aircraft  stability  was  that  of  Lanchester1  around 
1900,  followed  later  by  more  complete  and  rigorous  analysis  by  Bryan2,  which  is  still 
the  basis  for  much  of  the  dynamic  stability  work  today. 

Requirements  for  supplying  aerodynamic  characteristics  for  performance  analysis 
have  been  met,  to  some  extent,  by  relying  on  theories  supported  by  experimental  measure¬ 
ments  to  establish  the  validity  of  the  theories.  In  some  cases,  the  theory  has  only 
limited  application  especially  when  viscous  flow  has  a  pronounced  influence  on  the 
flow  field,  as  in  the  case  of  high  speed  flow  (hypersonic  speeds  and  above)  and  low 
Reynolds  numbers.  There  are  some  situations  for  which  the  phenomena  cannot  be  described 
adequately  by  theory.  In  other  cases  theory  is  not  applicable  because  the  configuration 
and  flow  field  are  so  complex  that  the  many  approximations  required  invalidate  a  wholly 
theoretical  approach.  Thus,  in  practice,  strong  reliance  is  placed  on  experimental 
results  for  use  in  analyzing  the  modes  of  motion  of  the  aircraft  or  missile. 

Early  experimental  techniques  making  use  of  oscillating  models  to  dynamically  simulate 
rigid  modes  of  motion  for  determining  stability  derivatives  are  well  summarized  hy 
Jones3.  It  is  of  interest  to  note  that  no  fundamental  changes  in  the  basic  experimental 
methods  have  taken  place  since  that  time,  1934,  although  much  work  has  been  done  toward 
refinement  of  the  methods,  equipment  and  instrumentation.  Little  emphasis  was  placed 
on  dynamic  stability  experiments  until  about  1940,  and  from  that  time  to  1945  experi¬ 
mentation  was  primarily  concerned  with  sui  sonic  flow  and  tests  were  made  for  small 
perturbations.  Wind  tunnel  experiments  and  flight  tests  within  the  early  forties  were 
directed  more  toward  static  stability  and  control  problems.  With  the  advent  of  high 
speed  aircraft,  missiles,  rockets,  and  re-entry  vehicles,  Increased  emphasis  has  been 
placed  on  performance  and  dynamic  stability  problems  since  missions  demand  precise 
control  not  only  in  small  disturbances  from  level  flight,  but  also  in  large  scale 
maneuvers.  As  a  result,  dynamic  stability  measurements  have  become  increasingly  important 
at  supersonic,  hypersonic  and  hypervelocity  speeds  and  at  large  perturbations.  The 
experimental  methods  suggested  by  earlier  experiments  have  been  used  in  different  forms, 
and  many  improvements  have  been  made  in  the  experimental  techniques  since  about  1950 
These  improvements  have  been  in  the  areas  of  balance  system,  model  and  support  system 
design,  instrumentation,  data  ac<*jisition,  and  data  reduction. 

In  1954  Valensi4  presented  a  review  of  the  techniques  primarily  in  use  in  Europe  for 
measuring  oscillatory  aerodynamic  forces  and  moments  on  models  oscillating  in  wind 
tunnels.  The  following  year  AGARDograph  11  (Ref.  5)  by  Arnold  was  published  on  the 
subject  of  dynamic  measurements  in  wind  tunnels.  The  most  recent  summary  of  techniques 
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for  measuring  oscillatory  derivatives  in  «ind  tunnels  was  published  by  Bratt6  in  1963 
as  a  part  of  the  AGARD  Manual  on  Aeroelasticity.  The  paper  discusses  the  basic  principles 
employe*  in  aeasuring  derivatives  and  gives  some  account  of  the  associated  instrumentation. 
Much  of  the  report  deals  with  the  methods  and  instrumentation  employed  in  England. 

Since  these  earlier  works,  the  emphasis  on  dynamic  stability  testing  has  increased, 
and  many  reports  have  been  prepared  on  the  subject,  so  there  exists  a  need  for  a  review 
and  summary  of  the  current  dynamic  stability  testing  techniques  in  use  at  supersonic 
speeds  and  above. 

There  are  many  experimental  techniques  which  are  suitable  for  measuring  dynamic 
stability  derivatives;  however,  only  those  techniques  that  are  in  most  common  use  in 
wind  tunnels  will  be  described.  In  addition,  many  different  systems  exist  for  measuring 
the  same  quantities  and,  since  they  differ  primarily  in  the  specific  application,  no 
attempt  will  be  made  to  describe  all  systems;  only  systems  considered  to  be  representa¬ 
tive  and  those  having  unique  features  will  be  discussed.  The  report  is  organized 
primarily  around  the  aspects  of  dynamic  stability  testing  in  wind  tunnels  at  supersonic 
speeds  and  above,  although  much  of  the  report  is  pertinent  to  aspects  of  testing  at 
transonic  and  subsonic  speeds.  Although  free-flight  testing  in  ground  test  facilities 
encoapasses  the  use  of  the  wind  tunnel  and  aeroballistic  ranee  for  measurements  of 
dynamic  stability  derivatives,  these  techniques  are  not  included  as  a  part  of  this 
AGARDograph  because  other  publications  on  these  subjects  have  been  published  and  are 
in  preparation7,  8. 
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SECTION  2 
STATEMENT  OF  THE  PROBLEM 


2. 1  DEVELOPMENT  TRENDS 

Aircraft  or  missiles  in  high  speed  flight  are  free  to  respond  to  disturbances  with 
motions  involving  many  degrees  of  freedom.  During  various  stages  of  aerodynamic  design 
of  a  vehicle,  as  the  configuration  is  developed,  dynamic  characteristics  must  be  known 
and  progressively  refined  to  a  high  degree  for  the  final  configuration  performance 
specification.  The  designer  must  necessarily  be  concerned  with  the  interpretation  of 
aerodynamic  stability  data  and  with  the  evaluation  of  the  influence  of  their  accuracy 
on  an  integrated  system  design. 

There  are  three  methods  of  obtaining  stability  data: 

(i)  Theory  and  empirical  data. 

(ii)  Wind  tunnel  and  aeroballistic  range  model  testing. 

(iii)  Sub-  or  full-scale  flight  testing. 

During  the  preliminary  design  stage,  it  is  most  expedient  to  rely  primarily  on  theory 
and  empirical  data  for  performance  evaluation  whereas,  in  the  later  stages,  design 
concepts  inevitably  require  evaluation  through  the  U3e  of  some  scale  model  tests.  As 
a  final  phase  in  the  design  development.,  sub-  or  full-scale  flight  tests  are  required 
to  verify  performance  and  obtain  in-flight  measurements  for  verification  of  theory  and 
prediction  methods. 

There  are  certain  speed  regimes  where  theoretical  results  may  be  inaccurate  due  to 
insufficient  or  inaccurate  basic  data.  For  example,  in  the  transonic  speed  regime 
(0. 9  <  <  15)  theoretical  methods  for  determining  dynamic  stability  derivatives  are 

very  limited,  and.  therefore,  reliance  is  placed  on  methods  involving  the  use  of 
empirical  data  on  similar  configurations  and  correlation  plots  showing  the  manner  in 
which  the  stability  derivatives  vary  with  Mach  number  for  similar  configurations. 

Within  the  supersonic  and  hypersonic  region  (1  5  <  Me  '  10).  theoretical  methods 
are  available  for  estimating  derivatives  for  bodies  of  revolution  and  relatively  simple 
lifting  configurations.  At  Mach  numbers  above  10  in  the  hypervelocity  speed  regime, 
theories  generally  fail  to  predict  the  derivatives  because  viscous  effects  have  a 
strong  influence  on  the  vehicle  aerodynamic  characteristics 

In  order  to  improve  the  accuracy  of  aerodynamic  derivatives  that  are  used  in  per¬ 
formance  predictions  and  to  validate  theoretical  estimates,  model  tests  in  ground  test 
facilities  are  recognized  as  being  essential  to  the  development  of  aircraft,  missiles 
and  re-entry  vehicles. 
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Study  of  vehicle  motions  in  flight  involves  consideration  of  many  degrees  of  freedom. 
When  the  vehicle  is  considered  as  a  rigid  body,  six  degrees  of  freedom  are  required. 
Since  the  vehicle  is  not  a  rigid  structure,  additional  relative  motions  of  components 
such  as  aeroelastic  deflection  of  lifting  surfaces  and  movable  controls  constitute 
additional  degrees  of  freedom. 

The  vehicle  motions  may  be  separated  into  high  frequency  and  low  frequency  motions, 
where  high  frequency  motions  studies  such  as  flutter  are  primarily  concerned  with 
structural  elasticity  and  unsteady  aerodynamics.  At  low  frequencies  the  motions  involve 
rigid  body  motions  and  are  studied  as  dynamic  stability  using  quasi- static  aerodynamics. 
The  mode  of  model  testing  discussed  in  the  present  AGARDograph  is  concerned  with  this 
latter  type  motion  and  the  determination  of  rotary  damping  derivatives. 

In  scaling  down  a  vehicle  for  model  tests  the  parameter  that  must  be  considered  in 
addition  to  the  Mach  number  and  Reynolds  number  is  the  angular  velocity  or  frequency 
of  model  oscillation  in  relation  to  full-scale  conditions.  The  scaling  parameter, 
called  the  reduced  frequency,  in  the  form  u.d/2Va<  contains  the  frequency  of  oscillation 
and  represents  the  ratio  of  some  characteristic  dimension  of  the  vehicle  to  the  wave 
length  of  the  oscillation.  This  type  of  scaling  parameter  also  applies  to  dynamic  roll 
tests  and  has  the  form  Pd/2V00  . 

Exper)  mental  techniques  for  measuring  dynamic  stability  derivatives  in  ground  test 
facilities  may  be  classified  in  many  different  ways.  The  classification  that  will  be 
followed  here  is  presented  in  Figure  1.  The  two  major  classifications  are  derived  on 
the  basis  of  whether  the  model  is  free  to  move  in  all  basic  degrees  of  freedom  or 
whether  it  is  restricted  in  several  degrees  of  freedom.  The  subtopics  shown  in  Figure  1 
for  discussion  were  selected  as  those  methods  currently  being  employed  for  most  high 
speed  testing;  however,  it  should  be  emphasized  that  the  outline  does  not  show  the  case 
of  the  model  fixed  in  the  test  section  with  perturbations  generated  in  the  flow.  As 
noted  in  the  Introduction,  only  those  methods  which  apply  to  captive  model  testing 
techniques  will  be  discussed  in  this  AGARDograph  since  free-flight  tests  in  the  wind 
tunnel  and  range  are  the  subject  of  AGARDographs  published  and  now  in  preparation7,8. 
Additional  general  information  on  methods  not  described  herein  may  be  found  in  Reference  9 


2. 2  EQUATIONS  OF  MOTION 

In  order  to  place  the  role  of  dynamic  stability  testing  in  perspective,  it  is  helpful 
to  examine  the  equations  of  motion.  Although  the  mathematical  treatment  of  vehicle 
dynamics  is  well  known  and  found  in  many  tests10*12,  a  brief  derivation  of  the  equations 
of  motion  is  included  to  illustrate  the  basis  for  the  experimental  techniques  used  in 
the  study  of  dynamic  stability. 

2. 2. 1  Axes  Systems 

The  equations  of  motion  will  be  developed  using  the  axes  systems  shown  in  Figures  2 
and  3.  Model  orientation  is  determined  relative  to  the  tunnel  fixed  axes  XTYTZT 
(inertial  reference).  The  XYZ  axes  are  fixed  in  the  model  with  the  X  and  Z  axes 
in  the  model’s  plane  of  symmetry.  When  the  origin  of  this  system  lies  at  the  center 
of  gravity,  these  axes  are  then  defined  as  the  body  axes.  The  angular  orientation  of 
the  XYZ  system  with  respect  to  the  XTYTZT  system  and  inertial  space  is  given  by  the 


5 


Euler  angles  ^  ,  6  ,  and  *>  .  Sometimes  it  is  convenient  to  use  the  nonrolling  axes 
system  (XYZ).  These  are  special  body  axes  whose  angular  orientation  in  space  relative 
to  the  XtYtZt  system  is  determined  by  the  angles  and  0  .  This  axis  system  can 
be  used  even  though  the  model  is  rotating  about  the  X  axis.  The  origins  of  all 
systems  arc  assumed  to  occupy  the  same  point  in  space  (Pig.  3). 

2.2.2  Inertial  Moments 

The  angular  momentum  of  the  model  may  be  expressed  as 

H  -  /( r  x  Vp)  dm  ,  (1) 

where  r  is  the  radius  vector  from  the  origin  of  the  rotating  XYZ  axes  to  the  particle 
p  of  mass  dm  and  Vp  is  the  velocity  of  the  particle  p  relative  to  inertial  space. 

The  radius  vector  r  is  measured  relative  to  the  body  fixed  XYZ  axes;  therefore, 
in  order  to  evaluate  the  velocity  of  p  relative  to  inertial  space,  it  is  necessary 
to  use  the  well-known  transformation  for  the  rate  of  change  of  any  vector  from  fixed 
to  rotating  axes  as  follows: 


where 


Ou 


dr  dr 

—  =  —  +  Zi  x  r  , 

dt  T  dt 

=  Pi  +  Qj  +  Rk  . 


Assuming  the  model  is  a  rigid  body,  the  equation  reduces  to 

Vp  =  Z  x  r  , 

since  r  is  then  not  a  function  of  time. 

Substituting  this  equation  into  Equation  (1)  yields 

H  -  Jr  x  (cl  x  r)  dm  . 

This  expression  is  then  integrated  over  the  vehicle  to  obtain 

H  = 

[j]  - 

[I]  - 

W)  = 


[jILIU  . 

Li  JkJ 

h  ”jxy  "Jxz 
-Jyx  Iy  -JyZ 
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Y 
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The  moment  acting  on  the  vehicle  is  equal  to  the  rate  of  change  of  angular  momentum 

dH  dH 

—  =  —  +  o>xH  =  2||  . 

dt  T  dt 

This  vector  equation  can  be  reduced  to  the  following  scalar  forms: 

IXP  +  (lz  -  I,)RQ  -  JXZ(R  +  CP)  +  Jxt(RP  -  Q)  +  J,Z(R2  -  Q2) 

1,4  +  (Ix  -  IZ)PR  -  JX,(P  t  CR)  +  JYZ(P«  -  R)  +  JXZ(P2  -  R2) 

1^  +  (I,  -  IX)(JP  -  J,Z(Q  +  PR)  +  JXZ(QR  -  P)  +  JX,(Q2  -  P2) 

Roll  pitch  and  yaw  rates  for  the  XYZ  system  may  be  obtained  from  Figure  3  as 

follows: 

P  =  -  ¥  sin  0 

Q  =  S  cos  $  +  ¥  cos  0  sin  $ 

R  =  4*  cos  0  cos  $  -  0  sin  $  . 

In  order  to  sinplify  the  nonlinear  differential  equations  given  above  to  a  form 
which  facilitates  analysis  of  the  motion,  vehicle  symmetry  will  be  assumed,  and  the 
motion  will  be  restricted  to  small  disturbances  from  a  reference  condition.  Let 
P-p0+p,  ¥  =  <//Q  +  «//  ,  etc. .  where  the  zero  subscripted  quantities  are  the  reference 
conditions  and  the  unsubscripted  lower  case  quantities  are  perturbations  from  this 
reference  condition. 

For  the  following  discussion,  the  reference  condition  is  selected  such  that 
q0  =  r0  =  0Q  =  \pQ  =  0  and  the  reference  quantities  p0  and  <£0  are  not  necessarily 
small.  Considering  vehicles  having  symmetry  about  the  XZ  plane  (Jxy  =  Jy2  =  0)  and 
neglecting  products  of  perturbation  quantities,  the  equations  of  motion  become 

IXP  -  Jxz(r  ♦  qp) 

I,q  ♦  (Ix  -  Iz)Pr  t  jx,p2 

Izr  +  (i y  -  Ix)Pq  ~  J 

The  roll,  pitch  and  yaw  rates  at  this  condition  are  given  as 

Po 
P 

q 


=  <aox 

=  (2m)  ,  (y, 

=  (2m)  z  . 


=  <S“*> 

=  (2m,)  (2) 

=  (Si,). 


r 
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For  models  having  triagonal  or  greater  symmetry,  the  nonrolling  coordinate  system 
XYZ  may  be  conveniently  used  when  these  axes  correspond  to  the  models’  principal 
inertia  axes.  Since  the  body  is  moving  with  respect  to  the  Y  and  Z  axes  (rotating 
about  the  X  axis),  the  restriction  of  symmetry  is  necessary,  otherwise  the  moments  of 
inertia  with  respect  to  the  XYZ  system  will  be  time  varying.  For  the  above  case, 
the  angular  momentum  expression  is  given  as 

H  =  IYPi  +  I  of  +  Irk  . 


I  =  ly  =  IZ 


q  -  q  cos  v  -  r 


r  =  r  cos 


$  +  q  sin  $ 


^  f  P(t)  dt  ,  |$|  »  |-0  sin  6\ 

Jo 


Note  that,  even  though  $  =  0  for  the  nonrolling  coordinates,  its  XY  plane  is 
rotating  with  respect  to  the  XTYT  plane  at  the  angular  rate  -0  sin  6  . 

The  equations  of  motion  referenced  to  the  nonrolling  coordinates  are  given  as 


—  +  n  x  h  =  . 

dt 


=  Oi  +  qj  +  rk  . 


These  equations  of  motion  reduce  to 


Iq  +  IYPr  = 


I7r  IvPq 


2.2.3  Aerodynamic  Moments 

The  general  procedure  for  expanding  the  aerodynamic  moment  relations  is  that  deve¬ 
loped  by  Bryan2.  It  is  assumed  that  the  moments  are  functions  of  the  steady- state 
reference  conditions  and  the  instantaneous  values  of  the  disturbance  velocities,  con¬ 
trol  angles,  and  their  time  derivatives.  A  Taylor  series  expansion  of  an  aerodynamic 
reaction  yields  this  form.  For  example,  the  aerodynamic  reaction  A  =  A(U,V,  ...)  is 
expressed  as 
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Aerodynamic  reactions  can  be  expanded  with  respect  to  the  variables  U  ,  V  ,  W  ,  P  , 
Q  ,  R  ,  the  control  deflection  angles,  and  time  derivatives  of  these  variables.  For 
a  given  configuration,  conditions  of  synmetry  and  past  experience  concerning  relative 
effects  of  these  derivatives  can  be  used  to  reduce  the  number  of  derivatives  needed 
to  describe  the  aerodynamic  moment  system.  For  a  spinning  vehicle  having  triagonal  or 
greater  symmetry,  the  moment  equations  can  be  reduced  to 


The  term  Ollr/3v)0  in  the  pitch  equation  and  the  term  OMj/Bw)0  in  the  yaw 
equation  are  usually  negligible  except  for  the  case  when  the  vehicle  is  spinning  rapidly. 
These  terms  are  caused  by  the  Magnus  effect.  They  are  given  in  a  more  familiar  form 
as  follows. 

Let 

3v,  0 


0 
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Taking  the  first  tern  as  an  example  and  assuming  that  it  is  primarily  a  function 
of  P  ,  its  Taylor  series  expansion  about  the  condition  P  =  0  is  given  as 


Myv0(P....) 


mYv0(0,...) 


(P  -  0)  +  . . . 


The  first  term  on  the  right  hand  side  is  usually  negligible;  therefore,  this  expression 
may  be  reduced  to 


Similarly, 


my?0(P) 


The  pitching-  and  yawing-moment  equations  are  then  given  as 


The  terms  Mxo  ,  Myo  .  and  Mzo  in  Equations  (6),  (7),  and  (8),  respectively,  are 
the  moments  about  the  X  ,  Y  and  Z  axes  while  the  model  is  at  the  reference  flight 
condition.  The  reference  condition  is  often  chosen  such  that  these  terms  are  zero. 


Equations  (3)  describe  the  motion  of  a  model  which  has  three  degrees  of  freedom  in 
rotation  about  a  point  fixed  in  space.  Unique  test  equipment  is  required  in  order  to 
obtain  these  motions  using  a  captive  model  test  technique.  Two  balances  which  have 
this  capability  are  described  later  in  Section  3. 1. 1.  In  most  captive  model  dynamic 
stability  tests,  the  model  is  restrained  to  a  single  degree  of  freedom,  i.e. ,  rotation 
about  a  single  axis.  An  example  of  a  typical  dynamic  stability  balance  would  be  a 
single-degree-of- freedom  system  using  a  flexure  pivot.  The  equation  of  motion  for  a 
system  of  this  type  will  be  developed  in  the  following  discussion. 


2.2.4  Single  Degree  of  Freedow 

The  velocity  components  in  the  body  axis  for  the  case  of  a  model  mounted  on  a  single- 
degree-of-  freedom  flexure  pivot  at  an  angle  of  attack  a  and  oscillating  symmetrically 
at  a  small  angle  6  are 

U  -  v*  (a  ♦  £)  =  cos  a  -  sin  a 

V  =  0  • 


.0 


W  =  V,,,  sin  (a  +  9)  =  9\go  cos  a  +  V*  sin  a ! 
P  =  0 

Q  =  q  =  £  ,  q  -  9 

R  =  0  . 


(9) 


The  reference  conditions  used  here  are  those  that  exist  when  the  model  is  at  the  angle 
of  attack  a  .  Under  these  conditions,  Equation  (2)  reduces  to 

I,Q  =  <»)Y  . 

Only  the  derivation  of  the  equation  for  a  single  degree  of  freedom  in  pitch  is  presented, 
since  the  derivations  for  oscillations  in  yaw  and  roll  are  similar  and  are  treated  in 
detail  in  many  of  the  references  included  (see,  e.g.  ,  Reference  4). 

An  elastic  restraint  is  provided  in  the  balance  to  aid  in  restoring  the  model  to  the 
reference  attitude  following  a  displacement.  When  the  aerodynamic  moment  is  expanded 
to  Include  all  significant  terms  and  the  restoring  moment  Kt(a  +  6)  and  the  damping 
moment  C of  the  restraint  are  added  to  this  to  obtain  the  total  moment,  the  equation 
of  motion  in  pitch  is  written  as 

IyQ  =  Myo  +  MyU  +  Mqq  +  M.W  +  M**  +  K,(0  +  a)  +  C ,£  . 

where 

=  (3Nr/du)0  .  Hq  =  (3my/3q)0 

h,  -  0«y/3*)#  ,  M*  =  Omy/3*)0  . 

After  including  the  relationships  for  the  linear  and  angular  components  of  velocity 
(Equation  (9)),  the  equation  above  reduces  to  the  form 

1y9  +  (-Mq  -  Vajll*  cos  a  -  Cj)<9  + 

+  (Voo  8in  a  -  V®  C0S  a  -  K1>£  “  MY0  +  Kia  =  0  .  ( 10) 

In  the  undisturbed  case,  conditions  of  equilibrium  give 

myo  "  Kia  • 

and.  for  small  values  of  a  ,  Equation  (10)  becomes 

I,fi  -  Mq  -  V,,!!*  -  C,)0  +  <-K,  -  M.VJ0  =  0  . 


(11) 
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Defining  the  damping- moment  parameter  and  the  restoring-moment  parameter  by  and 
Mg  ,  respectively,  where 


Mga  =  (Mq  +  V^)  ,  aerodynamic 

damping- moment  parameter  of 
elastic  restraint  considered 
as  a  tare 

aerodynamic 

static-moment  parameter  of 
elastic  restraint  considered 
as  a  tare. 

With  these  definitions,  Equation  (11)  becomes 

id  -  <Mjt  t  -  (Met  +  M9a)£?  =  0  . 

Following  the  procedure  presented  in  Reference  4,  similar  expressions  can  be  written 
for  the  equations  of  motion  in  the  body  axis  for  yaw  and  roll. 

2.2.5  Stability  Derivative  Coefficients 

The  stability  derivative  coefficients  are  obtained  in  this  section  using  the  reference 
condition  =  &Q  =  wQ  =  vQ  =  qQ  =  rQ  =  0  .  This  is  a  typical  condition  for  a  sym- 


metrical  model  oscillating  about  its  zero-angle-of-attack  trim.  The  NACA  system11  will 
be  used  here  with  the  time  constant  defined  as  t*  =  1/ uQ  .  The  aerodynamic  moments 
are  defined  in  this  system  as 

"x 

=  PVCJA1C, 

(12) 

my 

=  pv’aic. 

(13) 

"z 

=  ^ v|AiC„  . 

(14) 

where  Vc 

=  1  (U2  +  vJ  +  wJ) 

"h  =  ci  • 

M*a  =  “wv« 


Mgt  =  Kx  , 


and  it  is  necessary  that  the  reference  length  ( l )  be  defined  in  half  lengths;  e. g. , 
l  -  d/2  ,  l  =  b/2  ,  etc. 


The  coefficient  form  of  the  derivative  (BMy/Bw)0  is  obtained,  using  Equation  (13), 
as  follows: 

^)o  =  2pVc(5)0AlC'0+^A,(lf)0-  (15) 

Since  Ovc/B»)  =  w/Vc  ,  the  first  term  on  the  right-hand  side  of  Equation  (15)  is 
negligible. 
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SECTION  3 

CAPTIVE  MODEL  TESTING  TECHNIQUE 


3.1  GENERAL  DESCRIPTION  OF  TECHNIQUE 

Methods  for  obtaining  rotary  derivatives  in  pitch,  yaw  or  roll  involve  restraint  of 
the  model  in  several  degrees  of  freedom  and  measurement  of  moments  required  to  sustain 
an  oscillation  for  specified  conditions  of  frequency  and  amplitude  (forced  oscillation 
technique)  or  measurement  of  model  motion  after  it  has  been  disturbed  (free  oscillation 
technique). 

The  basic  units  in  the  captive  model  testing  system  are  the  pivots,  the  support 
system  and,  in  the  forced  oscillation  technique,  mechanisms  for  forcing  the  model 
oscillation.  The  pivot  is  a  very  important  part  of  a  dynamic  stability  test  system, 
and  the  selection  of  the  type  of  pivot  and  design  will  depend  on  several  factors,  some 
of  which  have  a  direct  bearing  on  the  quality  of  data  that  may  be  obtained  with  the 
system.  Also  of  importance  is  the  suspension  system  and  the  influence  it  may  have  on 
the  aerodynamic  measurements.  Detailed  descriptions  of  model  pivots,  suspension  systems, 
and  balance  tare  damping  characteristics  are  presented  in  Sections  3.  1.  1  to  3.  1.3. 

3.  1.  1  Model  Pivots 

Pivots  of  various  types  which  are  employed  in  dynamic  stability  balances  include 
ball  bearings,  gas  bearings,  knif*  edges,  cones,  and  crossed  flexures.  The  selection 
of  a  particular  type  pivot  will  depend  on  the  system  operating  requirements  and  the 
environment  in  which  it  will  be  used.  One  of  the  most  important  factors  that  must  be 
considered  in  the  design  of  a  dynamic  stability  balance  for  aerodynamic  damping  measure¬ 
ments  is  the  damping  that  the  system  contributes  to  the  measured  wind  on  damping. 
Aerodynamic  damping  decreases  with  increased  Mach  number  and  becomes  a  small  percentage 
of  the  measured  damping  at  the  hypersonic  Mach  numbers.  One  exception  that  will  be 
discussed  in  Section  3.1.3  are  gas  bearings  which  generally  have  negligible  tare 
damping.  Because  of  its  importance,  tare  damping  will  usually  be  one  of  the  principal 
considerations  for  selection  of  the  pivot  type  or  bearing  for  roll  damping  balance 
systems. 

3.1. 1.1  Flexure 

High  pivot  friction  may  be  avoided  by  employing  crossed  flexures  of  the  general 
type  shown  in  Figure  4.  These  bend  in  such  a  way  that  the  center  of  rotation  remains 
at  essentially  a  fixed  position.  The  crossed  flexure  provides  stiffness,  which  is 
necessary  to  overcome  static  moments  on  models  tested  at  angles  of  attack  by  the  free- 
oscillation  technique.  Flexure  pivots  are  used  in  many  applications,  in  particular 
scientific  instruments,  because  they  have  many  favorable  characteristics  They  are 
simple  in  design,  yet  rugged  in  construction  for  their  flexibility;  they  do  not 
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introduce  hysteresis  or  backlash  in  the  system;  they  do  not  require  lubrication;  and 
they  are  not  susceptible  to  wear.  Their  main  disadvantage  is  that  they  are  not  suitable 
for  large  angles  of  rotation  because  of  stress  considerations,  especially  under  com¬ 
bined  loadings.  The  design  of  flexure  pivots  for  a  particular  application  requires 
consideration  of  the  model  loads,  moments,  and  oscillation  frequency  for  dynamical 
similitude.  Experience  at  the  AEDC-VKF  has  shown  that  crossed  flexures  have  linear 
stiffness  characteristics  and  structural  damping  that  is  repeatable  even  though  the 
flexure  is  of  multipiece  construction.  This  method  of  construction  is  desirable  because 
the  flexures  can  be  machines  as  a  single  unit  and  then  cut  into  sections  for  the  in¬ 
dividual  flexures.  The  theory  for  design  of  symmetrical  crossed-flexure  pivots  is 
given  in  References  13  through  16. 

Examples  of  flexure  pivots,  including  cruciform  and  torsion  tube  pivots,  that  have 
been  used  in  various  applications  are  shown  in  Figure  5.  For  each  of  these,  the  pivot 
tare  damping  varies  inversely  with  the  frequency  of  oscillation.  Strain  gages  mounted 
to  one  member  of  the  crossed  flexure  are  employed  to  obtain  a  signal  voltage  which  is 
proportional  to  the  angular  displacement. 

For  some  applications,  such  as  tests  in  very  small  tunnels  or  low  density  tunnels 
which  have  a  small  test  core,  it  may  be  impractical  to  machine  flexure  pivots  of  the 
type  shown  in  Figure  5.  The  AEDC-VKF  has  developed  a  small  dynamic  stability  balance 
using  a  commercially  available  instrument  flexure  pivot  of  the  type  shown  in  Figure  6. 
The  characteristics  of  this  pivot  were  found  to  be  repeatable  and  linear.  Angular 
motion  time  histories  are  measured  with  an  angular  transducer  described  in  Section 
3.2.  2.  1. 

3. 1.1.2  Knife  Edge  and  Cone 

Knife  edge  and  cone  pivots  of  the  type  shown  in  Figures  7  and  8  have  recently  been 
used  in  dynamic  stability  balance  systems  because  they  are  nearly  frictionless. 

The  principal  disadvantage  of  either  the  knife  edge  or  cone  pivot  is  that  some 
method  of  restraint  must  be  employed  to  hold  the  pivots  in  contact  and  fine  adjustments 
are  necessary  to  eliminate  binding. 

The  single-degree-of- freedom  knife  edge  pivot  shown  in  Figure  7  was  used  in  a  single- 
degree-of- freedom  free-oscillation  balance  system  developed  at  DTIIB.  Here  knife  edges 
are  used  to  constrain  the  model  and  also  form  the  axis  of  rotation. 

A  similar  technique  has  also  been  used  by  DTMB  for  a  three-degree-of-freedom  balance 
system;  however,  in  this  case,  the  model  mounting  system  is  a  cone  on  the  end  of  the 
sting  inserted  in  a  conical  depression  contact  within  the  model  (Fig  8)  The  model 
is  free  to  roll,  pitch,  and  yaw  about  the  cone  and  is  held  on  the  pivot  by  model  drag. 

In  each  of  these  cases,  model  attitude  is  measured  from  photographic  records  of  model 
motion. 

3113  Ball  Bearing 

Ball  bearing  pivots  may  be  used  in  a  dynamic  stability  balance  system  for  pitch  and 
yaw  derivative  measurements  if  high  resultant  model  loads  arr  to  be  supported  and  if 
aerodynamic  damping  is  large  so  that  bearing  damping  constitutes  a  small  part  of  the 
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measured  damping.  Experience  has  shown  that  a  ball  bearing  pivot  is  likely  to  have 
about  two  orders  of  magnitude  higher  damping  than  that  of  a  flexure  pivot  designed  to 
carry  comparable  loads.  Although  the  ball  bearing  pivot  has  the  advantage  of  being 
suitable  for  large  angles  of  rotation,  bearing  damping  is  a  function  of  the  amplitude 
of  oscillation  and  is  often  nonrepeatable.  A  ball  bearing  pivot  used  in  an  AEDC-VKF 
free-oscillation  balance  is  shown  in  Figure  9. 

Ball  bearing  pivots  are  quite  frequently  used  in  forced-oscillation  balances  when 
high  amplitudes  of  oscillation  cannot  be  obtained  with  flexure  pivots.  Although  the 
bearing  damping  is  high,  its  effects  may  be  eliminated  in  the  measurements  by  measuring 
the  forcing  torque  on  the  model  side  of  the  bearings,  as  shown  in  Figure  10.  For  such 
a  balance  system,  angular  contact  bearings  are  used  with  preload  to  eliminate  free  play 
in  the  bearings. 

Ball  bearings  are  used  in  Magnus  force  and  moment  balance  systems  where  models  are 
spun  to  high  rotational  speeds.  In  this  application,  bearing  friction  is  not  important 
since  the  measurements  are  static  forces  and  moments  in  the  yaw  plane;  however,  bearing 
heating  may  be  a  problem  at  extremely  high  rotational  speeds. 

Until  recently,  roll-damping  balance  systems  have  been  designed  with  ball  bearings; 
however,  new  systems  in  development  are  designed  with  gas  bearings  to  overcome  bearing 
damping,  wear,  requirements  for  lubrication,  and  sometimes  cooling.  Angular  motion 
time  histories  may  be  measured  with  the  type  of  angular  transducer  described  in  Section 
3.  2.  2.1. 

3.1.1. 4  Gas  Bearing 

Although  gas  film  lubrication  is  an  old  concept,  it  is  only  since  about  1950  that 
the  study  of  gas  bearings  has  noticeably  accelerated.  The  reason  for  this  is  that 
special  bearing  requirements  can  often  be  best  satisfied  with  a  gas  bearing.  Gas 
bearings  were  used  in  special  wind  tunnel  drag-force  balances  as  early  as  1956;  however, 
they  were  not  applied  to  balances  for  dynamic  stability  in  pitch  and  yaw  measurements 
until  about  1961  (Ref. 18),  and  roll  damping  until  1964  (Ref.  19).  Applications  in  wind 
tunnel  testing  techniques  have  been  extensive  since  these  early  developments.  This  is 
understandable  since  all  of  the  desirable  features  of  a  pivot,  which  include  high  load- 
carrying  capability,  low  damping,  and  unlimited  rotation,  are  combined  in  a  gas  bearing 
pivot;  whereas  flexures  and  ball  bearing  pivots  have  only  some  of  these  characteristics. 
Gas  bearing  pivots  are  used  extensively  for  dynamic  stability  testing  in  the  AEDC-VKF 
at  hypersonic  speeds  where  the  aerodynamic  damping  is  very  low.  It  is  difficult  to 
measure  damping  accurately  with  a  flexi.e  pivot  balance  under  this  condition  because 
of  the  high  percentage  of  tare  damping. 

An  example  of  a  journal-type  gas  bearing  pivot  used  in  an  AEDC-VKF  one-degree-of- 
freedom  dynamic  stability  balance  for  measuring  pitch  or  yaw  derivatives  is  shown  in 
Figure  11.  The  inner  or  core  section  of  the  bearing  is  supported  on  each  side  by  a 
sting,  and  the  model  is  mounted  to  a  pad  on  the  outer  movable  ring.  Supply  gas, 
generally  dry  nitrogen,  is  supplied  through  the  core  to  a  radial  set  of  orifices  in 
the  center  of  the  bearing  Filters  are  provided  in  the  core  to  remove  foreign  particles 
that  will  cause  bearing  fouling.  A  closed  gas  supply  and  return  system  is  not  used 
at  the  AEDC-VKF.  since  tests  have  shown  that  the  small  flow  rates  have  negligible 
influence  on  base  pressures  and  wakes.  Bearings  have  been  developed  at  the  AEDC-VKF 
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with  two  radial  rows  of  orifices  spaced  to  improve  the  yawing-moment  carrying 
capability.  Bearings  are  usually  designed  for  specific  applications;  however,  com¬ 
mercial  gas  bearings  now  exist  that  may  be  used  in  one-degree-of- freedom  balance 
systems.  This  is  an  important  fact  since  design,  development  and  fabrication  of  gas 
bearings  can  become  very  time-consuming  and  costly. 

The  gas  bearing  shown  in  Figure  il  was  designed  with  inherent  orifice  compensation 
according  to  the  procedure  outlined  in  Reference  20.  Analysis  of  this  type  of  bearing 
is  a  problem  of  compressible  flow  with  friction  in  three-dimensional  passages  that 
change  area  with  loading.  Treatment  of  this  problem  consisted  of  reducing  the  complex 
analysis  of  the  bearing  to  the  analysis  of  a  simple,  flat  plate  model.  As  will  be 
shown  later,  these  simplifying  assumptions  lead  to  a  model  realistic  enough  to  yield 
acceptable  results  for  design. 

The  gas  bearing  shown  in  Figure  11  differed  from  the  design  model  of  Reference  20 
in  the  following  ways:  (i)  the  gas  plenum  was  located  in  the  fixed  center  portion  of 
the  bearing  to  supply  gas  through  radial  orifices  to  float  the  outer  moving  ring  and 
(ii)  plates  were  added  to  the  ends  of  the  bearing.  These  provided  a  self-centering 
feature  required  for  proper  operation  and  some  degree  of  lateral  restraint.  In  addition, 
the  end  plates  were  designed  with  shields  to  direct  the  exhaust  gas  toward  the  center 
of  the  bearing;  otherwise,  gas  exhausting  radially  would  impinge  on  the  inside  of  the 
model  and  possibly  influence  the  damping  measurements.  The  bearing  shown  in  Figure  11 
is  capable  of  carrying  high  radial  loads;  however,  the  lateral  load  capability  is 
small,  yet  adequate  for  most  purposes.  Where  side  loads  are  large,  additional  gas 
supply  orifices  and  pads  should  be  provided  on  the  ends. 

The  performance  characteristics  necessary  for  evaluating  the  journal  gas  bearing 
for  application  to  dynamic  stability  testing  include  the  radial  load  capacity,  bearing 
pressure,  and  the  tare  damping.  Curves  illustrating  the  relationship  of  radial  load 
to  critical  pressure  are  presented  in  Figure  12.  The  critical  pressure  is  defined  as 
the  pressure  loading  for  an  outer  shell  eccentricity  ratio  of  approximately  0. 5  of  the 
bearing  clearance.  Permissible  loading  is  nearly  a  linear  function  of  the  supply 
pressure  and  agrees  with  the  performance  predicted  by  theory. 

As  noted  previously,  the  one  most  desirable  characteristic  of  the  gas  bearing  for 
use  in  dynamic  stability  balance  systems  is  the  extremely  low  tare  damping.  Data  on 
the  damping  of  the  gas  bearing  in  Figure  11  are  presented  in  Figure  13  as  a  function 
of  frequency  of  oscillation  and  radial  load.  In  true  viscous  damping,  the  damping 
moment,  ,  is  defined  as  a  moment  proportional  to  Mg  which  is  invariant  with 

oscillation  amplitude  and  frequency.  This  was  not  found  to  be  the  case  for  the  bearing 
when  evaluated  over  a  range  of  frequencies  and  loads  at  11.5  deg  oscillation  amplitude. 
The  results  of  the  evaluation  of  the  bearing  damping  presented  in  Figure  13  shows  an 
inverse  relationship  with  frequency  as  is  t’  e  case  for  flexures  and  a  dependency  on 
radial  load.  The  damping  is  in  the  form 

=  Co,b  , 

where  C  and  b  depend  on  radial  load.  In  actual  practice,  these  small  values  of 
tare  damping  are  usually  neglected.  For  example,  the  damping  shown  in  Figure  13  is 
only  0.8%  of  the  aerodynamic  damping  for  a  typical  blunt  cone  re-entry  configuration 
at  Mach  number  10. 
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Gas  bearings  have  been  used  extensively  in  industry  for  application  to  machinery 
operating  at  high  speeds  where  bearing  friction  is  very  important.  Satisfactory  opera¬ 
tions  at  high  rotational  speeds  with  very  low  friction  makes  the  gas  bearing  very 
attractive  for  use  in  roll-damping  balance  systems.  Until  recently,  ball  bearings 
were  used  exclusively;  however,  they  have  many  disadvantages,  such  as  significant  tare 
damping,  wear,  lubrication  requirements,  and  heating,  that  can  be  overcome  with  the  use 
of  gas  bearings.  A  typical  example  of  gas  bearings  applied  to  a  roll-damping  balance 
is  shown  in  Figure  14. 

Bearing  pivots  described  up  to  this  point  only  allow  rotation  about  one  axis.  Many 
problems  encountered  in  flight  can  only  be  investigated  experimentally  using  a  balance 
which  has  rotational  freedom  about  three  mutually  perpendicular  axes.  A  prime  example 
of  a  flight  phenomenon  of  this  type  would  be  the  coupling  of  the  rigid  body  pitch  and 
roll  modes.  A  three-degree-of- freedom  balance  provides  this  capability  and  in  addition 
simultaneously  performs  the  functions  of  the  several  one-degree-of- freedom  balances 
as  it  incorporates  into  one  balance  the  ability  to  measure  pitch,  yaw,  and  roll  damping. 
The  single-degree-of-freedom  balances  still  have  an  application  as  they  can  be  used 
to  restrain  the  motions  to  pure  pitch,  yaw,  or  roll  in  situations  where  this  becomes 
necessary. 

An  inherently  compensated,  spherical,  gas  journal  bearing  shown  in  Figure  15  has 
been  developed  at  the  AEDC-VKF  for  use  as  the  pivot  in  a  three-degree-of-freedom 
dynamic  stability  balance.  The  inner  core  of  this  bearing  consists  of  four  spherically 
surfaced  pads  arranged  in  a  manner  which  allows  the  bearing  to  resist  loads  in  any 
radial  direction.  Gas  is  supplied  through  an  orifice  located  at  the  center  of  each 
pad.  This  configuration  is  capable  of  supporting  maximum  radial  loads  of  100  lb  using 
nitrogen  or  160  lb  using  helium.  Changes  in  the  bearing  core  design  will  allow  the 
bearing  to  support  radial  loads  of  approximately  350  lb.  The  damping  inherent  in  the 
spherical  bearing,  like  that  of  the  cylindrical  bearings,  is  extremely  low. 

Continuous  motion  histories  provide  the  data  required  to  determine  the  aerodynamic 
derivatives.  For  very  small  flexures  or  bearing  pivots,  it  is  desirable  to  measure 
the  motion  without  physically  contacting  the  moving  parts  of  the  system.  A  device 
referred  to  as  an  angular  transducer  was  originally  developed  by  the  AEDC-VKF  for  use 
in  dynamic  stability  balance  systems21.  The  variable  reluctance  angular  transducer, 
shown  mounted  on  a  cylindrical  gas  bearing  in  Figure  11,  consists  of  two  E-cores 
mounted  on  the  outer  movable  ring  of  the  bearing.  When  the  coils  contained  in  the 
E-cores  are  excited,  magnetic  paths  are  established  through  the  E-cores,  their  adjacent 
air  gaps,  and  the  eccentric  ring.  As  the  bearing  rotates,  the  reluctance  of  the 
magnetic  paths  through  the  E-cores  and  eccentric  ring  remain  constant,  while  the  re¬ 
luctance  of  the  air  gap  changes.  This  change  in  reluctance  produces  an  analog  signal 
proportional  to  the  angular  displacement. 

Instrumentation  for  the  spherical  gas  bearing  consists  of  three  mutually  perpendicular 
angular  transducers  (Fig.  15).  Ibe  pitch  and  yaw  transducers  located  on  the  forward 
portion  of  the  instrument  ring  operate  from  an  eccentric  which  is  concentric  in  roll 
and  eccentric  in  pitch  and  yaw,  whereas  the  roll  transducer  located  in  the  aft  portion 
of  the  ring  operates  from  an  eccentric  which  is  concentric  in  pitch  and  yaw  and  eccentric 
in  roll. 
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Gas  bearing  technology  is  growing  very  rapidly,  and  Deny  texts  have  been  prepared 
on  the  methods  for  design  and  analysis  (see  the  Bibliography).  Reference  2  Vol.I, 
contains  detailed  summaries  of  design  methods,  material  selection,  and  measu.  iment 
techniques,  and  Vol.  II  is  devoted  to  applications. 

3.1.2  Suspension  System 

The  design  of  the  suspension  system  will  depend  for  the  most  part  on  the  type  of 
model  to  be  tested  and  the  measurements  required.  Models  may  be  generally  categorized 
as  (i)  two-dimensional  models  which  span  the  test  section  and  are  supported  from  the 
side  walls,  (ii)  three-dimensional,  semispan  models  which  are  mounted  on  a  reflection 
plane  supported  by  the  test  section  side  wall,  and  (iii)  three-dimensional,  full-span 
models  which  are  supported  by  a  sting,  a  side  strut,  or  a  magnetic  suspension  system. 
With  the  exception  of  magnetic  suspension,  each  model  mounting  system  presents  a 
potential  problem  when  being  used  for  tests  of  some  types  of  configurations,  that  being 
the  influence  of  the  support  on  the  aerodynamic  measurements.  When  this  problem  cannot 
be  dealt  with  in  captive  model  testing  by  proper  support  design  or  evaluation  of  the 
influence  of  the  suspension  system,  free-flight  testing  in  the  wind  tunnel  and  range, 
and  testing  with  magnetic  model  suspension  is  used. 

Support  systems  for  the  three  types  of  models  commonly  investigated  in  wind  tunnels 
may  be  generally  grouped  into  two  categories:  (i)  the  well  support  method  for  the  two- 
dimensional  and  three-dimensional  semispan  models  and  (ii)  the  sting,  strut,  and 
magnetic  suspension  system  support  for  the  three-dimensional  full-span  and  body-of- 
revolution  models. 

3. 1.2. 1  W 'all  Support 

The  wall  support  method  for  testing  half  models  has  the  advantage  that,  for  a  given 
test  section  size,  a  larger  model  can  be  tested  and  thus  a  larger  Reynolds  number  can 
be  obtained  than  with  a  complete  model.  The  large  half  model  also  simplifies  the 
model  construction,  instrumentation  installation,  and  routing  instrumentation  wires  and 
pressure  tubes  out  of  the  model  through  the  large  contact  areas  at  the  tunnel  walls. 
Large  angles  of  attack  can  be  obtained  with  wall -supported  models  without  introducing 
interferences  from  sting  supports.  In  addition,  measuring  systems  and  instrumentation 
for  dynamic  stability  tests  of  the  complete  model  or  control  surfaces  can  be  located 
outside  the  tunnel,  so  that  equipment  size  is  not  an  important  factor.  In  many  tests 
of  wings  or  control  surfaces,  this  is  the  only  method  of  support  that  can  be  used  to 
acquire  dynamic  stability  data  without  altering  the  mode!  contour  for  an  internal 
balance.  An  example  of  a  two-dimensional  and  three-dimensional  model  mounted  on  the 
tunnel  walls,  with  equipment  for  measuring  dynamic  stability,  is  shown  in  Figures  16 
and  17  to  illustrate  the  principal  advantages  of  the  wall- support  method  described23. 

When  models  are  mounted  directly  on  the  tunnel  walls,  the  effect  of  the  tunnel 
boundary- layer  flow  hac  to  be  considered,  since  it  may  change  the  flow  field  over  the 
model.  As  noted  by  van  der  Bliek24,  the  wall  boundary  layer  influenced  the  lift  and 
drag  coefficient  on  delta  wings  and  wing-body  combinations  mounted  on  the  side  walls 
and  had  a  tendency  to  promote  separation  and  stall  near  the  wi.^g  root.  Although  similar 
conclusions  have  not  been  formulated  for  dynamic  stability  tests,  the  influence  on 
static  tests  results  should  be  sufficient  to  initiate  concern  about  the  influence  of 
the  test  section  wall  boundary  layer. 
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In  order  to  alleviate  the  influence  of  the  wall  boundary  layer,  several  methods 
have  been  used  with  success  at  subsonic  speeds  by  various  investigators.  These  include 
boundary- 1 ayer  removal  through  porous  walls  or  slots  ahead  of  the  model  and  the  in¬ 
stallation  of  V-shaped  vortex  generators  on  the  tunnel  wall  ahead  of  the  attachment. 
Formation  of  large  disturbed  areas  on  the  model  close  to  the  wall  may  be  prevented  by 
using  fences  attached  to  the  model  or  displacing  the  model  from  the  tunnel  wall  by 
inserting  shims  between  the  model  and  walls.  Each  of  these  methods  has  the  disadvantage 
of  altering  the  shape  of  the  model  configuration  and  not  completely  eliminating  the 
influence  of  the  wall  boundary  layer. 

Tests  have  shown  that  wall  boundary- layer  effects  are  largely  eliminated  by  using 
reflection  planes  of  the  types  illustrated  in  Figure  18.  With  this  method,  a  reflection 
plane  is  introduced  as  the  longitudinal  plane  of  symmetry  and  is  displaced  from  the 
wall  to  allow  passage  of  the  wall  boundary  layer  flow.  For  some  mounting  methods, 
contouring  in  the  flow  passage  may  be  necessary  to  eliminate  choking  of  the  flow  and 
the  resultant  ineffectiveness  of  the  reflection  plane.  The  reflection  plane  can  be 
fixed  to  the  wall  and  the  model  rotated,  or  the  model  can  be  fixed  to  the  plate  and 
rotated  as  a  unit.  When  the  model  is  rotated  on  the  reflection  plate,  as  is  generally 
done  in  dynamic  stability  tests,  a  gap  must  be  allowed,  which  may  introduce  leakage 
around  the  model  reflection-plane  junction  and  change  the  local  flow  field.  Tests  of 
very  blunt  models  on  reflection  planes  will  produce  significant  boundary- layer  separa¬ 
tions  and  thus  alter  the  model  flow  field.  Information  on  reflection  plane  design  is 
given  in  Reference  26. 

3. 1.2.2  Sting  Support 

Sting  supports  are  principally  used  in  tests  of  complete  vehicle  configurations  and 
are  inserted  into  the  model  base  through  the  center  of  the  wake  to  minimize  their  in¬ 
fluence  on  the  model  flow  field.  Careful  attention  must  be  given  to  the  design  of 
the  sting  to  insure  that  it  provides  a  stiff  support;  however,  modification  of  the 
model  contour  may  be  necessary  to  accommodate  the  sting.  Such  a  modification  is  a 
compromise  between  a  small  diameter  sting  of  constant  length  for  a  distance  behind 
the  model  and  one  which  will  carry  the  aerodynamic  loads  without  significant  movement 
of  the  pivot  axis.  Often  the  consequences  of  the  compromise  are  not  clear,  and  there¬ 
fore  tests  are  required  to  evaluate  the  influence. 

Static  tests^7  have  shown  that  both  the  sting  diameter  and  sting  length  have  an 
influence  on  force  and  base  pressure  data  for  most  configurations.  Some  forced- 
oscillation  systems  have  shaker  motors  in  enlarged  sections  of  the  sting  support,  to 
simplify  the  overall  design  and  reduce  the  complexity  of  the  driving  linkages.  In 
these  cases  it  is  important  that  the  enlarged  section  be  located  an  ample  distance 
downstream  of  the  model  base  to  permit  proper  formation  of  the  wake.  Very  few  experi¬ 
ments  have  been  conducted  to  investigate  specifically  the  influence  of  the  sting 
diameter  and  length  on  dynamic  stability  measurements. 

A  series  of  tests  have  been  conducted  at  the  AEDC-VKF  (Ref. 28),  over  a  range  of 
Reynolds  numbers,  to  investigate  the  influence  of  the  sting  diameter  and  length  on 
the  dynamic  stability  of  a  10-deg  half-angle  cone  with  a  flat  base.  To  determine  the 
relative  effects  of  sting  length,  disturbances  downstream  of  the  model  base  were 
introduced  by  placing  a  20-deg  conical  windshield  at  various  stations  along  the  sting 
from  0.75d  to  3d.  Relative  interference  caused  by  the  sting  diameter  size  was  created 
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by  repeating  the  tests  with  diameters  of  0.4,  0.6,  and  0.8d.  The  test  results  included 
dynamic  stability,  static  stability,  and  base  pressure  data,  since  base  pressures  are 
very  sensitive  to  the  degree  of  interference  that  the  support  presents  to  the  base 
flow.  Representative  results  are  shown  in  Figures  19,  20,  21,  and  22  for  Mach  numbers 
2.5,  3,  and  4  over  a  range  of  Reynolds  numbers  and  sting  lengths  using  three  sting 
diameters.  Evidence  of  interference  caused  by  the  support  le  lgth  may  be  judged  by  the 
deviation  of  the  base  pressure  ratio  from  the  reference  curve  corresponding  to  the 
longest  sting  length  (ig/ d  =  3). 

Boundary- layer  transition  is  generally  at  the  model  base,  at  a  Reynolds  number  which 
is  slightly  less  than  the  value  where  the  minimum  base  pressure  occurs.  An  analysis 
of  the  results  based  on  the  base  pressure  data  and  schlieren  photographs  for  lg/d  =  3 
shows  that  the  Reynolds  number  range  was  sufficient  for  a  transitional  and  fully  turbulent 
wake.  There  is  evidence  in  the  data  presented  in  Figure  19  to  confirm  that  some  support 
geometry  configurations  can  have  an  influence  on  the  dynamic  stability  derivatives. 
However,  it  appears  that  if  the  support  is  designed  using  criteria  for  static  force 
and  base  pressure  tests27,  it  is  unlikely  that  the  dynamic  stability  data  will  be  in¬ 
fluenced.  This  observation  does  not  apply  to  configurations  which  have  contoured  bases, 
since  these  are  more  critical  with  regard  to  the  influence  of  the  support  because  the 
wake  condition  has  a  direct  influence  on  the  pressure  distribution  around  the  base.  A 
further  complication  arises  in  this  case  because  the  base  must  be  partially  removed  to 
allow  movement  of  the  model  over  the  sting.  In  cases  where  the  relative  influence  of 
the  support  system  on  damping  results  cannot  be  assessed  by  captive  model  testing,  wind- 
tunnel,  free-flight,  and  aeroballistic  range  tests  are  necessary  to  obtain  comparative 
results. 

Some  examples  of  a  very  rigid  sting  and  support  for  a  high  frequency,  free-oscillation, 
dynamic  stability  balance  system  are  shown  in  Figure  23  and  a  more  conventional  support 
for  a  low  frequency  system  is  shown  in  Figure  24. 

3. 1.2.3  Strut  Support 

Where  limitations  on  model  amplitudes  are  imposed  by  mechanical  interference  between 
the  model  base  and  sting,  these  may  be  overcome  by  mounting  the  model  on  a  transverse 
rod  support  system  so  that  it  is  free  to  rotate  on  the  rod  with  no  restraints.  An 
example  of  one  such  system  used  in  the  AEDC-VKF  50- in.  Mach  10  tunnel  is  shown  in  Figure 

25.  The  balance,  consisting  of  a  gas-bearing  pivot,  model  lock,  and  position  trans¬ 
ducer,  is  contained  within  the  model;  nitrogen  gas  for  the  bearing  and  the  model  posi¬ 
tion  lock  and  instrumentation  leads  pass  through  the  transverse  rods. 

Extension  of  the  transverse  rod  into  the  model  will  affect  the  flow  over  the  model 
primarily  in  the  region  aft  of  the  rod  and  along  the  sides  of  the  model.  A  region  of 
flow  separation  will  exist  on  the  model  ahead  of  the  rod.  Since  the  interference  flow 
field  caused  by  the  rod  will  depend  on  the  boundary- layer  flow  and  the  position  of  the 
interference  on  the  model  surface,  it  is  to  be  expected  that  measurements  obtained  with 
a  transverse  rod  support  system  will  be  dependent  upon  the  oscillation  amplitude  and 
the  Reynolds  number  of  the  flow.  An  example  of  such  a  dependence  is  shown  in  Figure 

26.  where  data  obtained  with  a  sting  support  are  included  for  comparison.  It  is  apparent 
from  these  results,  for  a  10-deg  half-angle  cone  at  Mach  number  10,  that  both  Reynolds 
number  and  oscillation  amplitude  are  important  parameters.  It  should  be  noted  that 

when  the  Reynolds  number  is  large  the  influence  of  Reynolds  number  and  amplitude  dis¬ 
appear  and  the  results  obtained  with  a  sting  support  and  transverse  rod  support  are 
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in  agreement.  However,  at  low  Reynolds  numbers,  very  large  differences  can  exist  in 
results  obtained  by  the  two  methods,  thus  indicating  significant  adverse  effects  from 
the  transverse  rod. 

The  influence  of  the  size  of  the  transverse  rod  supports  at  Mach  5  is  shown  in 
Figure  27  for  a  blunt  cone  cylinder  model  with  and  without  a  base  flare.  A  cylindrical 
rod  and  double-wedge  rod  were  tested  with  roughness  on  the  model  nose  to  produce 
turbulent  flow  and  thus  to  eliminate  separation  ahead  of  the  base  flare.  The  results 
show  no  significant  effects  caused  by  rod  size  for  diameters  up  to  30%  of  the  model 
diameter  and  no  effects  caused  by  rod  configuration,  i.e. ,  a  cylindrical  rod  as  com¬ 
pared  to  a  double-wedge  rod.  A  significant  reduction  of  the  damping  occurred  when 
roughness  and  the  transverse  rods  were  eliminated  for  the  configuration  with  the  flared 
base,  because  of  the  rod  wake  influence  on  the  region  of  separated  flow  ahead  of  the 
flare. 

Although  very  few  systematic  investigations  have  been  conducted  to  investigate  the 
effects  of  transverse  rod  geometry  on  dynamic  damping  over  a  range  of  model  configura¬ 
tions,  Mach  numbers,  Reynolds  numbers,  and  amplitudes,  the  data  available  tend  to  in¬ 
dicate  that  (i)  rod  influence  becomes  negligible  at  the  high  Reynolds  numbers  when  the 
flow  is  turbulent  and  (ii)  rod  size  effects  are  small  provided  the  rod  is  small  com¬ 
pared  to  the  model  size.  Tests  with  transverse  rod  supports  should  be  supplemented 
with  tests  employing  a  sting  support  balance  system  to  provide  comparative  data  at 
low  amplitudes,  where  the  influence  of  the  transverse  rod  is  most  pronounced. 

J.i.2.4  Magnetic  Model  Support 

Magnetic  suspension  provides  a  means  for  suspending  a  model  in  a  test  unit  without 
any  physical  support,  and  thereby  eliminates  any  questions  about  the  experimental 
results  being  influenced  by  the  method  of  physical  support  (sting  or  strut).  The 
magnetic  suspension  system  consists  of  a  series  of  electromagnets  which  produce  magnetic 
fields  and  gradients  of  magnetic  fields  to  provide  forces  and  moments  on  a  model.  The 
model  has  a  magnetic  moment  which  may  be  a  part  of  the  model  if  a  permanent  magnet  is 
used  in  the  model  construction.  Other  methods  for  providing  a  magnetic  moment  are 
described  in  Reference  29.  Included  in  the  system  for  control  is  an  optical  system 
which  monitors  the  position  of  the  model  by  light  beams  passing  through  the  test 
section  and  focused  onto  sensors.  When  the  system  is  in  equilibrium,  forces  and 
moments  caused  by  gravity,  aerodynamics,  and  inertial  effects,  for  the  case  of  a  moving 
model,  balance  out  the  applied  magnetic  forces  and  moments.  Changes  in  model  position 
activate  the  control  system,  and  currents  flow  through  the  magnetic  coils  to  produce 
forces  and  moments  on  the  model  to  hold  it  in  equilibrium. 

A  magnetic  suspension  system  with  V-magnet  orientation  developed  at  the  AEDC-VKF 
(Ref.  30)  for  static  force,  moment,  and  wake  measurements  is  shown  in  Figures  28  and 
29.  A  similar  system,  but  with  L  orientation  of  the  magnets  has  been  developed  by 
MIT  (Ref.  31)  to  obtain  static  force,  moment,  and  dynamic  stability  measurements. 

3.1.3  Balance  Tare  Damping 

3.1.3. 1  Description  of  Damping 

Wind  tunnel  models  have  to  be  mechanically  suspended  on  a  bearing  or  spring  to  pro¬ 
vide  oscillatory  motion.  Exceptions  to  this  are  recent  developments  of  magnetic  model 
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suspension  systems  which  have  been  refined  to  the  degree  that  damping  data  have  been 
obtained29.  Mechanical  suspension  introduces  damping  which  must  be  accounted  for  in 
the  analysis. 

Aerodynamic  damping  as  a  rule  cannot  be  measured  directly  with  most  balance  systems, 
since  it  is  a  part  of  the  measured  total  damping;  therefore  it  must  be  obtained  in¬ 
directly  as  the  difference  between  the  total  and  tare  damping.  With  gas  bearing  pivots 
this  is  generally  not  necessary,  since  they  usually  have  a  negligible  amount  of  tare 
damping.  It  is  customary  in  writing  the  linear  differential  equations  expressing  an 
oscillatory  motion  to  assume  that  the  total  damping,  which  is  the  aerodynamic  plus  tare 
damping,  is  proportional  to  6  .  Therefore  the  motion  in  pitch  for  a  free  oscillation 
may  be  expressed  by 

IY0  -  (M^a  +  M 0t)O  -  (Mga  +  M<9t)6>  =  0  , 

where  the  subscript  a  identifies  the  aerodynamic  term  and  t  identifies  the  tare 
term. 

The  term  tare  damping  is  applied  here  to  denote  the  operation  of  non- aerodynamic 
influences  which,  by  resisting  the  motion  and  absorbing  energy  from  the  system,  reduce 
the  amplitude  of  motion.  A  constant- amplitude  oscillation  would  continue  to  exist  if 
damping  forces  of  various  kinds  did  not  exist.  Damping  forces  arise  from  the  mechanical 
system  damping,  damping  capacity  of  the  materials  used  in  the  system,  and  air  damping. 
Although  it  is  convenient  to  classify  the  contributing  factors  in  this  manner  in  dis¬ 
cussing  tare  damping  as  applied  to  dynamic  stability  testing,  the  mechanical  system 
damping  and  damping  capacity  of  the  materials  are  inseparable  in  the  measurements  and 
constitute  the  term  M^t  .  Air  damping  is  important  only  insofar  as  determining  the 
magnitude  of  M^t  from  air-off  tests. 

Mechanical  system  damping  involves  distinguishable  parts  of  the  system  and  arises 
from  slip  and  other  boundary  shear  effects  at  mating  surfaces,  interfaces,  or  joints. 
Energy  dissipation  here  may  occur  as  the  result  of  dry  sliding  (Coulomb  friction) 
lubricated  sliding  (viscous  forces)  or  both. 

The  damping  capacity  of  materials,  or  material  damping,  refers  to  the  energy  dis¬ 
sipation  that  occurs  within  the  material  when  the  structure  is  undergoing  cyclic  stress 
or  strain.  These  internal  energy  losses  may  be  considered  as  being  due  to  imper¬ 
fections  in  the  elastic  properties  of  the  material.  In  a  truly  elastic  member,  not 
only  is  strain  proportional  to  stress  but  also,  on  removal  of  the  stress,  the  material 
retraces  the  stress-strain  relationship  that  existed  during  the  application  of  the 
stress.  Even  if  materials  are  stressed  within  the  elastic  limit,  some  energy  is  lost. 
Energy  dissipation  in  a  material  during  a  stress  cycle  leads  to  a  hysteresis  loop  in 
the  stress- strain  curve  and  the  area  o*  the  loop  denotes  the  energy  lost  per  cycle. 

Air  damping  results  from  a  loss  of  energy  caused  by  restrictive  effects  as  a  mass 
moves  through  a  medium  such  as  air.  The  damping  may  result  from  dynamic  pressure, 
viscous  effects,  and  vortex  formation  and  shedding;  however,  the  relative  effects  of 
each  are  not  known.  Theoretical  studies  of  the  viscous  damping  force  on  a  sphere  and 
cylinder  oscillating  in  a  fluid32,33  show  that  each  will  experience  forces  proportional 
to  the  velocity  and  square  root  of  the  fluid  density.  Experimental  air-damping  data 
for  circular  and  rectangular  plates,  cylinders,  and  spheres  are  presented  in  Reference 
34  and  show  the  effect  of  pressure,  vibratory  amplitude,  frequency,  shape,  and  surface 


23 


area  on  the  air  damping  for  models  that  are  mounted  on  the  end  of  a  cantilever  beam. 

The  results  show  that  air  damping  associated  with  plates  is  directly  proportional  to 
the  ambient  pressure  and  amplitude  is  indicative  of  the  influence  of  dynamic  pressure. 
Low  amplitude  results  for  cylinders  oriented  perpendicular  to  the  plane  of  oscillatory 
motion  are  in  agreement  with  a  theory  by  Stokes33  and  are  proportional  to  the  square 
root  of  the  density  but  show  no  variation  with  frequency. 

It  is  preferable  to  design  dynamic-stability  balance  systems  in  such  a  way  that  the 
tare  damping  is  small  compared  to  the  aerodynamic  damping.  For  some  classes  of  models 
at  hypersonic  and  hypervelocity  speeds,  where  the  aerodynamic  damping  is  low,  the  tare 
damping  can  approach,  or  even  exceed,  the  aerodynamic  damping.  An  example  of  the 
contribution  of  the  mechanical  system  and  air  damping  relative  to  aerodynamic  damping 
is  shown  in  Figure  30.  It  is  apparent  from  this  example  that,  as  the  Mach  number  is 
increased,  extreme  care  must  be  taken  to  obtain  measurements  of  tare  damping  with  a 
high  degree  of  accuracy,  so  that  the  difference  constitutes  an  accurate  measurement 
of  the  aerodynamic  damping. 

3. 1.3.2  Evaluation  of  Still-Air  Damping 

Mechanical  and  material  damping  can  be  determined  by  tests  in  still  air  at  atmos¬ 
pheric  pressure;  however,  this  can  introduce  large  errors  in  some  cases,  since  the 
damping  thus  obtained  includes  air  damping.  To  determine  the  magnitude  of  air  damping, 
tests  are  made  in  a  chamber  of  the  type  shown  in  Figure  31  to  permit  variation  in  the 
ambient  pressure.  The  wind  tunnel  can  also  be  used  for  this  test;  however,  its 
environment  must  be  free  of  even  small  amounts  of  airflow  and  turbulence.  A  typical 
set  of  tare  damping  measurements  for  a  blunt  cone  model  are  shown  in  Figure  32  for  a 
range  of  pressures  from  atmosphere  to  vacuum.  The  level  of  damping  at  a  vacuum  re¬ 
presents  the  mechanical  system  and  material  damping,  since  the  air  damping  is  essen¬ 
tially  zero.  An  example  of  the  magnitude  of  the  air  damping  as  a  percentage  of  the 
aerodynamic  damping  is  shown  in  Figure  33  for  a  series  of  model  shapes.  At  supersonic 
speeds  the  air  damping  is  on  the  order  of  2  to  4%,  but  increases  significantly  with 
Mach  number,  even  approaching  100%  of  the  aerodynamic  damping  at  Mach  10. 

The  necessity  for  accurate  tare  damping  results  has  been  mentioned  and  specific 
reference  made  to  some  of  the  sources  of  damping  in  a  mechanical  system,  one  being 
energy  dissipation  in  joints  of  multiple  piece  systems.  To  ensure  that  tare  damping 
is  repeatable  and  consistent,  ",  number  of  tare  damping  recordings  are  usually  made 
following  removal  and  replacement  of  the  model  and  flexure  or  pivot  system  until  tare 
damping  is  shown  to  be  repeatable. 

In  cases  where  the  axial  and/or  normal  forces  are  large,  measurements  are  made  with 
the  pivot  loaded  to  determine  its  influence  on  the  mechanical  damping  and  pivot 
stiffness.  In  most  cases  the  pivot  static  loading  will  have  little  effect  on  the  tare 
damping  characteristics. 

Since  the  model  undergoing  test  will  oscillate  at  a  frequency  different  from  that 
of  the  tare  damping  tests  at  a  wind-off  condition,  it  is  necessary  to  establish  how 
the  energy  loss  or  tare,  M^t  ,  at  a  vacuum,  will  vary  with  frequency.  For  flexure- 
type  pivots  Smith35  notes  that  the  tare  damping  per  cycle  is  a  constant  fraction  of 
the  energy  stored  in  the  spring  at  the  beginning  of  the  cycle,  independent  of  the 
frequency.  This  is  equivalent  to  =  C/u.  ,  which  was  found  earlier  in  experiments 
by  Welsh  and  Ward36. 
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3. 1.3.3  Damping  Characteristics  of  Different  Metals 

Damping  characteristics  of  materials  exposed  to  cyclic  stress  are  affected  by  several 
factors,  which  oust  be  considered  in  the  selection  of  materials  for  flexure  pivots. 

Some  of  the  factors  are 

(i)  condition  of  the  material  -  composition  and  effect  of  heat  treatment, 

(ii)  state  of  internal  stress  -  effect  of  machining  and  changes  in  stress  due  to 
cyclic  motion  and  temperature  histories, 

(iii)  stress  imposed  by  service  conditions  -  type  of  stress,  magnitude,  stress 
variations,  and  environmental  conditions. 

Factors  which  include  the  magnitude  of  the  stress,  stress  history,  and  frequency  may 
be  significant  at  one  operating  condition  and  unimportant  at  another.  Final  selection 
of  a  material  should  be  based  on  tests  of  a  specimen  in  the  configuration  in  which  the 
material  will  be  used  and  under  simulated  operating  conditions. 

It  should  be  recognized  that,  in  most  systems  where  oscillations  or  vibrations 
occur,  the  damping  introduced  by  the  properties  of  the  materials  employed  is,  in  general, 
only  one  factor  contributing  to  the  total  damping.  Nevertheless,  a  knowledge  of  the 
characteristics  of  structural  materials  is  Important  and  the  choice  of  material  should 
be  guided  by  its  damping  characteristics. 

Materials  having  a  Young’ s  modulus  which  ranged  from  6.2  x  106  for  magnesium  to 
28.5  x  106  lb/in.2  for  Armco  17.4  stainless  steel  have  been  tested  in  the  VKF  to  provide 
material  damping  data  for  use  as  a  guide  in  selecting  materials  for  flexure-pivot 
designs.  An  example  of  the  material  damping  results  for  a  range  of  frequencies  and 
beam  stiffnesses  is  shown  in  Figure  34.  Additional  information  on  material  damping  may 
be  found  in  Reference  37. 


3.2  DYNAMIC  STABILITY  DERIVATIVES  -  PITCH  OR  YAW 

The  basic  methods  for  obtaining  dynamic  stability  measurements  are  well  known,  dating 
back  to  early  1920  when  British  reports  on  the  subject  began  to  appear.  It  was  found 
that  aerodynamic  pitch  damping  derivatives  consisted  of  two  parts,  Mq  and  MA  ,  and 
experiments  were  devised  to  account  for  these  terms  separately.  These  experiments  jrere 
successful;  however,  it  was  concluded  that  the  technique  was  not  satisfactory  for  de¬ 
ducing  M^  at  high  speeds  (transonic)  or  high  frequencies.  Some  modes  of  oscillation 
such  as  short-period  damping  can  be  evaluated  satisfactorily  if  the  sum  of  the  dancing 
derivatives,  for  example  (CBq  +  CBi),  is  known,  even  though  the  separate  values  of 
CBq  and  Cb£  are  unknown.  Instances  may  exist  at  supersonic  speeds  in  which  th* 
application  of  test  data  to  calculations  may  require  the  damping  derivatives  in  some 
way  other  than  the  combined  form;  however,  the  general  practice  has  been  to  measure 
the  derivatives  in  combined  form.  Wind  tunnel  techniques  that  are  in  most  general  use 
at  supersonic  speeds  and  above  can  be  summarized  as  the  free-oscillation  and  forced- 
oscillation  methods 

The  free-oscillation  method  is  one  of  the  earliest  used  and  is  considered  to  be  the 
simplest  from  the  viewpoint  of  equipment  and  instrumentation.  The  model,  restrained 
by  a  pivot,  is  deflected  to  an  initial  amplitude  and  released,  and  the  decaying  motion 
is  observed. 
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With  the  forced-oscillation  technique,  a  restrained  model  is  forced  to  perform 
simple  harmonic  motion  by  means  of  an  oscillator,  transmitting  forces  through  a  linkage 
to  a  torque  member  which  is  rigidly  connected  to  the  model.  Normally  the  structural 
damping  of  balance  systems  used  at  supersonic  and  hypersonic  speeds  is  very  low  and 
the  torque  member  stiffness  is  also  low.  This  requires  the  system  to  be  operated  at 
or  near  its  resonant  frequency.  In  order  to  obtain  data  over  a  wide  frequency  range, 
the  model  moment  of  inertia  or  the  stiffness  of  the  structural  restraint  may  be  varied. 

An  alternative  method  that  may  be  used  for  tests  over  a  wide  range  of  oscillation 
frequencies  is  referred  to  as  inexorable  forcing;  and  as  the  name  implies,  the  model 
is  driven  in  a  controlled  manner  by  rigid  linkage. 

Most  forced-oscillation  systems  employ  special  instrumentation  techniques  which 
utilize  a  feedback  control  system  for  controlled  oscillation  of  the  model  when  negative 
damping  is  present.  In  addition,  these  control  systems  provided  amplitude  control  for 
measuring  positive  damping. 

Both  free-oscillation  and  forced-oscillation  balances  may  be  further  categorized  as 
either  small- amplitude  or  high- amplitude  systems.  A  small -amplitude  balance  measures 
local  values  of  damping  and  pitching  moment  slope,  since  the  model  oscillates  at 
amplitudes  of  about  one  or  two  degrees,  and  is  often  used  to  obtain  data  at  angles  of 
attack  other  than  zero.  The  balance  stiffness  or  restoring  moment  for  large-amplitude 
balances  is  normally  zero  or  near  zero.  Forced-oscillation  systems  may  be  operated  at 
about  ±15  degrees,  as  compared  to  amplitudes  which  are  only  limited  by  the  balance 
support  for  free-oscillation  balances. 

The  mathematical  treatment  of  systems  with  free-  or  forced-oscillation  motion  may 
be  found  in  many  texts;  however,  the  methods  are  treated  in  the  next  section  for  the 
sake  of  completeness. 

3. 2. 1  Theory 

The  theory  for  linear  and  non-linear  one-degree-of- freedom  and  three-degrees-of- 
freedom  free-oscillation  systems  will  be  treated.  The  forced-oscillation  theory  will 
be  developed,  followed  by  a  brief  treatment  of  local  and  effective  values  of  damping 
and  finally  the  aerodynamic  transfer  equations  will  be  developed. 

3.2. 1.1  Free  Oscillation  (One -Degree -Of -Freedom) 

The  equation  of  motion  for  a  body  in  pitch  was  derived  in  Section  2.2  and  has  the 
same  general  form  as  the  equations  of  motion  for  a  body  in  yaw  and  roll.  For  free- 
oscillation  systems  having  structural  stiffness,  the  equation  of  motion  has  the  form 

I,r;  -  (M<sa  ♦  «9t)e  («„,  ♦  »$t)9  =  0  .  (17) 

After  introducing  new  nomenclature  for  simplification  of  the  derivation  of  the  solution 
to  Equation  (17),  it  becomes 


+  cc*  ♦  K 6  - 


0 


(18) 
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From  Equation  (21)  for  the  natural  damped  frequency,  the  static  moment  is 

K  =  Iyo^  t  Iy (C/2Iy) 2  .  (24) 

The  term  I y (C/21 y ) 2  is  usually  very  small  and  can  be  neglected. 

As  shown  in  Equations  (17)  and  (18),  the  dynamic  damping  and  static  moments  consist 
of  moments  due  to  the  aerodynamics  and  moments  produced  by  tares  in  the  balance  system. 
The  tares  are  evaluated  from  wind-off  measurements,  and  the  final  aerodynamic  dynamic 
damping  and  static  moments  are  given  by 

C  =  -«*»  =  -  <«*.  + 

M5a  ”  ^  ,  ~ 

*  =  -  *9,  =  -  <«9a  + 

M0a  "  ^M9a  +  • 

where  the  subscripts  w  and  v  denote  wind-on  and  vacuum  conditions,  respectively. 

Since  the  structural  damping  moment  parameter  varies  inversely  as  the  frequency  of 
oscillation,  incorporation  of  Equations  (23)  and  (24)  along  with  the  equations  above 
gives 


M 


0a 


=  2Iy  loge  R 


(f/CjR),  -  (f/C  '  ,  ■— 


MPa  r  -  td) 


In  cases  where  the  balance  system  has  no  structural  stiffness  (systems  using  a  gas 
bearing  or  ball  bearing  pivot),  the  tare  damping  may  be  evaluated  independently  of  the 
model.  Semi- empirical  means  may  be  needed  to  find  the  tare  damping  at  the  frequency 
and  the  motion  axial  load  that  will  be  experienced  in  the  wind  tunnel.  Equation  (23) 
is  simply  used  to  yield 

=  2ITf  loge  R/C,r 

The  tare  damping  is  noted  as  Mat  and 

M^a  =  21 y f  loge  R'CjR  -  Mflt  . 

The  aerodynamic  stiffness  is  obtained  from  Equation  (24)  as 

Me.  = 

The  derivatives,  when  reduced  to  a  dimensionless  form  for  pitching  motion,  become 

caq  *  caa  =  "^(^'q^Ad2) 

C.o 
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It  should  be  noted  that  the  above  equations  were  developed  for  a  linear  model.  It 
will  be  shown  later,  in  Section  3. 2. 1.2,  that  erroneous  values  of  damping  may  result 
if  the  linear  equations  are  applied  to  models  which  have  a  nonlinear  restoring  moment. 
It  will  also  be  shown  in  Section  3. 2. 1.2  that  Equation  (22)  may  be  generalized  as 


where 


t2 


‘V°> 


" (C/21 y)t 

e  cos 


[T(t>  -<£]  . 


T(t)  =  Mt)  dt  . 


and  Equation  (23)  now  becomes 


C  =  -2Iyf 


loge  {K(u,2/u,x)*]/G 


yR 


3.2. 1.2  Effects  of  Nonlinearities 

Many  aerodynamic  configurations  of  interest  have  damping  rates  and  oscillation 
frequencies  that  vary  with  time.  Quite  often  these  parameters  are  functions  of  the 
dependent  variable  and  its  time  derivatives;  therefore  the  differential  equation  which 
describes  the  angular  motion  is  nonlinear.  In  many  cases,  these  nonlinearities  are 
slight,  and  the  amplitude-time  history  may  be  broken  into  intervals  over  which  the 
motion  is  approximately  exponentially  damped.  It  has  often  pegumed  that,  by 
applying  the  linear  theory  over  these  approximately  exponentially  damped  segments,  an 
accurate  estimate  of  the  aerodynamic  damping  may  be  obtained.  This  is  true  only  when 
the  aerodynamic  stiffness  is  constant.  Nicholaides38  has  shown  qualitatively  how 
errors  in  the  damping  measurements  can  be  made  using  this  approach.  In  order  to  show 
this,  Equation  (18)  will  be  used  in  the  form 

6  +  D,(?  +  =  0  (25; 

where 

Dj  =  t) 

D2  =  D2  (<9. t)  . 


Multiplying  Equation  (25)  through  by  2^/D,  and  integrating  yields 


'2  =  (t’J/D2) 


r*  ♦  r  id, 

i * t  jt.  • 


+  2D1D2](6)/D2)2  dt  . 


Assuming  that  tl  and  t2  are  times  that  peaks  in  the  motion  occur,  this  expression 
reduces  to 


a?  _  ^2 
-  1 


=  Jtt?  LD2  +  2DjD2]  (6/D2)2  dt  . 


i 
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A  dyrnmic  stability  condition  is  obtained  from  this  expression  as 

[i>2  +  2D1D2]  >  G  . 

For  zero  viscous  damping  (D1  =  0),  the  convergerce  or  divergence  of  the  motion  is 
dictated  by  the  rate  of  change  of  D2  with  respect  to  time.  Thus  the  variable  aero¬ 
dynamic  stiffness  term  D?  supplies  a  mechanism  by  which  energy  may  be  fed  into  or 
taken  away  from  the  system.  The  linear  theory,  Equation  (18),  assumes  that  D2  is 
constant  aad  therefore  that  the  constant  term  determines  the  convergence  or 
divergence  of  the  system.  One  can  readily  see  that  if  the  linear  theory  is  applied  to 
a  situation  where  D2  is  a  variable,  the  value  of  Dx  obtained  will  be  erroneous. 

This  is  true  even  if  the  linear  theory  is  applied  to  the  nonlinear  oscillation  over 
small  approximately  exponentially  damped  intervals.  The  magnitude  of  this  error  is, 
of  course,  dependent  of  the  magnitude  of  D2  . 

When  the  aerodynamic  nonlinearities  are  large,  Equation  (25)  is  very  difficult  to 
solve.  Fortunately,  as  mentioned  previously,  there  are  many  aerodynamic  configurations 
of  interest  for  which  these  nonlinearities  are  small.  In  addition,  these  shapes  usually 
have  damping  moments  that  are  small  compared  to  the  restoring  moment.  For  such  shapes 
Equation  (25)  may  be  written  as 

8  *  Dr  (8, 8,  t)8  +  [d20  +  D2i(8,8,t)]8  =  8  ,  (26) 

where 

°2  ■  D20  +  °21 
D20£  » 

D2  0  >>  D2i  • 

This  equation  may  be  reduced  to  the  van  der  Pol  equation,  which  is  given  as 

8  +  d20£  =  e$(8,8,  t)  ,  (27) 

where  $(8,8,  t)  contains  the  nonlinear  terms  and  6  is  a  small  dimensionless  parameter 
which  characterizes  how  close  the  system  is  to  a  linear  conservative  one.  For  the 
above  example, 

e$(8,8,  t)  -  -  \S),(8,8,t)8  t  D2l(8,8,t)8]  . 

The  motion  may  be  approximated  for  short  time  intervals  by  the  solution  for  a  harmonic 
oscillator  (e  =  0)  which  is  given  as 

8  -  80  cos  (v  (D20)t  +  <p)  . 

Since  the  right-hand  nde  of  Equation  (27)  is  small  but  not  zero,  8Q  and  4>  are 
actually  slowly  varying  functions  of  time. 

It  is  shown  in  Reference  39  that  a  second-order  nonlinear  differential  equation, 
like  Equation  (27),  may  be  reduced  to  an  equivalent  linear  second-order  equation  with 
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variable  coefficients  which  has  a  first-order  solution  identical  to  that  of  Equation 
(27).  This  equivalent  linear  equatior  .?j  be  written  as 

e-^e  +  [v(Dj0)  +  j>mVe  =  o . 

where  the  ip  subscripts  indicate  mean  values  over  a  cycle  of  oscillation. 

The  equation  above  can  be  rewritten  as 

0  +  +  C2(t )<9  =  0  .  (28) 

In  order  to  solve  this  equation,  the  following  transformation  will  be  used: 

d(t )  =  u(t)A(t)  .  (29) 

Differentiating  this  expression  and  substituting  into  Equation  (28)  yields 
uA  +  2uA  +  iiA  +  CjUA  +  CjUA  +  C2uA  =  0  . 


Choose  u  and  A  such  that 

2uA  +  CxuA  =  0  . 

This  reduces  the  previous  equation  to 

uA  +  iiA  +  CjiiA  +  C2uA  =  0  . 

Dividing  Equation  (30)  by  A  and  integrating  yields 

-i/c  dt 
u  =  Ke  1 


(30) 


(31) 


(32) 


where  K  is  the  constant  of  integration. 

Differentiating  this  expression  and  substituting  into  Equation  (31)  yields 


A  +  B(t)A  =  0  ,  (33) 

where 

B(t)  =  [c2  -  ic]  -  jcj  . 

According  to  a  theorem  by  Winter40,  the  asymptotic  solution  of  Equation  (33)  is 

A (t )  =  [B(t)J'i[K1  cos  Y(t)  +  K2  sin  Y(t)]  ,  (34) 

Y(t)  =  jr/B(t)  dt  . 


where 
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This  solution  applies  after  a  finite  length  of  time  has  elapsed,  such  that  the  oscilla¬ 
tory  motion  approaches  being  periodic. 

It  is  interesting  to  note  that  for  certain  forms  of  B't)  ,  Equation  (33)  may  be 
reduced  to  Bessel’s  differential  equation. 

Equations  (29),  (32),  and  (34)  may  be  combined  and  manipulated  algebraically  to 
yield 

i  ,  . 

4  -ijc.d)  dt 

e  cos  [T(t)  +  <p]  ,  (35) 

where  £(0)  and  B(0)  are  the  values  of  envelope  amplitude  and  C(t)  respectively 
at  t  =  0  . 

For  the  case  being  considered,  where  Ci  veries  slowly  and  is  small  compared  to 
C2  ,  B(t)  given  in  Equation  (33)  may  be  reduced  to 

B  =  (J  ^  C2  . 


In  order  to  use  this  solution  to  determine  C;  ,  cc  and  y  must  be  known  as  functions 
of  time. 

3.2.1  3  Free  Oscillation  (Three  Degrees  of  Freedom) 

Solutions  of  the  equations  of  motion  referenced  to  the  nonrolling  coordinates  will 
be  obtained  in  this  section4.  For  this  development  it  is  assumed  that  the  model  has 
triagonal  or  greater  symmetry  and  that  the  XYZ  axes  correspond  to  the  principal 
inertia  axes  of  the  model.  The  effects  of  slight  asymmetries  of  the  model  will  be 
considered,  i.e.,  asymmetries  which  have  a  negligible  effect  on  the  above  assumptions. 

3.2. 1.3.1  Equation  of  Pitching  and  Yawing  Motions 

For  a  model  of  the  above  type  rotating  about  a  point  fixed  in  space,  Equations  (5) 
describe  the  motion. 


These  methods  have  been  extracted  from  References  38  and  41 
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The  complex  angle  of  attack  may  be  reduced  to 

v  +  iw 


b  - 


~  $  +  ia  , 


where  it  is  assumed  that  and  d  ^  w/Vc  .  Substituting  the  time  derivatives 


of  this  complex  angle  into  the  equation  of  motion, 

1 1  -  HxPj  =  iflcoAtUCj,  +  iCn)  . 


(41) 


Assuming  that  first-order  linear  aerodynamics  are  adequate  and  that  the  motion  is 
described  by  a,/3,5t,p,P,  q,  and  r  ,  the  pitching-  and  yawing-moment  co¬ 
efficients  are  given  in  the  body  coordinates  XYZ  as 


The  terms  Cn0  and  Cn0  were  included  in  this  expansion  of  the  coefficients  to  account 
for  the  assumed  slight  asymmetry  of  the  model*.  As  outlined  by  Murphy41,  the  Magnus 
coefficients,  CQpi  ,  Cop/j  ,  Cnpq  ,  and  Cnpr  ,  have  a  negligible  effect  on  the  motion; 
therefore  they  were  not  included. 


For  a  missile  with  triagonal  or  greater  symmetry t, 


~^mcx 

=  -C.a 

C„<3 

=  -C.d 

=  -C.i 

C 

=  C 

z  c„ 

n  r 

■  q 

nq 

npa 

C«p  B 

cupa 

Substitute  these  relations  into  the  yawing-moment  equation,  then  multiply  the  result 
by  i  and  add  it  to  the  pitching-moment  equation  to  obtain 

i<c.  +  ic„)  =  i(cB0  ♦  icn0)  t  c„0(^  +  la)  + 

+  (cm  +  cniH'=  +  ii)  +  C„pa(p  +  ia)P  , 

^VC 


*  This  approach  wa  first  used  by  Nicolaides42  in  his  development  of  the  "Tricyclic”  theory, 
t  See  Inferences  38  and  41  for  a  discussion  of  syvetry. 
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where 

P  =  Pd/2VC  . 

Substituting  the  complex  angle  of  attack  referenced  to  the  body  coordinates, 

£  2:  P  +  ia  , 


into  the  above  relation  yields 

i(C„  +  1CB)  =  1(CB0  t  icn0)  +  c-e£  + 


+  (c 


C.a)  —  +  CMpof  P  . 


Kq  1  „ 


(42) 


It  now  becomes  necessary  to  transform  these  coefficients  from  the  body  coordinates 
to  the  nonrolling  coordinates.  This  transformation  is  given  as 


i(C,  +  icn)e14  =  1(C„  +  1C„)  . 

The  transformation  for  the  complex  angle  of  attack  is  given  as 

ff.1*  =  }  . 


This  is  easily  proved  by  expanding  this  relation  and  equating  real  and  imaginary 
part 8  as  follows 

A  +  iw\  v  +  iw 

- )  (cos  $  +  i  sin  $)  -  - 

V  Vc  ) 

V  =  V  cos  $  +  w  sin  $ 
w  =  w  cos  $  -  V  sin  $  . 

These  relations  are  the  nonrolling  velocity  components41. 

Murphy41  used  the  above  transformation  to  transform  Equation  (42)  to  the  nonrolling 
coordinates-  The  result  may  be  written  as 

i(CB  ♦  Cn)  -  i(CB0  ♦  iCn0)e  ^  +  CMa£  + 

+  ♦  CM. )  il  ♦  icMpa|p 

c 

Substituting  this  relation  into  Equation  (41)  results  in 

l  +  (G  ♦  iH)|  +  (J  +  iK)|  =  LelPt  . 


(43) 
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subscript  j  =  1  or  2  . 

Using  boundary  conditions,  the  remaining  constants  are  evaluated  as 

K  .  j  - 

’  \.l  ~Kl.l  +  ‘K,!  ' 

The  solution  given  in  Equation  (44)  is  'Tricyclic"  in  that  it  is  described  by  three 
two-dimensional  vectors  or  arms  (Fig. 35).  For  models  spiuning  at  a  constant  rate,  the 
first  two  rotating  arms  are  damped,  and  the  third  is  a  constant  amplitude  rotating 
vector  which  results  from  the  forcing  function  given  in  Equation  (43).  This  rotating 
trim  vector,  due  to  the  asymmetry  of  the  model,  causes  the  axis  of  symmetry  of  the  model 
to  rotate  on  an  imaginary  conical  surface.  As  P  increases,  the  magnitude  of  this 
rotating  trim  is  decreased.  For  a  statically  stable  model  with  sufficiently  small 
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damping,  the  total  angle  of  attack  may  increase  as  P  approaches  .  In  this  case, 
the  possibility  of  resonance  exists  only  with  the  nutational  frequency  K)  since  the 
precessional  frequency  (o2)  carries  a  sign  opposite  that  of  P  (Ref. 38).  This  will 
be  shown  later  in  the  text. 

The  restrictions  on  the  above  solution  are: 

(i)  Configurations  with  triagonal  or  greater  symmetry  having  only  slight  configura¬ 
tional  asymmetries. 

(ii)  First-order  linear  aerodynamics. 

(iii)  P  =  constant. 

(iv)  Steall  oscillation  amplitudes. 

If  the  relations  for  the  modal  damping  rates  and  frequencies  given  in  Equation  (44) 
for  j  =  1  and  2,  respectively,  are  first  added  together,  then  multiplied  together, 
the  following  exact  relationships  result41. 

\  +  K:  -  = 

+  ^  =  -»  = 

\\2  -  =  J  = 

V2  +  =  *  = 

For  most  configurations  of  interest 

(-H2  -  4J)  »  (G2  t  2iGH  -  4iK) I  . 


OopAd2 

2VLI 


^CMq  +  caa^ 


Ix 

JL  p 


QooAd 


dq,  Ad 

2V„I 


"Mpa 


Therefore  the  expression  for  the  modal  frequencies  and  damping  rates  is  approximated  by 


Ki +  ‘‘j  -  - 


(G  +  iH)  1 


•  {-  (H2  +  4J)}  ♦ 


G2  +  2iGH  -  4iK 


2i  (H‘  -  4J) 


In  order  to  derive  the  expression  for  the  gyroscopic  stability  factor,  only  pre¬ 
dominant  aerodynamic  effects  are  considered  in  this  relation,  which  reduces  it  to 

\)  *  iij  ~  *  [-iH  *  •  ^_(h2  *  4J)>]  • 


For  an  aerodynamically  statically  unstable  model,  J  <  0  .  It  is  obvious  then,  for 
this  case,  that  unless  -(H2  +  4J)  <  0  the  motion  will  diverge,  since  >0  If 
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H2  >  -4J  ,  then  Aj  =  0  ,  the  constant -amplitude  periodic  motion  results.  The  gyro¬ 
scopic  stability  factor  is  defined  as 


-h2  _  i  /ixy  p2 

4J  4q00Ad  \lj  CHa 


(45) 


The  stability  condition  is  given  as 

1/Sg  <  1  • 

Static  stability  is  insured  for  all  aerodynamically  statically  stable  models  (CMa  <  0) 
and  for  aerodynamically  statically  unstable  models  (CMa  >  0)  when  Sg  >  1  . 

If  we  let 

r  =  (1  -  1/Sg)-’  , 

then  t  <  i  for  CMa  <  0  ,  and  r  >  l  for  CMa  >  0  when  the  roll  rate  of  the  model 
is  sufficiently  high  to  insure  static  stability. 


(48) 


It  is  shown  in  Equation  (47)  that  the  Magnus  moment  tends  to  increase  the  damping 
of  one  modal  vector  and  decrease  the  damping  of  the  other.  A  model  will,  therefore, 
diverge  when 

I,CMClT/V  5  +  <W(1  ♦*)!  • 

For  a  spinning  model,  the  modal  frequencies  are  unequal,  as  evidenced  in  Equation 
(48).  The  largest  angular  rate  (^)  is  known  as  the  nutation  frequency  and  the  smaller 
angular  rate  (^2)  is  the  precession  frequency.  The  nutation  frequency  ('^ )  always 
carries  the  same  sign  as  the  roll  rate  (P).  therefore.  Kj  rotates  in  the  same  direction 
as  K3  .  For  an  aerodynamically  statically  stable  model  (r  <  l),  the  frequencies  are 
unequal  in  sign  and  therefore  cause  the  modal  vectors  Kj  and  K2  to  rotate  in  opposite 
directions  (Fig.  35)  When  P  =  0  ,  t  =  0  ,  then  ^  -  -x>2  ,  Kx  -  \2  .  and  the  motion 
is  described  by  lines,  ellipses  and  circles  (Fig. 36).  As  P  increases  (r  >  0)  the 
motion  becomes  flowerlike  and  finally,  for  large  P  the  motion  becomes  like  that  of 
a  fast  spinning  gyroscopic  pendulum  (Fig. 37).  For  an  aerodynamically  statically  unstable 
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model  (t  >  1),  the  frequencies  are  equal  in  sign  and  therefore  the  modal  vectors 
and  K2  rotate  in  the  same  direction.  Motions  characteristics  of  a  fast  spinning 
aerodynamical ly  statically  unstable  model  are  given  in  Figure  38. 

In  order  to  obtain  the  damping  rates  and  frequencies  of  the  modal  vectors,  the 
usual  approach  is  to  fit  the  solution  given  in  Equation  (44)  to  the  observed  motion, 
using  a  least  squares  differential  correction  method.  For  a  description  of  this 
method  and  its  application,  see  References  41  and  42. 

3.2  1.3.3  Solution  for  Varying  Holl  Hate 

In  many  cases  of  interest,  the  roll  rate  is  not  constant  as  assumed.  If  the  region 
of  resonance  is  avoided  and  the  roll  rate  is  slowly  varying,  the  solution  given  in 
Equation  (44)  may  be  used  to  fit  the  data  if  it  is  applied  over  intervals  of  the  motion 
where  P  is  nearly  constant. 

For  roll  rates  near  resonance,  P  is  usually  small  enough  such  that  G  »  H  and 
J  »  K  .  Roll  rates  varying  near  resonance  will  therefore  have  a  negligible  effect 
on  the  homogeneous  solution  of  Equation  (43).  This  is  not  true  of  the  particular 
solution  which  undergoes  drastic  variations  at  this  condition;  therefore  Equation  (44) 
is  not  valid  for  this  condition. 

In  order  to  investigate  the  problem  of  roll  rates  varying  near  resonance,  we  will 
rewrite  Equation  (43)  as 

l  +  (G  ♦  +  (J  +  iK)|  =  Lemt)  ,  (49) 


where 

Wt)  2:  f1  P(t)  dt 

*o 

G  +  iH  -  constant 
J  ♦  iK  ~  constant 

The  complementary  function  of  the  differential  equation  (49)  will  be  the  same  as 
that  given  in  Equation  (44)  The  task  is  now  to  obtain  a  particular  integral.  This 
will  be  accomplished  using  the  method  of  variation  of  parameters,  as  Murphy  did  in 
Reference  41 

Assume  the  particular  solution  is  of  the  form 

(X.  ♦i«Ot  (X,  *i*,)t 

*p  1  A,e  1  ♦  A,e  2  ‘  ,  (50) 

where  the  coefficients  Aj  are  variable  unknowns  and  \j  and  jo.  are  known  constants 
(j  1  or  2). 


39 


Applying  the  method  of  variation  of  parameters  results  in  the  following  two  simul¬ 
taneous  equations: 


(X .♦io,  )t  .  (X  ♦  icu  )t 
A:e  1  1  +  A2e  2  2  =0 


A1(X1  +  i^)e^  1  ^  +  A2(\2  +  io2) 


(X,  ♦  icujt 


=  Le 


i$>(t. 


(51) 

(52) 


Solving  Equations  (51)  and  (52)  simultaneously  results  in 

-(X,  ♦  ico.)t  ♦  i$(t) 

Le  1  1 


1  (\  -  +  i(^j  -  ^?) 

-(X2  ♦  lco2)t  ♦  i$(t) 


Le 


(^2  ~  \)  +  i  (^2  ~ 


These  expressions  are  integrated  and  substituted  into  Equation  (50).  The  solution 
to  Equation  (49)  is  then  given  as 


(X.  ♦  i^H  ^  (X2  ♦  i«2)t 

"  V2C 


S  =  Kje  ‘  1  +  K,e 


+ 


L  p  r<^i  •».,)(»-»,) 

Jo  L 


(X,  ♦  ioO(t  - 


i 


i  f*  i  P(t)  dt 
e 


In  order  to  use  this  solution.  P  must  be  known  as  a  function  of  time. 


By  arbitrarily  including  the  term  e  in  the  K3  vector,  the  solution  given  as 
Equation  (44)  can  be  made  to  approximately  compensate  for  the  changes  in  the  rolling 
trim  angle*.  This  solution  then  does  not  require  knowing  P  as  a  function  of  t  and 
is  given  as 


(X,*ioOt  (X,MWj)t  (X,*iu>)t 

1  1  ♦  K2e  2  2  +  K3e  3  3 


(53) 


where  ^3  =  P  . 

It  should  be  emphasized  that  the  solutions  obtained  thus  far  are  for  slowly  varylog 
roll  rates.  If  the  roll  rate  is  large  an*  rapidly  varying,  these  solutions  are  not 
applicable. 


Ibis  method  was  related  to  Mr  L.R. lard  in  a  private  communication  with  Professor  Robert 
S.Eikenberry  of  the  Department  of  Aerospace  Engineering.  University  of  Notre  Dame.  Notre  Dame. 
Indiana. 
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3,2. 1.3. it  Solution  for  Variable  Coefficients 

The  effects  of  variable  aerodynamic  stiffness  and  damping  were  discussed  previously 
in  Section  3. 2. 1.2.  In  that  section,  the  need  for  correcting  the  damping  measurements 
foi  the  effects  of  variable  aerodynamic  stiffness  was  emphasized.  It  was  shown  that, 
for  slowly  varying  coefficients  and  small  damping,  a  second-order  nonlinear  differential 
equation,  like  the  homogeneous  portion  of  Equation  (49),  can  be  reduced  to  a  linear 
equation  with  variable  coefficients  and  retain  the  same  first-order  solution. 

For  the  nonlinear  case  (J  =  J(£,£, t),  etc.),  the  homogeneous  portion  of  Equation  (49) 
may  be  written  correctly  to  first  order  as  follows: 

Z  i.  M 

f  +  t0(t)  +  i.l(t)]f  +  [j(t)  +  iK(t)]£  =  0  .  (54) 

Murphy4 1  derived  a  correction  to  the  linear  complementary  function  by  assuming  the 
following  solution  to  Equation  (54): 

~  i  <t>.  i<t>~ 

5  =  K,e  1  +  K2e  2  .  (55) 

where  Kj  is  a  variable  and 


<pj(t)  =  J *  ojj  (t)  dt  (J  =  1  or  2)  . 

This  solution  is  similar  in  form  to  linear  epicyclic  complementary  functions  obtained 
for  constant  coefficients. 

The  following  set  of  equations  was  obtained  by  substituting  the  first  and  second 
time  derivatives  of  the  solution  (55)  into  Equation  (54)  and  separating  real  and  imaginary 
parts: 


—  [Kj/Kj]  +  [kj /Kj ] 2  -aj  +  Otkj/Kjj  -  HOj  +  J  =  0  (56) 

dt 

+  2[kJ/KJ]wJ  +  Htkj/Kj]  +  Gcdj  +  K  =  0  .  (57) 

Since  cjj  »  kj  /Kj  and  kj/K-  is  slowly  varying,  Equation  (56)  may  be  reduced  and 
manipulated  algebraically  to  yield 

^ =  •  (58> 

This  is  essentially  the  same  as  the  frequency  relation  obtained  for  constant  coefficients. 
Equation  (57)  can  be  rearranged  to  obtain 

1  » 

H 


Kj/Kj  = 


u/j  n 

4<&j 


(59) 


where 
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K  _  1  (GoJj  +  K) 

J  2  (oJj  +  H/2) 

-  coj  +  H/2  . 

This  relation  shows  that  variable  aerodynamic  stifiness  and  variable  roll  rate 
both  affect  the  convergence  or  divergence  of  the  motion.  If  the  linear  theory  were 
applied  directly  to  motion  of  this  type,  the  complete  right-hand  side  of  Equation  (59) 
would  be  obtained  as  the  damping  rate. 

Equation  (59)  is  now  integrated  to  obtain 


Using  Equations  (58)  and  (59),  and  the  condition  of  slowly  varying  parameters,  the 
equation  above  can  be  rewritten  as 


^(0)  +  H(0)/2  *  X.t 

Ki  (t )  =  KAO)  -J - -  e  J 

J  J  Uj(t)  +  H(t)/2 


where  it  is  assumed  that 

V,  »  H/{4(^j  +  H/2)}  . 

Substituting  this  into  Equation  (55),  the  complementary  function  becomes 


i  =  Kj(0) 


^(0)  +  H(0)/2 


*  X.t  ♦  !<*>(t) 
e  1  + 


cxJj(t)  +  H(t)/2 

P^2(0)  +  H(0)/2 


+  K,(0) 


;w,(t)  +  H(t)/2 


La 


■'*  X-t  ♦  102  (t  ) 
e 


This  solution  may  be  applied  over  portions  ot  the  motion  history  where  is 
approximately  constant  if  ^j(t)  and  P(t)  are  known. 

3.2. 1.3.5  Graphic  Data  Reduction  Method  (Approximate) 

A  graphic  data  reduction  method  may  be  used  to  determine  the  modal  damping  rates 
(\j)  and  frequencies  (tuj).  In  order  to  apply  this  method,  it  is  necessary  to  convert 
the  tricyclic  motion  to  epicyclic  notion.  This  is  accomplished  using  the  roll  rate 
history  as  follows: 

~  -lPt  i(«.  -P)t  i(a>  -P)t  i$Q 

<e  =  Be  1  +  Ce  2  +  K3e  0  , 
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where 


B  =  Kte 


xit 


c  - 


K2e 


\2t 


i*0 

K3e  =  constant  . 

The  lotion  is  transformed  to  body  fixed  coordinates  by  rotating  each  angle  through  its 
corresponding  roll  orientation  (Fig. 39).  This  results  In  eplcycllc  lotion  about  the 
point  Kje1^0  . 

At  laxlaa  of  the  eplcycllc  motion,  B  and  C  add,  and  at  ilnlia  they  subtract 
(Fig. 40).  Vein®  the  maximum  and  ilnliun  amplitude  points  and  their  corresponding 
times,  the  variation  of  B  and  C  with  respect  to  tlie  Is  obtained  (Fig. 41).  Taking 
the  ratio  of  the  natural  logarlthis  of  these  expressions  yields 


\  = 


— - 10,.  (3h; 


(60) 


Hie  modal  frequencies  are  obtained  using  Figure  40  as  follows: 

277+7 


"i  = 


"*  =  “ 


(27 7-7) 
*2  -*1 


Csa  <  0  • 


Combining  Equations  (45).  (46).  and  (48)  with  the  above  frequency  relations,  the 
pltchlng-ioient  slope  Is  given  as 


*MaB 


p*  -  ^  -  cj}y 


Assuilng  the  roll  rate  Is  sufficiently  siall  that  Magnus  effects  can  be  neglected. 
Equations  (47)  and  (60)  can  be  combined  to  yield 


(c«q  +  Ca^« 


2  V.  1  1  /B,\ 

- 9 - 1°C«  I-4  ) 

M1  (1  +t)  (t,  -  t,)  *  \B,y 


tCy.  +  Cy4]g 


_2_  _ 1 _ 1__ 

Ad*  (1  -  r)  (t,  -  t,) 
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The  effect  of  balance  tare  damping  on  the  total  damping  rates  must  be  accounted  for 
before  the  aerodynamic  coefficients  are  obtained.  The  use  of  a  spherical  gas  bearing 
pivot  eliminates  this  difficulty  for  almost  all  cases  of  practical  interest,  due  to 
its  negligibly  small  tare  damping. 

3.2.1. 4  Forced  Oscillation 

Hie  motion  of  a  model  balance  system  which  is  rigidly  suspended  in  an  alrstreem  and 
forced  to  oscillate  in  one  degree  of  freedom  (pitch)  as  a  slnusodlal  function  of  time 
may  be  expressed  as 

+  "*>0  -  ("fla  +  uet)e  =  M  cos  ■  (61) 

wfcere  the  subscript  a  refers  to  the  aerodynamic  values  and  the  subscript  t  refers 
n  the  tare  contribution  of  the  model  balance  system 

Simplifying  the  nomenclature  by  letting 

-  <“9a  +  =  c  <32> 

and 

-  (“9.  +  ■«>  =  *  •  («3> 

Equation  (35)  becomes 

IT0  +  C8  +  K0  =  II  cos  cut  .  (64) 

The  complete  solution  of  this  equation  consists  of  the  sum  of  the  complementary  and 
particular  solutions.  It  can  be  shown  that,  when  M  =  0  ,  the  complementary  solution 
is  Equation  (19)  for  the  case  of  the  free-osclllation  system. 

To  solve  for  the  particular  solution  which  is  necessary  for  the  steady-state  case 
after  the  decay  of  the  initial  transients,  it  is  assumed  that  the  solution  is  of  the 
form 

8  -  8Q  cos  (cot  -  0)  , 

which  may  be  expanded  to  yield 

8  -  6^  (cos  cot  cos  0  +  sin  cot  sin  0)  . 

The  velocity  and  acceleration  are  then  found  as 

8  -  6Qco(-Bincot  cos  0  ♦  cos  cot  sin  0) 

8  -  -8Qco2  (cob  cot  cos  <p  +  sin  c^t  sin  8)  . 

Substituting  in  Equation  (64)  and  equating  coefficients  yields 

_ M _ 

(K  -  I^2)  cos  <P  +  Cco  sin  0  [ 
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tan  <P 


Ccj 

(1C  -  Iff) 


I 

I 


(65) 


Equation  (65)  shows  that  for  constan;  amplitude  motion,  when  the  forcing  frequency  is 
equal  to  the  natural  frequency  of  the  system  (<p  =  90  deg),  the  Inertia  term  balances 
the  restoring-moment  term  and  the  damping  moment  is  precisely  equal  to  the  forcing 
torque, 


C0ou  =  M  .  (66) 

Under  these  conditions,  it  is  only  necessary  to  measure  the  forcing  moment  II  ,  the 
angular  frequency  o>  and  the  system  amplitude  to  measure  the  damping  parameter  C. 

It  is  not  always  necessary  to  operate  at  resonance,  since  electronic  resolvers  can 
be  used  to  separate  the  moment  into  components  in-phase  and  90  deg  out-of-phase  with 
the  displacement.  Application  of  resolvers  is  covered  later  in  this  section. 

When  operating  at  resonance  the  displacement,  torque,  damping  and  frequency  are 
related  by 


c  =  u/eQu) . 

Introducing  the  relation j  for  C  and  K  (Equations  (62)  and  (63)),  the  aerodynamic 
damping  moment  and  the  equation  for  M9a  become 

*9.  =  +  «*).  -  ("«)»  <87) 

"9*  =  (“9.  +  «9t>,  -  («9t>Y  <88) 


where  v  and  w  indicate  the  vacuum  aad  wind-on  conditions,  respectively.  In  terms 
of  the  measured  quantities,  moment,  amplitude  and  frequency 


“*a 


M 


9a 


1 

% 


Iy 


The  derivatives  when  reduced  to  a  dimensionless  fora  become 

c.,  +C.a  =  H9»(»«/q»MJ) 

C.a  =  “9./«-8d  . 

A  further  examination  can  be  made  of  the  forced  oscillation  system  with  viscous 
damping  by  the  use  of  vector  diagrams  which  are  convenient  for  deriving  the  equations 
used  In  a  system  having  resolvers  to  permit  testing  near  a  resonance  condition. 
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The  equation  of  motion  with  viscous  damping  and  an  external  source  of  energy,  written 
earlier  in  the  form 

lYe  +  C0  +  K 6  =  M  cos  at  ,  (64) 

has  a  particular  solution  that  may  be  expressed  by 

6  -  &Q  cos  (cot  -  <p)  ,  (69) 

where  dQ  is  the  amplitude  of  the  steady-state  oscillation  and  0  the  phase  angle 
by  which  the  motion  lags  the  forcing  moment.  The  angular  velocity  and  acceleration 
functions  obtained  from  Equation  (69).  are  given  as 

w  -  -CO 0Q  sin  (cot  -  4>)  -  CO0Q  cos  (cot  -  0  +  v/2)  (70) 

d  =  -c cos  (at  -  0)  =  co260  cos  (at  -  0  +  v)  .  (71) 


Substli  .ti^  Equations  (69\  (7  ),  and  (71)  into  Equation  (64)  gives  the  following 

vect'*.  velatjun: 


'<  9 0  cc'  /  „  -  0)  -  C co6q  cos  (cot  -  0  +  v/2)  - 

-K0O  cos  (at  -  0)  +  M  cosat  =  0  . 

The  vector  relation  is  shown  in  Figure  42  for  forced  oscillation  with  viscous  damping. 
From  the  diagram  (Fig. 42),  we  obtain 


6  - 

0  ^[(K  -  IYa2)2  +  (Ca)2] 


tan  <t>  = 


Ca 


K  -  IjO2  ’ 


which  are  also  shown  as  Equations  (65)  in  the  earlier  derivation.  Referring  to  Figure 
42,  the  system  damping  is 


C  = 


where 


M  sin  0 

co  9„ 


sin  0  =  uin  (at  -  y) 


-  sin  at  cos  y  -  cos  at  sin  y 


and  hence 


C  =  (sin  at  cos  y  -  cos  at  sin  y) 
^0 


(72) 


46 


Tlie  in-phase  and  out-of -phase  moaents  and  angles,  indicated  by  Mr  and  M1  and  #r 
and  Si  .  respectively,  in  Figure  42.  are  related  by 

* 

sin  o)t  =  ICi/M  .  sin  y  = 

(73) 

cos  vt  -  iip/y  ,  cos  y  =  0r/#o  • 

Combining  Equations  (72)  and  (73)  leads  to  the  following  relationship  for  the  daaplng 
function: 

c  =  "lgr  -  “A 

o>(8\  +  e\)  ’ 


Introducing  the  relationship  for  C  ,  Equation  (72)  becoaes 


+  "4t 


“A  -  «A 

<*(e\  +  e\)  ' 


Since  the  structural  daaplng  aoaent  paraaeter  varies  inversely  with  the  frequency,  tho 
equation  for  the  aerodynaalc  danplng.  Equation  (67)  is  given  as 


Ml*r  -  Ml*r  -  “A 

+  e\)  w  vie*  ♦  e\)  Y  a,w  1 


where  v  and  w  Indicate  vacuua  and  wind-on  conditions. 

Proa  the  vector  representation  shown  in  Figure  42(a).  the  equation  for  the  static 
aoaent  say  be  expressed  as 


cos  0  =  cos  (o>t  -  y)  - 


(K  -  1^)6, 


K  -  l  jo2  +  —  cos  (cot  -  y)  . 


Ooablnlng  Equations  (73)  and  (74)  yields 


K  =  Io/  +  [V*  1  "fll  . 

1  le}+oi\ 

Using  Equations  (73)  and  (74).  the  aerodynealc  aoaent  derivative  becoaes 


In  a  case  where  the  phase  angle  between  the  reference  and  the  aodel  dlsplaceaent  is 
aero  (7=0,  figure  42),  the  resolution  of  the  forcing  aoaent  is  soaewhat  simplified. 
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Prom  the  phase  diagram  in  Figure  42(a),  we  have 


C<A  =  M  st=  0 

or 

C  =  («/Cj90)  sin  0  . 

and 

K 0Q  -  IY<*>2dQ  +  II  cos  0 


or 


K 


i*“ 


0  . 


The  aerodynamic  derivatives  now  become 


These  static  and  dynamic  derivatives  may  be  expressed  in  dimensionless  form  as 

C.Q  +  C.d  = 

c.a  =  Wl/fl-M)  • 

3.2. 1.5  Effective  and  Local  Values  of  Coefficients 

Methods  have  been  described  tor  obtaining  values  of  aerodynamic  damping  and  aero¬ 
dynamic  stiffness  from  models  in  wind  tunnels.  In  the  strict  sense  of  the  word,  all 
coefficients  that  are  obtained  from  these  measurements  should  be  classified  as  effective 
coefficients  (subscript  E);  however,  if  the  oscillation  amplitude  is  low  (6m  2:  l  deg), 
the  coefficients  may  be  termed  the  instantaneous  or  local  value  (subscript  L). 

Effective  values  of  the  coefficients  are  obtained  as  a  function  of  the  oscillation 
amplitude  6  as  the  model  oscillates  about  its  trim  attitude.  Local  values  are 
obtained  as  a  function  of  angle  of  attack  a  ,  and  the  data  are  obtained  with  the  model 
oscillating  at  small  amplitudes  (usually  by  means  of  a  forced-oscillation  balance 
system)  at  various  angles  of  attack.  There  appears,  however,  to  be  some  question  as 
to  the  validity  of  this  so-called  local  value  technique  for  obtaining  damping  data  and, 
at  this  writing,  the  question  has  not  been  satisfactorily  resolved. 

Redd  et  al.  * J  have  developed  a  method  which  relates  the  damping  derivatives  measured 
hy  the  two  methods  so  that  derivatives  can  be  expressed  as  a  function  of  instantaneous 
singula*  displacement  which  is  the  form  generally  used  in  trajectory  programs. 
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If  the  local  damping  coefficient  t(&)  can  be  expressed 

f  (d)  =  A  ♦  A/2  ♦  A^O"  ♦  A 696  ♦  Ae^e  +  . . . 

and  the  effective  damping  coefficient  f($0)  can  be  expressed  as 

f  (<90)  =  C  ♦  C2<9q  +  C„(9J  ♦  Cc£j  +  +  . . .  , 

the  coefficients  are  related  by 

C  =  A 
C2  =  (1/4)A2 

C4  =  (1/8)AU 

C6  =  (5/64)A6 

Ce  =  (7/128)Ae  . 

t 

From  these  equations.  It  can  be  shown  that  the  effective  damping- in-pltch  coefficient 
over  one  full  cycle  of  oscillation  designated  as  (C  +  cia^E  is  related  to  the  co¬ 
efficient  at  an  Instantaneous  angular  displacement  (?Bq  +  Cb6i)l  by 

<C.q  +  C.i>L  =  rS?  (A  +  A^J  +  +  A‘*‘  +  *>e' +  •  •  • ) 


Therefore.  If  values  of  (CBq  ♦  CB^)L  have  been  determined  in  small  amplitude 
forced-oscillation  tests,  the  amplitude  time  history  can  be  determined  by  using  the 
logarithmic  decrement  relation 


.  (C.^C.iHo^d^V.Hat/ZI) 

vr0Me 


Additional  details  on  the  method  and  the  derivation  of  the  series  coefficients  are 
given  In  Reference  43. 


3.2. 1.6  Transfer  Equations 

Aerodynamic  analysis  frequently  requires  the  transfer  of  derivatives  measured  at 
one  Mis  location  to  another  and  may  Involve  the  transfer  of  derivatives  along  the  Z 
am  well  as  the  X  axis.  Although  equations  for  transferring  results  along  the  X 
axis  are  given  In  several  texts  anc  reports  (Refs. 4.  10,  44  and  45).  few  references 
Include  equations  for  the  case  Involving  a  transfer  of  resuits  along  the  Z  axis  as 
well,  for  this  reason,  the  derivation  of  the  transfer  equations  are  presented  next 
to  Illustrate  the  procedure  for  the  case  of  longitudinal  stability  derivatives. 
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The  quantities  referred  to  the  initial  axes  location  do  not  have  subscripts,  and 
the  new  axis  location  is  designated  by  the  subscript  1. 


Sketch  (a)  shows  the  position  of  axes  XYZ  and  X1Y1Z1  .  The  transfer  distances  x 
and  t,  are  positive  in  the  senses  shown  in  the  sketch. 


The  equations  for  transferring  forces  and  moments  front  axes  XYZ  will  be  derived 
by  first  writing  the  equations  for  the  linear  velocity  components  and  the  angular 
velocities.  Considering  the  linear  velocity  components  of  .  they  become 


=  U  ♦  fe 
=  V  -  Pz  ♦  Rx  ► 


(75) 


Since  the  axes  the  data  are  being  transferred  to  are  parallel  to  the  Initial  axes,  the 
angular  velocities  are 


p,  =  p 

a,  =  Q 

R,  =  R  . 
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3.2. 1.6.  t  Stability  Derivatives 

The  transformation  equations  are  derived  here  by  first  expressing  the  force  and 
moment  component  derivatives  in  a  Taylor  series  expansion  and  employing  the  relation¬ 
ships  derived  for  linear  and  angular  velocities  and  forces  and  moments.  It  is  con¬ 
venient  here  to  introduce  the  standard  short  form  notation: 

X,Y,Z,  Components  of  the  external  forces 
parallel  to  axes  X  ,  Y  ,  Z  . 

L.li.N,  Moments  about  axes  X  ,  Y  .  Z  . 


3.2. 1.6.2  Longitudinal  Stability  Derivatives 

The  general  expansion  of  the  aerodynamic  force  components  along  the  body  XZ  axes 
and  the  pitching  moment  may  be  expressed  as  follows: 

X  =  x0  +  X„u  ♦  X-ii  +  X,q  *  X-q  ♦  X,w  +  Xji  +  ... 

Z  -  Z„  +  ZyU  +  Z;ii  +  Z^q  +  Z-q  +  Zt»  +  Zji  +  . . .  .  (78) 

M  =  K0  ♦  Muu  *  H-u  *  ♦  H-q  ♦  *  Up  +  . .  . 
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The  terms  XQ  ,  ZQ  ,  MQ  are  the  forces  and  moments  for  tin  reference  steady- flight 
condition. 

After  neglecting  all  derivatives  with  respect  to  acceleration  except  2^  and 
and  neglecting  the  derivative  Xq  ,  Equation  (78)  becomes 


X  =  X0  +  Xuu  +  Xww 


Z  =  Z0  +  Zuu  +  Zqq  +  Zww  +  Z^w 
M  =  M0  +  Muu  +  Mqq  +  M„w  +  M^w  . 


(79) 


If  we  next  consider  the  equation  for  the  force  along  the  body  X  axis,  the  de¬ 
rivatives  with  respect  to  velocities  u  and  w  for  station  1  are 


BXt  _  BFXl  _  BXt  Bu  Bxt  Bw 
xui  “  ^  ~  +  “  ' 


Bu  Bu  Buj  Bw  But 


Since  =  X  . 


Bu  Bw 

x-  =  +  ^ 


and,  using  Equations  (75)  and  (77), 


Xui  =  ^ 


(80) 


Similarly, 


Xwi  = 


BX,  BF„ 


Bxx  Bu  BXx  Bw 


Bwj  Bwj  Bu  Bwj  Bw  Bwt 


and 


Bu 


Bw 


and,  using  Equations  (75)  and  (77), 

By  following  the  same  procedure,  it  can  be  shown  that,  for  the  force  derivatives, 

=  ZU 
z*t  = 

=  Zi 

ZQ 1  =  z, +  V  -  V 


(81) 


(82) 
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Por  the  pitching-moment  derivatives,  neglecting  M^w  , 


"ui 


_  BMX1  _  3llj  Bu  Bllj  Bw  Blll  Bq 

Buj  BUj  Bu  Buj  Bw  But  Bq  But 


Using  Equations  (50), 


BM  Bu  BZ  Bu 
+  x 


Bu  Bu, 


BX  Bu 
BuBu,  Bu  Bu, 


BM  Bw  BZ  Bw  BX  Bw 

+ - +  x - z — - —  + 

Bw  Bu,  Bw  Bu,  Bw  Bu, 


i 


BM  Bq  BX  Bq 
+ - ♦  x  • 


Bq  But  Bq  Bux 


After  simplifying, 


Bu  Bu  Bu 


Bw 


Bw 


Bw 


Bu, 


Bui 


Bq  Bq 

+  Mn -  +  xZn - 

QBu:  Q  BUj 


Using  Equations  (77), 


“ui  =  “u  +  xZu  -  zXu  • 


(83) 


By  following  the  same  procedure,  it  can  be  shown  that  the  other  pitching-moment  derive 
tives  transfer  equations  are 


H.1  =  «,  ♦  XZ,  -  zX, 

*il  =  Mi+xZj-zXi 

V  =  "q  +  x(zq  +  “»>  -  z“u  +  +  *»>  +  • 

The  dimensional  form  of  the  derivatives  in  the  transfer  equations  are 

' 

^  =  -CaP.VA 

X*  =  -C^P.VA/2  , 

I 


(84) 


S3 


Zu  =  -  CjjPaVA  I 
\  =  -CNa^VA/2 
Zi  =  -CNiP„VA/4 

•  (85 

Z,  =  -CN0P„Vi./4 
Mu  =  C^VAd 
“w  =  C.^VAd/2 
“q  =  C^VAdJ/4 
“i  =  CB^„VAd2/4 

i 

After  the  equations  are  nondimensionalized,  Equations  <80)  through  (84)  for  trans¬ 
ferring  results  at  one  station  (XZ)  to  another  station  (XZ)1  become 


Following  the  same  procedure,  transfer  equations  for  lateral  stability  derivatives  may 
be  derived. 
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It  should  be  noted  that  the  signs  for  the  transfer  distances  from  the  initial  axes 
location  were  selected  as  positive  for  transfer  in  the  positive  sense  of  the  axes, 
i.e. .  transfer  distances  down  and  forward  are  positive.  References  are  not  consistent 
on  this  point,  as  may  be  seen  in  Reference  44  where  transfer  distances  aft  along  the 
X  axis  are  considered  positive. 

An  alternative  method  of  deriving  the  transfer  equations  for  the  case  of  transfer 
of  the  derivatives  (CBq  +  <W  along  the  X  axis  may  be  found  in  Reference  44.  Meyer 
and  Seredinsky45  point  out  that  the  transfer  equations  normally  found  in  the  current 
literature  are  invalid  at  moderate  to  large  angles  of  attack  and  where  large  transfc; 
distances  are  involved.  Transfer  equations  applicable  to  large-angle-of-attack  data 
may  be  found  in  Reference  45.  Wehrend44  conducted  an  investigation  at  Mach  0. 25  to 
2.2  to  test  the  method  of  computing  damping- in-pitch  derivatives  at  an  arbitrary 
location  from  derivatives  measured  at  two  locations.  A  similar  type  of  investigation 
was  conducted  by  Nicolaides  and  Eikenberry46.  These  investigations  showed  that  the 
transfer  equations  were  satisfactory  whether  the  configuration  was  stable  or  unstable. 

In  order  to  compute  the  damping- in-pitch  derivatives  at  another  body  station  from 
measurements  at  one  station,  the  transfer  equation  shows  that  not  only  the  damping 
derivatives  at  the  initial  station,  but  also  the  damping  normal-force  and  static 
characteristics,  must  be  known.  Usually  static  stability  data  are  known  from  static 
force  tests;  however,  damping  normal-force  characteristics  must  be  estimated.  This 
uncertainty  in  the  transferred  data  is  eliminated  by  measuring  damping- in-pitch  de¬ 
rivatives  at  two  body  stations,  which  eliminates  the  damping  normal  force  from  the 
transfer  equation.  When  this  is  done,  the  following  equation  is  obtained  using 
Equations  (86)  and  letting  z/d  be  zero: 


C.,  4  C.a 


d  (Cwq  +  C«a^2  d  (Ciq  +  Cia^  l 


-  2  ^  ^  C, 
d  d 


Na 


where 


The  reference  length  l  used  in  Reference  44  has  been  changed  to  d  to  be  consistent 
with  the  nomenclature  used  throughout  the  report. 

Damping  norm il- force  derivatives  may  also  be  computed  from  damping- in-pitch  measure¬ 
ments  at  two  body  stations  from  the  equation 
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CNq  +  CNa 


(Cbq  +  C»a^  2  “  ^Cmq  +  Cna^  l  "  ^  ^  ^  Ci 


■a 


*Na 


where 


C 


■a 


The  sul  rripts  2  and  1  on  the  damping  derivatives  refer  to  the  station  where 
the  derivatives  are  measured,  and  the  terms  xt  and  x2  refer  to  the  distance  between 
moment  centers.  The  constant  2  on  the  static  terms  in  the  equations  is  valid  only  when 
the  damping  derivatives  are  based  on  2VC  . 

3.2.2  Representative  Systems  for  Testa  of  Complete  Models 
3.2.2. 1  Free -Oscillation  Technique 

The  free-osci llation  testing  technique,  by  virtue  of  its  simplicity  because  coig>lex 
control  systems  are  not  required,  provides  the  experimenter  with  a  wide  latitude  in  the 
approach  to  balance  system  design. 

Although  the  merits  of  the  free-osci llation  technique  are  simplicity  in  construction 
and  use,  which  Implies  economy,  it  has  some  basic  deficiencies.  These  consist  of 
(1)  the  difficulty  in  obtaining  data  under  conditions  of  dynamic  instability,  since 
the  model  will  diverge  and  may  damage  the  crossed  flexures  if  they  are  amplitude  limited 
without  mechanical  stops,  (il)  a  reduction  in  model  angular  displacement  due  to  the 
static  restoring  moment  in  the  pitch  or  yaw  plane,  and  (ill)  the  inconvenience  of  re¬ 
ducing  data  Influenced  by  nonlinear  amplitude  and  frequency  effects. 

Different  techniques  have  been  used  for  displacing  models  mounted  on  a  free-osci llat ion 
balance  system  These  generally  consist  of  some  form  of  (a)  mechanical  release  from  a 
displaced  position,  (b)  application  of  a  mechanical  Impulse,  and  (c)  excitation  at 
resonance  and  interruption  of  the  forcing  source.  Methods  (a)  and  (b)  provide  positive 
control  over  the  initial  amplitude  but  have  the  disadvantage  of  requiring  large  forces 
to  obtain  high  initial  amplitudes  when  the  elastic  stiffness  of  the  flexure  is  large. 
Method  (c)  overcomes  this  difficulty  but  requires  more  complicated  equipment  to  force 
the  model  to  oscillate  at  resonance.  Aside  from  using  electrical  excitation  for  re¬ 
sonance,  excitation  can  be  provided  by  using  alternating  air  jets  impinging  on  the 
inside  of  the  model. 


i 
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Many  free-osclllation  dynamic  stability  balance  systems  have  evolved  in  the  United 
States  and  in  NATO  countries  to  obtain  information  experimentally  for  support  of  vehicle 
development  and  theoretical  studies.  To  show  the  basic  features  of  most  of  the  systems 
in  current  use,  examples  of  some  systems  will  be  described  which  are  designed  with 
crossed- flexure  pivots,  and  gas  bearing  pivots. 

To  perform  dynamic  stability  tests  at  Mach  numbers  from  18  to  21  in  the  AEDC-VKF 
100  in.  diameter  hypervelocity  tunnel  (Tunnel  F),  which  has  a  useful  run  duration  of 
about  msec,  a  high  frequency,  free-oscillatlon  balance  system  is  used.  To  obtain  up 
to  about  15  cycles  of  data  during  a  run,  high  model  frequencies,  minimum  inertia  models, 
and  high  crossed- flexure  stiffness  are  required.  A  flexure  system  which  provides  a 
model  frequency  of  150  c/s  at  oscillation  amplitudes  up  to  1  3  degrees  is  shown  in 
Figure  43  (Ref. 47). 

To  displace  the  pivot  to  an  initial  amplitude  a  couple  rather  than  a  moment  is  used 
to  minimize  a  beat  oscillation  that  would  be  imposed  on  the  oscillation  envelope  and 
reduce  the  total  stress  experienced  by  the  pivot. 

In  operation,  the  model  and  pivot  system  is  displaced  to  a  desired  amplitude  by  the 
pneumatic  operation  of  two  displacement  arms  which  apply  a  couple  about  the  pivot  axis. 
Once  the  des'red  amplitude  is  reached,  two  cams  operated  by  the  releasing  piston  make 
contact  with  the  displacement  arms  and  instantaneously  release  the  applied  couple. 

A  low-frequency,  free-oscillation,  flexure  pivot  balance  system  also  in  use  in  the 
AEDC-VKF  is  shown  in  Figure  44  to  illustrate  the  simplicity  of  the  design  and  a  method 
for  displacing  the  model.  Two  jets  located  180  deg  apart  are  mounted  to  the  sting 
just  inside  of  the  model  base  and  are  alternately  supplied  with  high  pressure  air 
through  an  oscillating  servo  valve. 

Forces  generated  hy  lmplngment  of  th u  Jets  on  the  inside  of  the  model  cause  the 
model  to  oscillate.  Amplitude  is  controlled  by  the  jet  pressure,  and  frequency  is 
controlled  by  the  servo  valve.  The  model  may  be  forced  to  oscillate  at  its  natural 
frequency  and  when  the  required  amplitude  is  achieved  the  air  excitation  is  terminated 
and  the  model  is  free  to  oscillate. 

A  free-oscillatlon  flexure  pivot  system  developed  by  ONERA45  is  shown  in  Figure  45 
to  illustrate  a  mechanical  system  for  displacing  the  model.  For  most  flexure  pivot 
balance  systems,  strain  gages  are  bonded  to  the  cross  flexures  to  provide  a  signal 
output  for  recording  the  amplitude  time  history  and  model  oscillation  frequency. 

Free-osclllation  balance  systems,  having  either  ball  bearings  or  gas  bearings  for 
the  pivot,  nay  not  have  mechanical  stiffness  to  provide  a  restoring  moment.  For  these 
systems  the  model  is  displaced  by  a  mechanical  system  and  released.  If  the  tunnel  has 
a  model  Injection  system,  the  model  can  be  preset  to  the  required  amplitude,  locked, 
and  then  released  after  injection  into  the  airstream.  An  example  of  an  AEDC-VKF  free- 
oscillatlon  balance  system  used  in  this  mode  is  shown  in  Flgur?  46. 

Most  of  the  free-oscillation  balance  systems  described  thus  far  have  been  developed 
for  use  in  relatively  small  wind  tunnels.  A  free-oscillation  balance  developed  by 
the  AEDC-PWT  (Ref .  49)  for  tests  in  the  16  ft  transonic  and  supersonic  tunnels  is  shown 
in  Figure  49.  Hie  balance  has  a  single  vertical  flexure  in  the  center  and  two  longi¬ 
tudinal  flexures  at  the  sides  for  a  pivot.  A  hydraulic  mechanism  is  provided  for 
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initially  displacing  the  flexure  and  model,  and  a  hydraulic  cylinder  provides  a  means 
for  locking  the  system.  Strain  gages  on  the  flexures  are  used  to  obtain  a  measure  of 
the  amplitude  decay  history. 

A  free-oscillation  dynamic  balance  system  has  been  developed  by  the  University  of 
California  Institute  of  Engineering  Research  to  measure  dynamic  stability  derivatives 
in  a  low  density  wind  tunnel50.  Figure  50  shows  the  components  of  the  balance  system. 

The  model  is  mounted  on  a  transverse  rod  supported  in  the  vertical  plane  by  bearings. 

The  bearings  have  an  electrically  vibrated  inner  bearing  race  to  minimize  the  friction 
and  provide  repeatable  bearing  tares.  System  restoring  moment  is  provided  by  a  quartz 
fiber  attached  to  the  tip  of  the  transverse  rod  support.  A  cross  shaft  is  mounted  on 
the  end  of  the  transverse  rod  support  to  displace  the  model  by  the  use  of  an  electro* 
magnet  and  also  interrupt  a  beam  to  a  photocell.  One  arm  casts  a  shadow  on  a  screen 
to  provide  direct  oscillation  angle  measurements. 

A  small  scale,  free-oscillation,  dynamic  balance  developed  by  the  AEDC-VKF  for  tests 
in  a  small  low  density  tunnel  is  shown  in  Figures  51  and  52.  The  commercial  flexure 
pivot  is  0. 125  in.  in  diameter.  The  fixed  portion  of  the  pivot  is  supported  by  two 
supports,  and  the  movable  center  portion  of  the  flexure  is  clamped  in  the  balance 
centerpiece.  As  the  movable  member  is  rotated,  the  air  gap  between  it  and  the  E-core 
changes  to  produce  a  change  in  signal  level  proportional  to  the  angle. 

Angular  position  histories  for  the  gas  bearing  pivots  are  recorded  with  the  E-core 
transducer  described  in  Section  3. 2. 2.  1.  Aerodynamic  damping  from  a  free-oscillation 
system  is  obtained  through  an  analysis  of  the  record  of  motion  recorded  on  magnetic 
tape,  film,  or  oscillograph  recordings  to  obtain  the  logarithmic  decrement.  Reduction 
of  the  data  from  a  magnetic  tape  may  be  accomplished  with  a  high  speed  computer  to 
curve  fit  the  motion  data,  read  the  peak  amplitudes,  and  finally  curve  fit  data  in  the 
form  of  log  of  the  amplitude  versus  cycles  of  oscillation.  If  the  data  are  recorded 
on  film  or  an  oscillograph,  many  manual  operations  are  involved  which  are  tedious  and 
reduce  the  accuracy  of  the  data  reduction. 

Hie  ablation  of  a  vehicle' s  thermal  protection  system  can  result  in  motions  which 
have  a  significant  effect  on  the  vehicle  performance.  Study  of  these  effects  experi¬ 
mentally  has  been  somewhat  restricted  because  present  capabilities  of  wind  tunnel 
facilities  do  not  provide  a  complete  duplication  of  re-entry  flow  conditions.  Ablation 
affects  dynamic  stability  primarily  by  changes  in  body  shape,  mass  distribution,  and 
the  mixing  of  gaseous  products  of  ablation  with  the  boundary  layer.  At  the  present 
time,  study  of  the  effects  of  ablation  have  proceeded  along  the  lines  of  using  low 
temperature  ablators  and  porous  model  blowing  techniques.  This  is  desirable  in  a 
study  of  shape  change  and  its  effects  on  stability;  however,  the  results  have  limited 
application  for  important  parameters  such  as  material  properties  and  phase  angle. 

Mass  addition,  on  the  other  hand,  offers  possibilities  for  simulating  many  of  the 
important  parameters  and  enables  measurement  and  systematic  variation  of  these  parameters. 
Some  of  these  parameters  include  mass  injection  rate,  distribution,  phase,  wave  form, 
gas  properties,  and  frequency  of  oscillation. 

LTV  has  pioneered  the  development  of  techniques  and  equipment  to  Investigate  the 
effects  that  mass  addition,  through  porous  models,  has  on  dynamic  stability.  T%ese 
investigations  have  been  conducted  in  the  LTV  hypervelocity  wind  tunnel,  which  is  a 
modified  hotshot-type  facility  with  relatively  long  run  times  (250  milliseconds  at 
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Mach  17)  at  near  constant  flow  conditions.  Dynamic  stability  test*,  are  conducted  with 
a  high-frequency  flexure  pivot  balance  to  provide  a  sufficient  number  of  cycles  of 
oscillation  for  adequately  determining  the  aerodynamic  damping. 

The  porous  model  and  balance  developed  in  the  investigation  are  shown  in  Figures  53 
and  54.  In  addition  to  supporting  and  oscillating  the  model  at  the  desired  frequency, 
the  flexure  must  provide  a  means  of  transferring  injection  gas  from  the  fixed  sting  to 
the  oscillating  model.  The  flexure  pivot  designs  shown  in  Figure  54  are  well  suited 
for  these  purposes  since  they  provided  a  means  for  bringing  gas  to  the  model  through 
passageways  in  the  flexures  without  increasing  the  tare  damping  of  the  system.  Relative 
to  the  rectangular  flexure  cross-section  design,  the  cruciform  design  has  greater  stiff¬ 
ness  in  the  axial  and  nonal  components;  however,  the  rectangular  design  has  the  advan¬ 
tages  of  easier  fabrication  and  longer  supply  passages  for  a  given  spring  constant. 

The  mass  injection  balance  was  fabricated  from  one  piece  of  aluminum,  and  the  flexures 
were  designed  to  operate  at  very  low  stresses,  thus  ensuring  low  tare  damping.  Model 
motion  histories  are  recorded  with  outputs  from  strain  gages  bonded  to  the  balance 
flexures. 

Details  of  the  gas  supply  system,  test  procedures,  and  typical  test  results  are 
presented  in  References  51  and  52. 

Some  of  the  devices  that  have  been  used  in  the  acquisition  and  reduction  of  data  to 
damping  coefficient  form  are  described  next. 


3.2.2. 1.1  Daatpometer 

A  coserclal  electronic  apparatus  which  nay  be  used  to  evaluate  automatically  the 
logarithmic  decrement  and  frequency  of  decaying  oscillations  is  the  damponeter  which 
waa  developed  by  Olsson  and  Orlik-Ruckemann53  (Fig. 55).  The  design  of  the  instrument 
is  baaed  on  the  fact  that  a  damped  oscillation  nay  be  represented  by  a  rotating  vector 
on  the  screen  on  a  cathode- ray  tube;  the  length  of  the  vector  is  proportional  to  the 
oscillation  amplitude  and  is  thus  a  measure  of  damping.  The  curve  described  by  the 
motion  of  the  end  of  the  vector  is  a  logarithmic  spiral  illustrated  in  Figure  56. 

For  a  decaying  harmonic  oscillation,  the  input  voltage  to  the  damponeter  nay  be 
written  in  the  manner  of  Equation  (22)  as 

- (C/2 I T ) t 

El  -  E0e  cos  o>dt 

and  the  ordinate  of  the  envelope  of  the  harmonic  oscillation  is 


where  voltage  Ex  is  proportional  to  the  length  of  the  rotating  radius  vector  on  the 
cathode- ray  screen. 

The  screen  of  the  cathode-ray  tube  is  covered  by  a  circular  disc,  as  shown  in  Figure 
56.  which  has  a  series  of  radial  slots  equally  spaced.  As  the  spot  representing  the 
end  of  the  vector  moves  in  a  spiral  path  around  the  screen,  it  passes  slots  in  the 
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disc  which  allow  light  to  pass  to  a  photocell.  The  photocell  triggers  an  electronic 
counter  (damping  counter)  which  gives  the  total  number  of  slot  passages,  n  ,  the 
vector  passes  in  going  from  the  outer  radius  rl  to  the  inner  radius  r2  ,  which 
contains  s  number  of  radial  slots. 


Denoting  the  time  at  which  the  spot  enters  the  slots 
by  t2  ,  the  expression  for  the  loarithmic  decrement  A 


tj  and  the  tine  it  leaves 
can  be  derived: 
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and 


A  =  ttC/I^  =  s/n  loge  (rt/r2)  . 


The  frequency  of  oscillation  is  measured  with  a  second  electronic  counter  by  counting 
the  number  of  cycles  of  a  crystal  oscillator  during  the  period  in  which  the  damping 
counter  is  in  operation  from  a  preset  number  nt  to  another  preset  number  n2  .  If 
the  frequency  of  the  crystal  oscillator  is  fc  and  the  reading  of  the  counter  nc  , 


Different  discs  having  a  range  of  radius  ratios  and  number  of  slots  are  used  to 
obtain  good  resolution  for  A  and  f  .  As  noted  in  Reference  53,  the  dampom%ter  can 
be  used  to  obtain  an  accuracy  of  <  2%  for  logarithmic  decrements  smaller  than  2; 
for  higher  values  of  the  decrement,  the  accuracy  decreases. 

3.2.2. 1.2  Travel  Summation  Method  of 

Free-Oscillation  Data  Analysis 

There  may  be  instances  when  performing  wind  tunnel  dynamic  stability  studies  where 
very  large  amplitudes  of  oscillation  will  be  of  interest.  In  such  cases,  where  it 
may  be  difficult  to  obtain  an  analog  output  of  the  displacement  tine  history,  a  light 
source  on  the  oscillating  model  may  be  used  to  indicate  model  position.  Orlik-Ruckemann 
presents  a  data  analysis  technique5 *  to  obtain  instantaneous  reduction  of  free-oscillation 
data  from  this  type  of  experimental  arrangement. 

In  the  simplest  form,  the  experimental  setup  consists  of  an  oscillating  light  source 
seen  through  a  stationary  window  (Pig. 57)  which  is  positioned  such  that  the  light  source 
can  be  seen  only  during  a  portion  of  a  cycle.  The  total  visible  travel  of  the  light 
source  in  relation  to  the  window  is  recorded  by  passing  the  light  through  a  large 
number  of  transverse  narrow  slots  to  a  photomultiplier  tube  to  record  the  number  of 
pulses. 

With  the  travel  summation  method,  the  logarithmic  decrement  can  be  expressed  as 
A  =  ( 2/s) ( 1  -  b  ♦  b  loge  b) 

^xo  &  ~  ~  (xo  ~ 


where 


8 


(87) 


60 


and  b  is  a  constant  with  convenient  values  around  0. 5,  depending  on  the  experimental 
arrangement.  Equation  (86)  thus  shows  that  the  logarithmic  decrement  is  simply  obtained 
by  dividing  a  constant  by  the  dimensionless  sum  of  the  travel  segments.  By  proper 
selection  of  b  and  a  slotted  window,  the  total  error,  including  the  approximations 
involved  in  retaining  terms  in  the  derivation  of  Equation  (87),  can  be  kept  under  1% 
for  most  cases. 

An  alternative  method  for  obtaining  the  logarithmic  decay  Involves  the  summation  of 
the  time  intervals  for  the  consecutive  cycles  of  oscillation  during  which  time  the  light 
source  is  between  x0  and  xa  .  This  method  is  less  accurate  than  the  travel  summation 
method,  because  of  approximations  used  in  the  derivation,  but  may  be  preferable  since 
it  permits  simpler  experimental  apparatus.  These  data  reduction  methods  are  described 
in  detail  in  Reference  54. 

3.2.2. 1.3  Logarithmic  Circuit 

A  recording  oscillograph  is  frequently  used  to  record  the  decaying  signal  of  model 
motion  from  a  displacement  pickup.  These  data  which  appear  as  a  harmonic  oscillation 
with  damping  must  be  replotted  in  terms  of  the  logarithm  of  the  amplitude  versus  a 
time  function  or  cycles  in  order  to  extract  the  damping  constant  from  the  results.  A 
method  which  effects  a  considerable  saving  in  time  in  the  analysis  of  the  data  traces 
consists  of  employing  a  logarithmic  distorting  circuit6,  which  converts  the  exponential 
envelope  of  oscillation  into  a  linear  envelope  which  is  proportional  to  the  number  of 
cycles  of  oscillation.  The  damping  constant  is  then  determined  from  the  slope  of  the 
envelope. 

3. 2. 2.1.4  Integration  Method 

Signals  from  displacement  pickups  recorded  on  magnetic  tape  may  be  analyzed  quickly 
without  the  aid  of  a  computer  by  making  use  of  the  integration  method  described  in  more 
detail  in  Reference  6.  To  employ  the  method,  a  series  of  complete  cycles  is  selected 
from  the  magnetic  tape  recording  and  the  signal  is  rectified  and  fed  into  an  analog 
Integrator  which  presents  the  integral  as  a  reading  R:  on  a  d.  c.  voltmeter.  The 
process  is  repeated  for  a  second  series  of  complete  cycles  to  obtain  the  second  Integral 
represented  by  R2  and,  if  n  is  the  number  of  cycles  between  sample  R1  and  R2  . 
the  logarithmic  decrement  is  given  by 


Hie  Integration  method  has  the  advantage  of  giving  weight  to  every  cycle  of  oscillation 
in  the  range  covered,  as  opposed  to  measuring  the  number  of  cycles  to  damp  to  a  given 
amplitude  ratio.  If  noise  exists  on  the  displacement  signal  the  integration  method 
tends  to  reduce  the  effects  it  may  have  on  the  damping  measurements. 

3.2.2.15  Other  Methods 

A  device  called  a  damping  trace  reader  has  been  developed  by  Leach65  to  measure  the 
logarithmic  decrement  from  oscillograph  traces  with  about  the  same  accuracy  as  for 
hand  data  reduction 
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J.2.2.2  Forced-Oscillation  Technique 

Testing  facilities  have  different  requirements  for  dynamic  stability  balance  systems. 
This  naturally  has  resulted  in  the  development  of  a  variety  of  different  systems  now 
in  use  in  the  United  States  and  in  European  countries.  Space  will  not  permit  des¬ 
criptions  of  all  systems  and  therefore  only  a  few  representative  one  will  be  ^iscussod 
to  show  important  features  and  techniques. 

3. 2.2. 2.1  AEDC-VKF 

Forced-oscillation  dynamic  stability  balance  systems  have  been  developed  at  the 
AEDC-VKF  for  operation  in  three  continuous  flow  tunnels,  which  cover  the  Mach  number 
range  from  1.5  to  12. 

The  dynamic  balance  system  (Fig. 58)  consists  of  a  crossed-flexure  pivot  mounted  on 
a  sting  and  connected  through  a  torque  member  and  linkage  to  an  electromagnetic  shaker 
housed  in  the  enlarged  aft  section  of  the  sting.  Details  of  the  crossed  flexures, 
input  moment  linkage,  and  model  locking  device  are  shown  in  Figure  59.  The  sting  and 
shaker  housing  contain  passages  for  water  cooling  which  is  necessary  for  tests  at 
hypersonic  speed.  Cooling  of  the  crossed  flexures  is  accomplished  through  the  use  of 
cool  air  sprays  which  are  turned  off  while  data  are  being  recorded.  Water-cooled 
jackets  which  surround  the  crossed  flexure  are  also  used  in  seme  applications. 

The  VKF  forced-oscillation  balance  is  equipped  with  an  amplitude  control  system 
which  allows  testing  of  both  dynamically  stable  and  unstable  models.  Testing  i$  normally 
accomplished  at  the  resonant  frequency  of  the  model  balance  system;  however,  the  control 
system  does  provide  a  means  of  resolving  the  in-phase  and  out-of-pLaae  (with  model 
displacement)  components  of  the  forcing  torque.  Because  of  the  low  structural  damping 
level  of  the  balance  and  the  aerodynamic  damping  levels  present  at  supersouic  and  hyper¬ 
sonic  speeds,  the  balance  system  is  always  operated  at  resonant  or  near- resonant 
conditions. 

It  was  shown  in  Section  3.  2.  1.  4  that  when  a  model  is  oscillated  at  a  constant 
amplitude  and  at  the  undamped  natural  resonant  frequency  of  the  system,  the  inertia 
moment  exactly  balances  the  restoring  moment  and  the  forcing  moment  is  equal  to  the 
damping  moment. 

A  block  diagram  of  the  control  and  data  acquisition  system  utilized  at  the  VKF 
system  is  given  in  Figure  60.  The  control  system  can  be  operated  in  either  an  open- 
look  mode  for  testing  dynamically  stable  configurations  or  in  a  closed- loop  mode  for 
testing  either  dynamically  stable  or  unstable  configurations  In  either  mode  of  opera¬ 
tion  the  frequency  of  oscillation  is  controlled  by  the  frequency  of  the  reference 
oscillator. 

Open-loop  operation  is  obtained  by  using  the  amplified  reference  oscillator  signal 
to  drive  the  electromagnetic  shaker  motor.  Resistance  ctrain-gsge  bridges  located  on 
the  input  torque  member  and  on  a  cross  flexure  (model  displacement  provide  analog 
signals  of  the  forcing  torque  to  the  model  and  the  model  displacement.  The  electrical 
output  of  the  bridges  is  amplified  by  a  carrier  Amplifier  to  provide  appropriate  ampli¬ 
tudes  of  the  torque  and  disp^ement  signals  necessary  for  the  data  acquisition  system. 
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The  closed-loop  mode  of  operation  operates  on  the  amplified  error  principle.  The 
reference  signal  and  the  displacement  signal,  which  is  phase-shifted  approximately 
90  degrees,  are  compared  in  a  difference  amplifier.  This  difference  or  error  signal 
is  amplified  and  used  to  drive  the  shaker  motor.  The  error  signal  will  remain  in 
phase  with  the  reference  signal  for  all  conditions  requiring  a  positive  application 
of  torque.  If  the  model  is  aerodynamically  unstable,  the  error  signal  will  be  180 
degrees  out  of  phase  with  the  reference  signal  and  a  negative  torque  will  be  supplied 
to  the  model. 

Due  to  high  frequency  vibration  in  the  wind  tunnels,  signal  conditioning  is  required 
to  remove  the  noise  from  the  torque  and  displacement  signals  before  recording  on  the 
data  acquisition  system.  Two  bandpass  filters  tuned  to  the  frequency  of  the  reference 
oscillator  are  provided.  The  filters  and  other  frequency-sensitive  circuits  track  with 
the  frequency  selected  by  the  reference  oscillator.  The  output  of  the  bandpass  filters 
is  amplified  to  provide  the  current  necessary  for  the  recording  instruments. 

The  data  acquisition  system  consists  of  an  oscillograph,  resolvers,  on-line  data 
system,  analog  magnetic  tape  system,  and  Beckman  210  analog-to-digital  converter  system. 
The  torque  and  displacement  signals  are  applied  to  the  inputs  of  all  the  recording 
instruments  and  to  an  oscilloscope  di  aying  a  Lissajous  pattern.  The  Lissajous 
pattern  provides  an  indication  of  the  90  degree  phase  angle  relationship  between  these 
two  signals  necessary  to  the  condition  of  resonance  of  the  model  balance  system. 

Under  certain  tunnel  conditions,  and  for  certain  model  configurations,  it  is  bene¬ 
ficial  to  operate  at  a  frequency  displaced  from  the  resonant  frequency.  The  resolver 
system  is  utilized  to  determine  the  in-phase  and  quadrature  components  of  the  forcing 
torque  with  respect  to  che  model  displacement. 

The  in-phase  component  is  obtained  by  multiplying  the  torque  and  displacement  signal 
with  an  analog  multiplier.  The  quadrature  component  is  obtained  by  phase-shifting  the 
torque  signal  90  degrees  and  then  multiplying  the  phase-shifted  signal  with  the  dis¬ 
placement  signal.  The  output  signals  of  the  multipliers  contain  a  d.c.  component  and 
a  double  frequency  component.  By  filtering  out  the  double  frequency  component,  a  d.c. 
value  is  obtained  which  is  proportional  to  the  in-phase  and  quadrature  torque  components. 

The  in-phase  torque  component  is  zero  at  resonance,  and  the  quadrature  component  is 
a  maximum  value.  A  d.c.  voltmeter  is  utilized  to  monitor  the  in-phase  component  and 
serves  as  an  indication  of  a  resonance  condition.  Both  torque  components  are  recorded 
on  analog  magnetic  tape,  the  Beckman  210  converter,  and  the  on-line  data  system. 

Data  recording  systems  utilized  for  forced-oscillation  testing  are  the  oscillograph, 
analog  magnetic  tape,  Beckman  210  converter,  and  on-line  data  system.  The  oscillograph 
and  analog  magnetic  tape  provide  permanent  recordings  of  the  data  as  obtained  in  the 
wind  tunnel.  The  Beckman  210  analog-to-digital  converter  provides  digital  conversion 
of  the  analog  data  and  records  the  digital  values  on  magnetic  tape,  which  is  then 
processed  by  a  high  speed  computer  to  obtain  final  data. 

Reduced  data  are  desirable,  and  in  some  cases  necessary,  while  the  test  is  in 
progress.  The  on-line  data  handling  system  provides  immediate  examination  of  the 
damping  derivatives  Here  the  term  on-line  signifies  that  the  test  data  are  fed 
directly  from  test  instruments  into  a  computer.  Therefore,  the  on-line  data  system 
includes  obtaining  analog  representation  of  parameters  to  be  measured,  conversion  of 
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these  parameters  to  equivalent  digital  values,  accumulation  in  a  digital  computer, 
computations  of  the  data,  and  the  presentation  of  the  calculated  data  in  derivative 
form  in  a  data  tabulation.  A  block  diagram  of  the  on-line  data  system  is  shown  in 
Figure  61.  A  console  containing  all  of  the  control  systems,  resolvers,  and  part  of 
the  on-line  data  system  is  shown  in  Figure  62. 

Included  iu  the  console  is  an  oscilloscope  for  monitoring  the  torque  and  displacement 
signals  in  either  normal  or  Lissajous  patterns  to  determine  the  resonant  condition. 

An  oscilloscope  is  also  used  to  set  up  both  the  automatic  control  system  and  the  re¬ 
solver  network,  and  to  monitor  the  reference  signal,  error  signal,  or  output  voltage 
to  the  shaker  motor.  A  root -mean- square  and  peak-reading  voltmeter  is  used  to  monitor 
the  output  voltage  and  thus  the  amount  of  power  or  watts  being  delivered  to  the  shaker 
motor.  It  can  also  be  used  to  obtain  either  root -mean- square  or  peak  voltage  of  the 
torque  and  displacement  signals.  The  ratio  of  these  signal  peaks  is  used  to  determine 
the  damping  derivatives. 

The  complete  resolver  system  mounted  in  the  console  includes  all  the  amplifiers, 
filters,  etc. ,  required  for  resolving  the  in-phase  and  out-of-phase  torque  components. 

One  centerscale  d.c.  voltmeter  is  used  to  read  the  in-phase  component  of  the  torque 
resolver  network.  Additional  information  on  the  system  is  given  in  Reference  56. 

J.2.2.2.2  AEDC-PWT 

The  forced-oscillation  balance  system  in  use  in  the  AEDC-PWT  is  shown  in  Figure  63. 

The  balance  has  a  bearing  pivot  with  a  hydraulic-controlled  mechanism  capable  of  pro¬ 
ducing  a  controlled  pitching  motion  to  the  model.  Angular  displacement  is  sensed  by 
strain  gages  on  the  flexure  pivot,  and  the  moment  is  measured  between  the  pivot  and 
the  actuating  mechanism. 

The  constant  amplitude,  sinusoidal  oscillation  is  controlled  by  feedback  systems 
in  which  the  displacement  signals  are  compared  with  a  command  signal,  and  frequency 
is  controlled  by  the  hydraulic  pump  motor  speed.  Figure  64  is  a  block  diagram  showing 
the  control  and  readout  instrumentation.  Oscillation  amplitude  and  frequency  are 
selected  by  settings  a  two-phase  oscillator  which  is  used  as  a  phase  reference  for 
model  motion  and  input  moment.  Deviations  from  a  set  amplitude  generate  an  error 
signal  which  is  a*nplified  and  results  in  a  system  drive  signal. 

In  addition  to  in-phase  signals,  the  two-phase  oscillator  also  generates  an  out-of¬ 
phase  signal  which  leads  the  in-phase  signal  by  90  degrees.  The  in-phase  and  out-of- 
phase  oscillator  signals  provide  readout  gage  excitations  for  the  moment  gage  and 
position  gage  to  form  a  sinusoidal  oscillation  of  the  four  readout  gage  outputs  about 
a  d.c.  level. 

Each  gage  output  wave  form  is  integrated  over  a  number  of  cycles,  the  a. c.  component 
is  averaged  out,  and  the  d.c.  component  remains  as  a  signal  datum.  Model  position  and 
applied  torque  are  proportional  to  the  signal  data  and  when  used  with  a  phase  diagram 
provide  relationships  for  determining  the  dynamic  stability  and  static  stability 
derivatives.  The  system  readouts  are  fed  to  a  computer  for  data  reduction  and  display 
of  the  aerodynamic  derivatives  in  tabular  or  graphical  form  while  the  test  is  in  progress. 

The  basic  theory,  description  of  the  physical  components,  theory  of  the  control  and 
readout  system,  and  derivation  of  the  associated  stability  parameters  are  presented  in 
References  57  and  58. 
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3. 2.  2. 2.3  NASA  -  LRC 

The  forced-oscillation  balance  system  used  by  NASA-LRC  is  shown  in  a  exploded  and 
assembled  view  in  Figure  65.  Sketches  of  the  balance  system  are  presented  in  Figures 
66  and  67.  Sinusoidal,  constant-amplitude  motion  about  an  oscillation  axis  is  produced 
by  a  mechanical  Scotch  yoke  and  crank  which  is  driven  by  a  variable- frequency  electric 
motor  in  the  downstrean  portion  of  the  sting.  The  stiff  strain-gaged  beam  located 
between  the  model  attachment  bulkhead  and  pivot  axis  is  used  to  measure  the  forcing 
moment  required  to  oscillate  the  model.  Because  of  its  location  the  signal  output  is 
not  influenced  by  the  friction  or  mechanical  play  in  the  system.  Model  displacement 
is  measured  with  a  strain-gaged  cantilever  beam  mounted  between  the  fixed  sting  and 
the  oscillating  model  support.  Springs  of  different  stiffnesses  are  used  to  provide 
a  range  of  resonant  frequencies. 

Outputs  from  the  displacement  and  forcing-moment  strain  gages  pass  through  coupled 
sine-cosine  resolvers  which  are  geared  to  the  drive  shaft  and  thus  rotate  at  the  funda¬ 
mental  drive  frequency  (Fig. 66).  Oscillation  frequency  is  measured  with  a  signal  genera¬ 
tor  coupled  to  the  drive  shaft.  In  operation  the  resolvers  transform  the  moment  and 
amplitude  functions  into  orthogonal  components  which  are  measured  with  d.c.  micro¬ 
ammeters.  From  these  components  the  resultant  -.v'plied  moment,  displacement,  and  phase 
angle  are  found  and.  with  the  oscillation  frequency,  the  aerodynamic  static  and  damping 
moments  are  computed.  A  block  diagram  of  the  electronic  circuit  used  to  measure  the 
applied  moment  is  shown  in  Figure  68  and  a  similar  circuit  is  used  to  measure  the 
angular  displacement. 

Use  is  made  of  automatic  digital  equipment  to  acquire  data  quickly  and  eliminate 
operator  errors  due  to  incorrect  manual  recording  or  reading.  All  of  the  system  vari¬ 
ables  are  automatically  digitized,  and  a  high  speed  computer  is  used  to  perform  the 
data  reduction  and  present  results  in  coefficient  form. 

A  complete  description  of  the  system  electrical  relations  and  methods  involved  in 
the  operation  and  data  reduction  are  found  in  References  59  and  60. 

3.2.2.  2.4  RAE -Bedford 

The  forced-oscillation  testing  technique  has  been  employed  by  the  RAE,  Bedford,  in 
the  development  of  a  system  for  measuring  stability  derivatives  and  damping  derivatives 
due  to  angular  velocity  in  roll  and  yaw  on  sting-mounted  models62.  Derivatives  can  be 
obtained  for  motion  in  vertical  or  horizontal  translation  (heave  or  sideslip);  however, 
it  is  noted  in  Reference  62  that  the  apparatus  was  not  designed  for  this  purpose  and 
the  accuracy  of  the  measurements  is  less  than  that  required  for  most  tests.  Emphasis 
in  the  design  has  been  toward  a  small  size  system  that  can  be  used  for  testing  slender 
wings  of  fiber  glass  construction. 

In  some  respects  the  balance  system  is  similar  to  those  described  in  Section  3. 2. 2. 1 
in  that  the  model  is  mounted  on  a  sting  by  means  of  a  flexure  pivot  allowing  one  degree 
of  freedom  in  rotation.  Continuous  excitation  at  or  uear  the  resonant  frequency  of 
the  system  is  provided  by  a  moving-coil  permanent -magnet  vibration  generator  mounted 
on  the  sting  and  coupled  to  a  driving  rod  passing  through  the  sting. 

A  schematic  arrangement  of  various  parts  of  the  balance  for  yaw-sideslip  oscillation 
is  shown  in  Figure  69.  The  balance  is  mounted  on  the  end  of  a  sting  by  means  of  a 
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flexure  pivot  to  allow  rotation.  A  view  of  the  flexure  pivot  and  cylindrical  section 
of  the  balance  which  mounts  inside  of  the  model  is  shown  in  Figure  69,  and  photographs 
are  presented  in  Figure  70.  Included  in  the  flexure  pivot  assembly  is  the  drive  linkage 
to  convert  the  longitudinal  force  in  the  driving  rod  into  a  couple  acting  between  the 
end  of  the  sting  and  the  model. 

The  balance  for  measuring  oscillation  in  roll  is  illustrated  schematically  in  Figure 
69.  The  unit  is  of  one-piece  construction  and  is  instrumented  with  yaw  and  roll 
accelerometers.  The  system  has  one  degree  of  freedom,  and  excitation  in  roll  is  pro¬ 
vided  by  applying  torque  through  the  driving  rod  by  means  of  linkage  from  the  vibration 
generator. 

Although  the  flexure  pivot  is  located  near  the  center  of  pressure  of  the  wing  models, 
in  order  to  avoid  large  deflections  under  steady  aerodynamic  loads,  the  flexibility  of 
the  sting  will  allow  some  translation  of  the  pivot  axis  which  may  be  appreciable  for 
some  cases.  To  allow  for  displacement  of  the  pivot  axis,  for  example,  when  measuring 
pitching  derivatives,  heaving  derivatives  are  measured  by  exciting  a  vertical  oscilla¬ 
tion  (which  will  be  accompanied  by  some  angular  movement)  and  measuring  the  changes 
in  natural  frequency,  excitation  force,  and  angular  movement. 

The  motion  of  the  model  on  the  balance  system  consists  of  a  pitching  component, 
generated  in  part  by  the  rotation  of  the  model  about  the  flexure  pivot  and  in  part  by 
the  bending  of  the  end  of  the  sting.  These  two  components  are  present  in  different 
proportions  in  the  two  natural  modes  of  the  system.  The  pitching  mode  includes  a 
certain  amount  of  heave  motion  of  the  pivot  axis,  and  the  heave  mode  includes  an 
appreciable  amount  of  pitching  motion  and  similarly  for  yaw  and  sideslip. 

As  was  noted  in  Reference  62,  it  may  not  be  obvious  how  an  exciting  couple  acting 
between  the  model  and  the  end  of  the  sting  can  produce  no  deflection.  The  reason  is 
that  in  the  natural  mode  of  oscillation  in  heave,  for  example,  there  is  a  certain 
amount  of  relative  angular  movement  between  the  model  and  the  end  of  the  sting  (even 
when  there  is  no  pitching  component  in  the  oscillation).  It  is  thus  possible  for  the 
excitation  couple  to  feed  energy  into  the  system  in  this  model. 

The  analysis  of  the  system  is  simplified  by  treating  the  system  as  having  two  degrees 
of  freedom  with  two  natural  modes  of  vibration  -  one  predominantly  pitching  and  the 
other  predominantly  heaving.  The  heaving  derivatives  appear  as  a  by-product  of  the 
analysis  to  obtain  the  damping- in-pitch  derivatives.  Yawing  derivatives  are  measured 
in  the  same  manner,  with  sideslip  derivatives  appearing  as  a  by-product. 

The  basic  measurements  made  with  the  balance  system  include  (i)  the  excitation 
frequency  which  is  measured  by  an  electronic  counter,  (ii)  the  excitation  force  which 
is  measured  in  terms  of  the  exciter  current,  and  (iii)  the  balance  displacements  which 
could  be  obtained  as  measurements  of  strains  or  relative  displacements;  accelerometers 
were  selected  foi  the  design  since  they  are  not  affected  by  steady  loads  and  are  com¬ 
paratively  free  from  interactions. 

The  balanc *  is  designed  (Fig. 70)  with  strain  gages  mounted  on  sections  for  measuring 
static  forces  aLd  moments,  except  for  rolling  moment  which  is  more  conveniently  measured 
with  the  flexure  unit  designed  for  roll  oscillations. 


66 


The  automatic  measuring  and  recording  equipment  developed  by  the  RAE  for  experiments 
concerned  with  oscillatory  derivative  is  a  digital  amplitude  and  phase  measuring  system 
(DAPHNE).  It  was  designed  to  use  the  wattmeter  method6  to  resolve  the  input  signals 
into  in-phase  and  quadrature  components,  and  accurate  measurements  of  frequency  are 
made  with  an  electronic  counter. 

3. 2. 2. 2.5  NASA -ARC 

A  method  and  a  balance  system  is  described  in  Reference  63  for  measuring  the  lateral 
and  dynamic  stability  derivatives.  The  features  of  the  system  include  single-degree- 
of-freedom  oscillations  which  are  used  to  obtain  components  of  pitch,  yaw,  and  roll  by 
varying  the  axis  of  oscillation.  Amplitude  control  is  achieved  with  a  velocity  feedback 
loop.  The  system  is  equipped  with  a  system  of  strain  gage  data  processing  in  which 
electronic  analog  computer  elements  are  used  in  measuring  the  amplitude  and  phase 
position.  Moment  derivatives  which  arise  from  angular  motions  can  be  measured  with 
the  balance  system.  These  include  the  rotary  damping  derivatives  Cnp  ,  CBQ  +  C>6t  . 
and  Cnr  -  Cn/§  ;  the  cross  derivatives  Cnp  and  Cjr  -  C ^  ;  and  the  displacement 
derivatives  CB(X  ,  and  C(j9  . 

The  forced-oscillation  system  has  two  oscillation  mechanisms  -  one  for  pure  pitching 
or  yawing  motion  in  which  the  oscillation  axis  is  perpendicular  to  the  longitudinal 
axis  of  the  sting,  and  one  for  combined  rolling  and  pitching  or  rolling  and  yawing  in 
which  the  axis  of  oscillation  is  inclined  45  degrees  to  the  longitudinal  axis  uf  the 
sting.  The  electromagnetic  shaker  is  housed  in  an  enlarged  section  of  the  sting  down¬ 
stream  of  the  model  and  connected  to  the  flexure  pivots  by  push  rods.  Strain  gages 
were  used  to  measure  the  oscillation  amplitude  and  supply  inputs  to  the  analog  computer. 

The  theory  of  the  system  operation,  feedback  control,  analog  computing  system,  and 
description  of  the  apparatus  is  presented  in  Reference  t?3. 

3. 2. 2. 2. 6  University  of  Michigan 

Results  obtained  at  hypersonic  and  hypervelocity  speeds  have  shown  the  importance 
of  the  reduced  frequency  of  oscillation.  With  the  use  of  the  free-oscillation  testing 
technique,  the  frequency  of  oscillation  can  be  varied  only  by  changing  either  the 
moment  of  inertia  or  the  flexure  spring  constant.  Thus  coverage  of  a  wide  range  of 
frequencies  requires  several  time-consuming  changes  in  the  test  equipment. 

The  University  of  Michigan  has  developed  an  electromagnetic,  forced-oscillation 
method  which  enables  a  range  of  oscillation  frequencies  to  be  investigated  with  ease64. 
The  method  utilizes  the  principle  of  the  moving  coil  galvanometer.  An  alternating 
current,  passing  through  a  coil  embedded  within  the  test  model,  interacts  with  a  static 
magnetic  field  Imposed  by  a  coil  system  mounted  in  one  plane  of  the  test  section.  This 
causes  the  model  to  oscillate  in  a  sinusoidal  manner  following  the  existing  torque  but 
out  of  phase  with  it.  Prom  the  phase  lag  of  the  motion  relative  to  the  existing  moment 
and  wind-on  and  wind-off  results,  it  is  possible  to  determine  the  aerodynamic  damping. 

A  schematic  view  of  the  test  setup  showing  the  external  coil  is  presented  in  Figure 
71.  A  circular  exciter  coil  was  mounted  coaxially  within  the  Plexiglas  model  which  was 
mounted  on  a  cruciform  flexure  pivot  shown  in  Figure  72. 
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The  methods  have  been  used  for  measuring  the  position  of  the  model.  One  of  the 
methods  utilized  the  output  of  the  strain  gage  bridge  shown  in  Figure  72  and  the  other 
was  based  on  the  response  of  a  second  circular  coil  within  the  exciter  coil  located  in 
the  model.  As  the  model  oscillates  in  the  magnetic  field,  an  electromotive  force  is 
induced  in  the  secondary  coil  and  may  be  used  for  sensing  the  cone  position.  Although 
this  technique  hes  been  successful,  it  has  the  drawback  that  the  data  reduction  pro¬ 
cedure  is  very  lengthy,  which  may  preclude  on-line  reduction.  Use  of  strain  gages  to 
measure  the  model  position  eliminates  the  difficulties  associated  with  the  two-coil 
techniques. 

The  technique  of  electromagnetic  excitation  and  using  the  phase  angle  to  determine 
the  damping  has  limitations  which  may  prohibit  its  use  for  some  configurations.  Although 
the  structural  damping  was  small,  the  phase  angle  measurements  at  Mach  number  2.5  on  a 
25  degree  half-angle  cone  were  less  than  one  degree,  which  makes  accurate  measurements 
of  aerodynamic  damping  at  this  and  higher  Mach  numbers  very  difficult  but  not  impossible. 
Improvements  in  the  accuracy  of  the  damping  measurements  will  result  if  the  model  size 
is  increased,  the  pivot  axis  is  moved  as  far  forward  as  possible,  and  a  gas  bearing  is 
used  to  eliminate  amplitude  restrictions  and  structural  damping  due  to  the  flexure. 

In  addition  the  gas  bearing  would  eliminate  sting  vibrations  that  influence  damping 
measurements.  Details  on  the  instrumentation,  data  reduction  procedures,  and  sample 
test  results  obtained  with  electromagnetic  excitation  are  presented  in  Reference  64. 

3. 2.  2. 3  Magnetic  Model  Suspension  System. 

Systems  for  captive  model  aerodynamic  testing  require  support  of  the  model  by  use 
of  a  sting  into  the  model  base  region,  a  side-mounted  strut,  or  even  thin  wires.  Each 
of  these  techniques  may  introduce  interference  effects  which  can  in  some  cases  have  a 
significant  influence  on  the  results.  Support  interference  effects  must  be  studied 
experimentally  to  gain  some  idea  of  the  magnitude  of  the  influence,  since  the  problem 
has  not  been  reduced  to  theoretical  analysis. 

Magnetic  suspension  of  wind  tunnel  models  provides  one  means  for  eliminating  the 
influence  of  the  support  system.  The  first  successful  suspension  of  a  wind  tunnel 
model  that  enabled  the  measurement  of  interference-free  drag  was  accomplished  at  ONERA, 
France,  in  1957  (Ref. 65).  Since  that  time,  several  laboratories  in  France,  England, 
and  the  United  States  have  been  engaged  In  the  development  of  magnetic  suspension 
systems  for  different  purposes.  The  potential  for  magnetic  suspension  is  very  great 
when  it  is  considered  that  not  only  is  it  an  interference- free  support  system,  but  it 
is  a  balance  system  for  static  forces  and  moments  and  dynamic  measurements,  since  model 
motion  can  be  controlled. 

For  example,  control  signals  may  be  introduced  into  a  feedback  control  circuit  to 
produce  model  motion  without  alteration  of  the  basic  system.  In  principle,  models  can 
be  subjected  to  rigidly  forced  motion  either  in  sinusoidal  or  random  modes  or  rotational 
motion  with  the  translational  degrees  of  freedom  fixed.  With  the  latter  mode  of  opera¬ 
tion.  dynamic  damping  data  in  pitch  or  yaw  are  obtained  by  either  the  free-  or  forced- 
oscillation  techniques.  Roll  damping  measurements  may  be  made  on  axisyumetric  and 
finned  bodies  by  using  the  spin  decay  technique. 

The  development  of  magnetic  suspension  systems  has  progressed  very  rapidly  in  the 
past  several  years;  however,  it  has  not  reached  the  stage  that  the  suspension  systems 
compete  favorably  with  the  more  conventional  techniques  for  obtaining  static  and  dynamic 
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stability  data.  Magnetic  suspension  systems  are,  however,  used  routinely  for  wake 
studies  and  some  drag  measurements.  Although  research  to  develop  techniques  for 
making  dynamic  stability  measurements  with  a  magnetic  model  suspension  system  is  in 
the  formative  stages,  the  status  of  the  work  accomplished  to  date  is  briefly  summarized 
to  point  out  the  major  areas  of  development.  Treatment  of  the  equations  of  motion 
and  data  reduction  have  been  omitted,  since  each  depends  more  or  less  on  the  specific 
system;  these  topics  are  covered  in  the  references. 

The  suspension  system  in  use  at  the  Gas  Dynamics  laboratory  at  Princeton  University66 
was  developed  to  support  models  for  experimental  research  on  hypersonic  wakes  of  axisym- 
metric  bodies.  Since  the  configurations  of  interest  ranged  from  spheres  to  slender 
bodies,  a  three-degree- of- freedom  system  and  a  one-point  suspension  was  used  (Fig. 73). 
The  system  is  mounted  vertically  to  use  the  gravity  force  to  oppose  drag,  to  help  the 
magnetic  suspension.  Models  are  made  with  ferrite  spheres  as  the  magnetic  element  to 

be  acted  on  by  the  magnetic  field.  A  conical  model  with  a  ferrite  sphere  embedded  in 

the  body  is  shown  in  Figure  74.  When  only  spheres  are  tested,  the  ferrite  is  used 

without  any  other  material.  Although  the  Princeton  system  is  not  used  for  dynamic 

stability  measurements,  such  data  can  be  obtained  with  the  system  and  model  shown  in 
Figure  74  by  making  use  of  the  free-oscillation  technique,  provided  position  sensing 
is  incorporated  in  the  system  for  model  motion  histories.  Completely  free  rotational 
model  motion  is  also  possible  with  the  balance  system  developed  by  the  University  of 
Virginia67,  provided  position  sensing  is  incorporated. 

The  MIT  magnetic  model  suspension  system  described  in  Reference  31  and  shown  in 
Figure  75  has  been  used  to  obtain  dynamic  stability  data  from  wind  tunnel  measurements 
at  supersonic  speeds.  The  forced-oscillation  technique  was  employed  to  measure  damping- 
in-pitch  derivatives  on  a  model  undergoing  a  forced  motion  by  measuring  the  change  in 
the  phase  of  the  motion.  Use  of  this  technique  requires  that  the  slope  of  the  aero¬ 
dynamic  moment  curve,  M^  ,  be  known  and  sinusoidal  motion  imposed  by  the  forcing 
function.  Although  no  attempt  was  made31  to  measure  damping  due  to  plunging,  a  method 
is  described  in  the  reference  whereby  such  measurements  could  be  obtained. 

The  RAE  magnetic  model  suspension  system  in  use  at  Mach  9  is  shown  in  schematic 
form  in  Figure  76,  and  a  view  of  the  system  installed  in  the  tunnel  is  shown  in  Figure 
77.  The  system,  now  in  the  final  stages  of  development,  was  designed  to  measure  static 
forces  and  moments  and  dynamic  stability.  Additional  details  on  the  system  may  be 
found  in  References  68  and  69. 

The  magnetic  suspension  system  at  Southampton  University70,  71  makes  use  of  three 
horseshoe-shaped  electromagnets  and  one  hollow  solenoid- like  coil  for  control  of  the 
model  in  six  degrees  of  freedom;  the  model  itself  contains  a  permanent  magnet.  A 
schematic  arrangement  of  these  magnets  is  shown  in  Figure  78,  and  a  model  subjected  to 
a  forced  oscillation  in  roll  while  suspended  in  the  six-component  balance  is  shown  in 
Figure  79.  The  system  has  been  designed  for  the  measurement  of  aerodynamic  derivatives, 
and  emphasis  has  been  placed  on  tests  of  winged  models  which  requires  proper  control 
over  the  model  roll  attitude.  Although  there  are  several  different  methods  that  have 
been  investigated  for  controlling  model  roll  attitude,  shaped  cores  within  the  model 
have  been  investigated  most  extensively.  The  usual  practice  in  magnetic  suspension 
systems  is  to  have  a  magnetic  core  arranged  to  lie  along  the  model  axis.  These  cores 
are  generally  circular  in  cross  section  and  thus  do  not  impart  any  resistance  to  dis¬ 
turbances  in  roll.  However,  if  the  core  is  constructed  with  a  noncircular  cross  section, 
the  core  will  adopt  one  or  more  preferred  attitudes  when  suspended  in  the  field  of  the 
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model  supporting  magnet  B  shown  in  Figure  78.  As  pointed  out  in  Reference  71,  a  model 
containing  a  shaped  core  can  be  induced  to  run  as  a  synchronous  motor  by  varying  currents 
to  the  lateral  magnets;  however,  the  model  would  have  to  be  spun  up  to  the  synchronous 
speed  to  allow  lock  on.  Speeds  greater  than  1000  r.p.m.  have  been  produced  on  a  0.75 
in.  diameter  body  of  revolution  at  the  University  of  Southampton.  This  technique  could 
be  used  for  measuring  roll  derivatives. 

Reference  should  be  made  to  Reference  29  for  information  on  magnets,  optics  systems, 
control  circuits,  operating  techniques,  equations  of  motion,  data  reduction,  and  problem 
areas. 

3. 2. 3  Representative  Systems  for  Tests  of  Control  and 
Lifting  Surfaces 

Methods  for  experimentally  obtaining  dynamic  stability  data  on  models  of  missiles 
and  aircraft  may  also  be  applied  to  investigations  of  the  aerodynamic  characteristics 
of  control  surfaces,  such  as  rudders,  elevators,  and  ailerons.  Although  aerodynamic 
theory  has  been  developed  to  analyze  control  surface  behavior,  recourse  to  experi¬ 
mental  verification  is  necessary  because  of  nonlinear  effects  which  may  be  a  function 
not  only  of  amplitude  and  frequency  but  also  of  Reynolds  number  as  associated  with 
flow  separation  ahead  of  the  control  surface. 

Experimental  investigations  to  date  have  been  designed  primarily  to  measure  the 
aerodynamic  characteristics  (hinge  moment  and  damping)  of  control  surfaces  oscillating 
at  transonic  and  low  supersonic  speeds  where  the  theories  show  a  strong  influence  of 
oscillation  frequency  and  Mach  number72. 

Since  the  specific  details  on  the  design  of  the  model  control  surface  and  dynamic 
balance  system  will  depend  on  the  ingenuity  of  the  designer,  no  attempt  will  be  made 
to  treat  detailed  problems  peculiar  to  a  particular  system.  Reference  will,  however, 
be  made  to  typical  systems  to  describe  the  general  techniques  that  may  be  used  to 
obtain  control  surface  damping  and  hinge-moment  data.  Representative  systems  for 
obtaining  damping  data  on  lifting  surfaces  are  also  presented. 

3.2.3. 1  Free-  and  Forced-Oscillation  Techniques 

The  methods  employed  and  the  data  reduction,  in  general,  are  the  same  as  those  des¬ 
cribed  in  Section  3.2.1  for  the  free-  and  forced-oscillation  techniques  applied  to 
control  surfaces  and  half  models.  Applications  of  these  techniques  to  a  control  surface 
mounted  on  the  trailing  edge  of  a  triangular  wing  (Fig. 80)  are  presented  in  Reference 
73.  The  control  surface  was  attached  to  the  wing  by  means  of  a  flat  spring  running 
the  full  length  of  the  control  surface.  Excitation  was  supplied  with  a  torque  motor 
mounted  in  the  model,  and  control  surface  position  was  measured  with  strain  gages 
mounted  on  cross  flexures.  These  were  used  to  increase  the  natural  frequency  of  the 
system.  To  operate  in  the  forced-oscillation  mode,  the  control  surface  was  driven  at 
the  resonant  frequency  of  the  system  where  the  torque  required  to  sustain  a  given 
amplitude  of  oscillation  is  directly  proportional  to  the  damping  of  the  control  surface. 
Thus  the  damping  jould  be  calculated  from  a  measurement  of  the  armature  current  of  the 
torque  motor. 
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For  the  free-oscillation  tests,  the  system  was  driven  at  the  resonant  frequency  of 
the  system  to  a  given  amplitude.  The  armature  circuit  of  the  torque  motor  was  then 
opened  and  the  ensuing  decay  in  amplitude  was  measured  to  obtain  the  damping  from  the 
log  decrement. 

Although  oscillating  at  resonance  to  obtain  the  initial  displacement  has  the  advantage 
that  the  torque  to  deflect  the  control  surface  is  small,  it  has  the  disadvantage  of 
increasing  the  mass  moment  of  inertia  of  the  system.  Initial  amplitudes  for  free- 
oscillation  motion  may  be  obtained  with  a  static  or  impulsive  deflection  of  the  surface 
to  the  desired  amplitude,  but  this  technique  has  the  disadvantage  that  the  localized 
loading  may  lead  to  some  distortion  of  the  control  surface. 

Forced-oscillation  systems  used  for  damping  measurements,  or  as  a  means  to  achieve 
initial  amplitudes,  generally  require  considerably  more  space  compared  to  the  require¬ 
ments  for  the  free-oscillation  system  with  a  static  displacement  mechanism  and  therefore 
are  restricted  to  large  models.  Orlik-Ruckemann74  has  applied  the  forced  oscillation 
technique  to  obtain  initial  amplitudes  for  free-oscillation  tests  on  small  aerodynamic 
control  surfaces  by  developing  an  experimental  apparatus  for  use  with  half  models. 

The  apparatus  for  elevator  oscillation  tests  on  a  half- airplane  model  is  shown  in 
Figures  81  and  82.  In  operation,  the  control  surface  is  oscillated  by  a  pair  of  driving 
colls  to  a  desired  amplitude  and  at  the  resonance  frequency  of  the  oscillating  system 
and  then  the  circuit  is  interrupted  to  allow  a  free  oscillation.  For  this  design,  the 
control  surface  motion  is  measured  by  strain  gages  located  on  the  control  surface 
hinges  which  are.  in  this  case,  cantilever  beams.  The  signal  from  the  strain  gages, 
which  is  proportional  to  the  angular  displacement  of  the  co  rol  surface,  was  processed 
with  a  dampometer  to  obtain  the  characteristics  of  the  oscillation,  and  thus  the  aero¬ 
dynamic  data. 

Ilartz75  has  used  the  free-oscillation  technique  on  rocket-powered  free-flight  models 
to  obtain  hinge-moment  characteristics  for  a  full-span,  trailing-edge  control  surface 
on  a  delta  wing.  This  same  technique  can  also  be  used  in  wind  tunnel  tests.  The  con¬ 
trol  surface  was  mounted  to  the  wing  through  a  cantilever- type  flexure  hinge  which 
was  instrumented  with  strain  gages  to  measure  the  control  surface  deflection.  Initial 
displacement  of  the  control  surface  was  achieved  by  plucking  periodically  at  the  root 
chord  by  means  of  a  motor-driven  cam  shown  in  Figure  83. 

Orlik-Ruckemann  and  Laberge25  have  used  electromagnetic  excitation  for  free-oscilla¬ 
tion  dynamic  tests  on  a  series  of  delta  wings.  An  example  of  the  system  is  shown  in 
Figure  84.  This  system76  has  been  developed  to  the  staga  that  its  operation  is  com¬ 
pletely  automatic,  that  is.  the  output  signal  from  the  model  suspension,  which  is  a 
strain-gaged  cruciform  elastic  member,  and  the  input  signal  to  the  electromagnetic 
oscillator  form  a  closed-loop  circuit  in  which  the  frequency  is  automatically  adjusted 
to  the  natural  frequency  of  the  system.  After  the  circuit  is  interrupted,  the  model 
oscillates  freely;  data  are  recorded  on  the  dampometer.  and  the  cycle  is  again  repeated. 
This  automatic  feature  allows  acquisition  of  a  large  amount  of  data  within  a  short 
period  of  time  and  is  particularly  useful  for  blowdown  tunnels  which  operate  for  very 
short  time  periods  where  manual  adjustment  of  instrumentation  is  not  possible  without 
reducing  the  amount  of  data  that  may  be  acquired.  A  schematic  block  diagram  of  the 
system  control  and  data  acquisition  in  shown  in  Figure  85. 

A  half-model,  fixed-amplitude,  derivative  measurement  system  used  in  the  Aerodynamic 
Research  Laboratory  of  Hawker  Siddeley  Dynamics  Ltd.  is  shown  in  Figure  86.  The  system 
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is  described  in  detail  in  Reference  23.  The  apparatus  provides  measurements  of  the 
wing  and  control  surface  derivatives  due  to  models  performing  sinusoidal  motion  in 
pitch  and  vertical  translation.  A  nonresonant  method  is  used  to  reduce  the  electronic 
and  mechanical  syctem  complexity. 

The  balance  system  is  built  on  a  nominally  rigid  platform  which  is  connected  to 
earth  by  three  supports  and  force  detecting  transducers.  Angular,  simple  harmonic 
motion  is  achieved  with  a  variant  of  the  well-known  Scotch  Yoke  mechanism  which  converts 
constant  angular  motion  into  linear  harmonic  motion.  This  motion  ic  then  imparted 
through  linkage  to  the  model  shaft.  The  model  is  driven  by  an  electric  motor  through 
a  variable  speed  gear  box  so  that  the  frequency  may  be  varied  between  5  and  25  cycles 
per  second  at  amplitudes  of  2  degrees. 

The  oscillating  sy*.em  is  inertia  balanced  so  that  the  difference  between  the  wind-on 
and  wind-off  transducer  outputs  in  terms  of  amplitude  and  phase  relationship  with  respect 
to  model  motion  are  measured  for  oscillations  about  two  model  axes.  In-phase  and 
quadrature  components  are  resolved  from  these  measurements  to  give  the  stiffness  and 
damping  derivatives  due  to  pitch  and  plunge. 


3.3  DYNAMIC  ROLL  DERIVATIVES 

As  a  vehicle  having  aerodynamic  lifting  and  control  surfaces  rolls  about  its  longi¬ 
tudinal  axis  with  an  angular  velocity  P  .  points  on  the  surfaces  follow  a  helix  path. 
Variations  in  spenwise  local  angle  of  attack  of  the  vehicle  aerodynamic  surfaces  alter 
the  local  normal  and  drag  force  loading  distribution  and  inclination  of  the  lift  vector 
In  the  case  of  a  winged  vehicle  having  tail  surfaces,  rolling  at  various  rates  produces 
a  change  in  symmetry  of  the  trailing  vortex  sheet,  and  thus  influences  the  aerodynamic 
characteristics  of  the  tail  surface.  These  factors  have  an  influence  on  derivatives 
due  to  rolling. 

3. 3. 1  Theory 

3.3.1. 1  Free  Rotation 

The  damping- in-roll  derivative,  C{p  .  can  be  measured  using  the  free-rotation 
method  where  the  model  is  free  to  rotate  under  the  influence  of  an  initial  rolling 
velocity,  and  the  subsequent  decay  in  the  roll  r~*e  is  measured.  Another  version  of 
the  free  rotation  technique  consists  of  forced  continuous  rotation  by  a  known  moment 
imparted  by  deflecting  a  model  fin. 

Although  some  of  the  equations  describing  free  rotation  are  included  in  Section 
3.2. 1.3.  the  derivation  is  included  here  for  completeness  and  in  a  different  format. 

For  a  freely  spinning  axisymmetric  body,  the  differential  equation  of  the  free-rolling 
motion  may  be  derived  by  equating  the  rate  of  change  of  the  angular  momentum  to  the 
applied  damping  moment: 


where  the  roll -damping  moment.  MXpP  ,  consists  of  the  contribution  due  to  aerodynamic 
damping  and  the  tare  damping  due  to  friction  in  the  bearing  system.  It  is  assumed  here 
that  these  moments  are  linear  functions  of  angular  velocity.  The  mechanical  damping 
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is  evaluated  by  a  series  of  roll-decay  tests  in  a  vacuum  and  the  results  are  subtracted 
from  the  wind-on  data  to  obtain  the  aerodynamic  roll  damping.  In  applying  this  method, 
the  tacit  assumption  is  made  that  the  decay  mechanism  during  the  vacuum  tests  is 
associated  entirely  with  friction  damping  and  is  not  dependent  on  aerodynamic  loading 
and  model  dynamics.  With  tare  damping  included.  Equation  (88)  becomes 

IyP  “  t(MXp)a  *  '  (8®) 

This  equation  is  a  linear  differential  equation  of  the  first  order  with  the  solution 

P  =  Aect  .  (90) 

where  A  is  a  constant  of  integration  and 


C 


(Mxp^a  f  (Mxp\ 


The  roll -damping-moment  parameter  may  be  obtained  by  taking  the  logarithm  of  the 
ratio  of  the  solution  (90)  at  any  two  points.  1  and  2.  of  the  P  versus  t  history 
as  follows: 


For  an  exponentially  damped  system  of  the  type  described  in  Equation  (90).  MXp 
may  be  determined  using  the  above  relation  over  any  time  interval  (t2  -  tx)  of  the  P 
versus  t  history.  If  the  damping-moment  parameter  is  not  constant,  as  assumed  in 
Equation  (90).  the  P  versus  t  curve  deviates  from  being  exactly  exponentially 
damped.  For  cases  where  this  deviation  is  not  large,  the  above  solution  may  be  applied 
over  portions  of  the  P  versus  t  history  which  are  approximately  exponentially 
damped. 

The  aerodynamic  portion  of  the  total  damping- moment  parameter  is  determined  as 
follows: 

‘°«e  <Wl 

.  J. 

In  coefficient  form,  the  roll-damping  derivative  may  be  written  as 
C,p  =  3C,/3(pd/2V)  =  (MXp),(2V„/  qja2) 

The  spin  decay  technique  can  be  used  to  determine  both  the  roll-damping  moment  and 
the  aerodynamic  rolling  moment  due  to  deflected  model  fins,  grooves  in  the  model  surface, 
etc.  For  a  finned  model,  the  equation  of  moticn  is  obtained  by  modifying  Equation  (89) 
as  follows: 

IxP  =  V*"^*4*  • 

■here  MxS*d*  Is  the  rolling  aoment  due  to  (In  deflection  t* 


log,  (P;/P .)' 

.  *»  -  *i  J» 


(91) 
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For  simplicity  in  the  derivation,  no  distinction  will  be  made  here  between  aero¬ 
dynamic  damping  moments  and  damping  moments  due  to  friction,  since  the  wind-off  cor¬ 
rections  for  tare  damping  may  be  introduced  in  the  final  equations.  The  solution  for 
Equation  (91).  where  Pi  is  the  initial  rolling  velocity  and  tj  is  the  initial 
time  is  given  as 


+  (Mvn/Iv)(t  -  tn) 

MXpP  ♦  -  (MXpP,  ♦  MxS,!»*)e  **  X 


(92) 


For  the  steady-state  roll  condition,  p  =  0  and  P  -  Pc  ,  Equation  (91)  yields 

5 

P  =  -  (Mx;*/MXp)o* 


(93) 


Using  this.  Equation  (66)  becomes 


P  -  P.  -  (P,  -  P„)e 


Here  again,  taking  the  logarithm  of  the  ratio  of  the  above  solution  at  any  two 
points,  1  and  2,  of  the  P  versus  t  history  yields 

MXp  =  Ix  loge  UP,  -  P8)/(P,  -  Ps)]/(tj  -  t,)  .  (94) 


If  the  roll-rate  history  deviates  from  being  exponentially  damped,  the  solution  (94) 
may  be  applied  over  intervals  near  the  steady-state  condition  (Pg)  which  are  approxi¬ 
mately  exponentially  damped. 

The  definitions  of  the  roll-damping  coefficient  derivative,  Clp  ,  and  the  rolling- 
moment  coefficient  dernativ**  due  to  fin  deflection,  5*  ,  are  respectively 


?C,  3(Pd  2V)  =  (M„p)t(2-V'Q „M;) 

(95) 

i*  =  ?C,  =  M^/q^Ad 

(96) 

Substituting  Equation  (93)  into  Equation  (95). 

c!p  T  fix  <*:  -  ‘I*'  ifaAd2]  loge  !(Pj  -  Ps)  (P.  -  P5)] 
and,  from  Equations  (93)  and  (96). 

CM.  -  (P5  2V„)C|p  . 

The  steady-state  rolling  velocity  is  the  asymptotic  value  achieved  over  a  period  of 
time  Bearing  friction  is  evaluated  by  obtaining  air-off  roll  histories  under  vacuum 
conditions,  and  the  value  of  M  measured  under  these  conditions  is  subtracted  from 
the  values  obtained  during  tunnel  air-on  conditions 


3  3.1  2  Forced  flotation 

Derivatives  due  to  rolling,  such  as  CJp  .  may  be  measured  using  the  method  of 
forced  rotation  which  is  similar  to  the  forced-oscillation  method  A  steady  rolling 
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notion  is  imparted  to  the  model  by  means  of  a  drive  system,  and  the  damping  in  roll 
can  be  measured  between  the  drive  system  and  model.  Since  p  =  0  for  steady  rolling 
motion,  the  roll-damping  moment  is  given  by 

«Xp  =  -«x/p,  • 

where  H%  is  the  input  torque  and  Ps  is  the  steady-state  roll  velocity.  In  co¬ 
efficient  form,  the  roll-damping  derivative  is  given  as 

C!p  =  -  WxV'P.^Ad2  . 

As  was  noted  in  Section  3.1.3,  tare  damping  due  to  friction  in  the  system  is  evaluated 
with  air-off  tests  to  obtain  a  correction  for  the  air-on  data. 

Another  version  of  the  forced-rotation  technique  which  is  convenient  for  obtaining 
the  derivatives  Cjp  and  Cnp  is  the  method  of  forced  oscillation  in  roll  at  constant 
amplitude  with  freedom  to  respond  in  yaw.  The  equations  of  motion  for  yaw  and  roll 
which  describe  the  motion  are  derived  in  Reference  77.  Because  they  are  for  a  very 
specialized  case,  they  are  not  presented  here;  however,  the  system  will  be  described 
in  Section  3.  3.  2. 


3.3.2  Representative  Systems 
3.3.2. 1  Free  flotation 

A  system  used  by  NOL70,  which  has  a  roll  drive  system  external  to  the  model,  is 
shown  in  Figure  87.  The  motor  drive  consists  of  a  sliding  vane  air  motor  which  is 
coupled  to  the  sting  by  a  magnetic  clutch.  After  the  model  is  up  to  speed,  the  clutch 
is  released  and  the  decay  history  recorded  by  a  transistor-type  tachometer.  Since  the 
roll  damping  is  influenced  by  the  bearing  system,  special  consideration  is  given  to 
bearings  design  to  minimize  friction  and  thus  Improve  the  accuracy  of  the  system.  Gas 
bearings  discussed  in  Section  3. 1.1.4  are  particularly  well  suited  for  this  application 
because  of  the  small  friction  in  relation  to  ball  bearings79. 

The  external  bearing  and  drive  system  (Fig. 88)  has  tie  advantage  that  it  can  be 
easily  used  on  small  models  and  models  of  various  sizes  and  shapes.  The  disadvantage 
is  that,  because  of  its  location,  the  gas  bearing  must  be  designed  for  moment  as  well 
as  radial  load. 

The  AEDC-VKF  free-rotatlon,  roll-damping  system,  which  is  located  within  the  model, 
is  shown  in  Figure  89.  The  model  rotates  on  a  gas  bearing  and  is  driven  up  to  speed 
with  an  air  turbine.  After  the  model  has  been  driven  to  the  required  spin  rate,  an 
air-activated  clutch  is  released  and  the  model  is  free  to  rotate.  Rotational  speed 
is  measured  with  a  magnetic  pickup. 

The  magnetic  model  suspension  systems  described  in  Section  3. 2. 2. 3  may  be  used  to 
obtain  damping- 1..  roil  measurements  not  influenced  by  the  presence  of  a  support  system. 
Tliis  is  dene  by  forcing  the  model  to  roll  at  a  known  rate,  using  an  alternating  current 
for  generating  the  roll  forcing  function.  Damping  in  roll  is  determined  from  the 
logarithmic  decrement  following  termination  of  the  roll  forcing  function.  The  wind-off 
data  Includes  eddy  current  damping,  which  must  be  subtracted  from  the  wind-on  roll¬ 
damping  measurements. 
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3. 3.  2. 2  Forced  Rotation 

An  example  of  a  system80  for  forced,  continuous-rotation,  damping- in-roll  tests  is 
shown  in  Figure  90.  The  balance  has  an  electric  motor  which  drives  the  model  through 
a  planetary  gear  system.  The  stationary  ring  gear  mounted  in  the  gear  system  is  re¬ 
strained  from  rotation  by  an  instrumented  cantilever  beam.  Torque  to  overcome  model 
aerodynamic  damping  and  friction  damping  in  the  system  is  proportional  to  the  restraint 
provided  by  the  ring  gear.  Tests  in  a  vacuum  provide  tare  values  for  correcting  the 
wind-on  data. 

A  two-degree-of- freedom  system,  in  which  oscillations  in  roll  are  forced,  has  been 
used  by  NASA  to  obtain  rolling  derivatives  and  directional  stability  and  damping-in¬ 
yaw  data7'.  The  system  has  been  used  for  yaw  amplitudes  up  to  2  degrees  and  forced 
at  oscillation  frequencies  up  to  20  c/s.  A  schematic  and  pictorial  representation  of 
the  system  is  shown  in  Figures  91  and  92  to  illustrate  the  modes  of  model  motion.  By 
use  of  flexures  and  ball  bearings,  the  model  is  allowed  two  angular  degrees  of  freedom  - 
rolling  and  yawing  -  rolling  being  a  forced  oscillation  and  yawing  being  free.  Harmonic, 
forced  oscillation  in  roll  is  produced  by  means  of  rotary  motion  of  an  eccentric  and 
oscillatory  motion  of  the  drive  flexure  indicated  in  Figure  91.  Variations  in  the 
amplitude  of  the  oscillatory  rolling  motion  are  provided  by  the  amount  of  eccentricity 
in  the  drive  system.  It  is  transmitted  to  the  torque  tube,  which  is  supported  at  the 
front  end  by  roll  flexures  and  the  aft  end  by  a  ball  bearing.  Thus  the  only  contact 
between  the  model  and  the  support  structure  is  through  the  roll  flexure  system.  This 
was  necessary  because  of  the  need  to  avoid  the  undesirable  damping  characteristics  of 
ball  bearings. 

Strain  gages  in  the  system  at  locations  shown  in  Figure  92  are  used  to  measure  angular 
roll  and  yaw  positions  of  the  model,  and  the  roll  input  torque  applied  to  the  model  to 
maintain  forced  oscillations  in  roll  at  a  constant  amplitude.  The  amplitudes  and  phase 
angles  of  the  strain  gage  outputs  are  obtained  by  passing  the  signals  through  a  sine- 
cosine  resolver  located  at  the  drive  motor  and  gear  driven  by  the  motor  shaft.  The  re¬ 
solver  outputs  are  directly  proportional  to  the  in-phase  and  out-of-phase  components 
of  the  strain  gage  signals.  By  vectorial  summation  of  these  quantities,  oscillation 
amplitudes  are  obtained  as  well  as  the  relationships  necessary  to  determine  the  forcing 
torque  and  yaw  position  phase  angles.  An  electronic  counter  operated  by  a  portion  of 
the  oscillating.  r"»ll  position,  strain  gage  signal  is  used  to  measure  oscillation 
frequency. 

With  the  model  rigidly  restrained  against  motion  in  yaw,  measurements  of  the  roll 
posi4  .  torque  required  for  constant  amplitude  roll  oscillations  and  the  roll-input - 
torque  phase  angle  are  made  at  various  frequencies  by  means  of  the  resolvers.  These 
data,  obtained  for  both  wind-on  and  wind-off  (near  vacuum)  conditions,  provide 
clp  +  *C ip  at  zero  angle  of  attack  and  at  angles  of  attack. 

With  the  model  unrestrained  against  yawing  motion,  the  model  is  oscillated  in  roll 
at  the  yaw  natural  frequency  to  increase  the  yawing  amplitude.  When  the  desired  ampli¬ 
tude  is  attained,  the  drive  motor  is  stopped  abruptly  and  the  time  history  of  the  de¬ 
caying  yawing  motion  is  recorded.  Wind-on  and  wind-off  data  provide  Cnr  -  CD/g  and 

Cn/3  ' 
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Although  the  system  can  be  operated  in  two  degrees  of  freedom  to  obtain  cross  de¬ 
rivatives  Cjr  -  and  C,*  ,  the  experience  related  in  Reference  77  indicates  that, 
where  the  system  was  used  in  this  mode,  ihe  results  were  unreliable. 

The  detailed  derivations  and  data  reduction  equations  used  with  this  system  for 
measuring  dynamic  lateral  stability  derivatives  are  presented  in  Reference  77. 


3.4  MAGNUS  EFFECT 

The  “Magnus  force  and  moment”,  “Magnus  effect”,  and  "Magnus  characteristics"  are 
synonymous  terms  which  are  given  to  the  aerodynamic  forces  and  moments  that  result  in 
the  plane  perpendicular  to  the  flow  when  a  body  is  spinning  in  the  cross  flow  produced 
at  angles  of  attack.  The  name,  Magnus,  was  given  to  the  phenomenon  in  honor  of  the 
Berlin  physicist,  Professor  G. Magnus,  who  around  1850  proved  with  a  series  of  experi¬ 
ments  the  existence  of  a  Magnus  force81.  Luchuk82  gives  a  thorough  historical  account 
of  the  Magnus  effect. 

Magnus  force  data  have  become  more  important  with  the  more  general  use  of  spin  to 
produce  stabilization  and  to  overcome  manufacturing  asymmetries  of  finned  projectiles. 
Dynamic  instability  may  result  on  a  spinning  projectile  if  the  Magnus  moment  becomes 
too  large83.  Also,  the  usual  assumption  that  Magnus  effects  may  be  neglected  because 
the  spin  rate  is  small  is  sometimes  erroneous84. 

3.4.1  Types  of  Magnus  Forces 

There  are  several  different  types  of  Magnus  forces  and  moments  produced  by  various 
aerodynamic  mechanisms.  The  different  Magnus  effects  are  usually  classified  as  those 
caused  by  the  body,  by  canted  fins,  by  driven  fins,  by  base  pressure,  and  by  vorticity 
from  forward  wing  panels  acting  on  the  tail  surfaces.  These  different  types  of  Magnus 
effects  will  be  discussed  briefly  and,  where  practical,  illustrated  sketches  are  given 
(Fig. 93). 

The  body  Magnus  effect  at  low  angles  of  attack  is  caused  by  the  shifting  of  the 
thicker  boundary  layer  on  the  leeward  side  in  the  direction  of  the  spin85  (Fig. 93(a)). 
Thus  the  body  shape  is  effectively  altered  and  is  asymmetric.  At  large  angles  of 
attack  where  the  leeward  side  flow  is  separated,  the  spin  shifts  the  separated  portion 
to  an  asymmetric  shape  (Fig. 83(b)).  The  body  Magnus  force  both  at  small  anL  large 
angles  of  attack  is  in  the  negative  Y  direction  for  positive  (clockwise  looking 
upstream)  spin  rates.  Generally,  the  force  acts  on  the  rear  part  of  the  body. 

On  positive-spinning  canted  fins  (Fig. 93(c))  at  an  angle  of  attack,  each  fin,  as  it 
rotates  through  the  leeward  side,  will  lose  part  of  its  force  because  of  the  blanketing 
of  the  fin  by  the  leeward  wake86.  At  small  angles  of  attack  there  will  be  a  net  nega¬ 
tive  side  force  which  will  tend  back  toward  zero  as  the  angle  of  attack  is  increased. 
Also  on  canted  fins,  there  is  a  Magnus  moment,  which  has  been  pointed  out  by  Benton87, 
caused  by  the  axial  component  of  fin  normal  force. 

On  driven  fins  at  zero  cant  angle,  each  fin  loses  part  of  its  lift  distribution  as 
it  rotates  through  the  leeward  wake  (Fig. 93(d)).  Thus  there  is  a  net  side  force  in 
the  positive  Y  direction86. 
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A  Magnus  moment  can  be  established  on  tne  fins,  caused  by  unequal  fin  base  pressures. 
This  moment  is  small  and  can  only  exist  if  the  boundary  layer  on  the  fins  is  laminar. 
Platou84  gives  a  thorough  explanation  of  this  fin  base  pressure  effect. 

Benton88  shows  that  wing- tail  interference  can  cause  large  Magnus  effects  on  a  finned 
missile  whose  wings  are  deflected  inco  an  aileron  setting.  Basically,  the  side  loads 
are  caused  by  the  vortices  from  the  wind  interacting  with  the  tail  surfaces. 

3.4.2  Test  Procedure 

Determining  the  Magnus  effects  by  wind  tunnel  tests  necessitates  the  measurement  of 
the  loads  on  the  model  while  it  is  spinning  at  an  angle  of  attack.  It  is  a  well-known 
fact,  as  illustrated  in  Figure  94  (Ref. 89),  that  the  Magnus  effects  are  dependent  on 
the  spin  parameter  (Pb/2V).  Therefore,  in  addition  to  the  usual  parameters  such  as 
Mach  number,  Reynolds  number,  and  model  attitude,  the  spin  parameter  (Pb/2V)  must  also 
be  monitored. 

The  spin  parameter  for  the  ground  tests  must  be  matched  with  the  full-scale  flight 
tests.  For  most  wind  tunnel  facilities,  and  especially  high  Mach  numbers,  the  model 
size  is  small,  and  the  velocity  is  much  lower  than  flight  velocity  because  enthalpy 
is  not  simulated.  Therefore  much  higher  roll  rates  are  required  for  the  wind  tunnel 
model  than  for  the  full-scale  vehicle  to  properly  simulate  the  flight  spin  parameter. 
These  required  high  model  roll  rates,  as  will  be  noted  later,  often  introduce  special 
problems  in  the  experimental  work. 

Tests  by  Platou89  at  Mach  2  at  angles  of  attack  up  to  5  degrees  showed  that  the 
Magnus  force  increased  with  grit  size  (Fig. 95).  Other  investigations84’  90  also  have 
shown  that  the  Magnus  effects  are  dependent  on  the  nature  of  the  boundary  layer,  i.e. 
whether  it  is  laminar,  turbulent,  or  transitional.  Therefore,  since  it  is  usually 
not  possible  to  keep  the  Reynolds  number  of  the  ground  tests  the  same  as  the  full- 
scale  flight  test  because  of  tunnel  limitations,  additional  consideration  of  the  in¬ 
fluence  of  Reynolds  number  and  boundary  layer  transition,  to  aid  in  the  application 
of  ground  test  results  to  flight  data,  is  required. 

The  test  method  employed  usually  depends  on  the  drive  system  of  the  test  u.it.  The 
different  test  units  will  be  discussed  in  a  later  section.  For  a  motor  drive  system, 
the  data  can  be  recorded  at  the  various  angles  of  attack  while  holding  the  rotational 
speed  constant,  or  the  speed  range  can  be  covered  for  each  angle  of  attack.  With  a 
turbine  drive  system,  the  data  can  be  recorded  continuously  as  the  spin  decays  while 
holding  other  test  conditions  constant.  In  the  blowdown  tunnels,  where  the  data  re¬ 
cording  time  is  limited,  data  have  been  recorded  while  increasing  the  rotational  speed 
with  the  turbine91.  When  using  the  Mpower-up”  technique,  careful  examination  should 
be  made  to  ensure  that  the  data  are  not  influenced  by  the  air  exhaust  from  the  turbine. 

Problems  arise  because  of  gyroscopic  effects  when  data  are  taken  continuously  while 
varying  the  angle  of  attack  with  the  model  spinning.  Although  attempts  have  been  made 
to  eliminate  the  gyroscopic  forces  present  by  corrective  data  reduction82,  it  is  best 
to  avoid  this  test  technique  if  possible  because  of  the  uncertainty  involved. 

3.4.3  Representative  Systems  and  Drive  Units 

Basically,  the  majority  of  the  Magnus  force  test  systems  have  the  model  mounting 
on  ball  bearings  whose  inner  races  are  fixed  to  an  inner  shell  which  mounts  directly 
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or  indirectly  to  a  four-component  balance.  A  drive  unit,  either  an  electrical  motor 
or  an  air  turbine,  is  used  to  spin  the  model  and  normally  is  an  integral  part  of  the 
system  enclosed  by  the  model.  Tachometer  units,  consisting  of  permanent  magnets  mounted 
on  the  rotating  unit  and  a  wire-wound  pole  mounted  on  the  stationary  portion  of  the 
assembly,  are  usually  used  to  measure  the  model  spin  rates.  As  the  model  rotates,  the 
combination  provides  a  signal  which  is  proportional  to  the  spin  rate.  Typical  systems 
for  conducting  Magnus  force  tests  are  shown  in  Figures  96  to  100. 

Rotational  speeds  on  the  order  of  100,000  r.p.m.  are  sometimes  required  to  simulate 
the  full-scale  flight  spin  rate  arameter  (Pb/2Vaj).  Such  speeds  are  very  difficult, 
if  not  impossible,  to  achieve  wich  a  reasonable  size  electric  motor.  Motors  which  are 
available,  consistent  with  model  size  and  power  requirements,  are  very  expensive  and 
require  positive  internal  cooling.  Cooling  provides  added  complications  to  the  system 
because  of  routing  the  coolant  lines.  If  the  motor  is  the  connecting  link  between  the 
balance  and  model,  extreme  care  must  be  taken  to  ensure  that  the  coolant  and  power 
lines  are  not  stiff  enough  to  take  some  load  off  the  balance  or  to  create  small  inter¬ 
ference  side  loads  that  would  influence  the  accuracy  of  the  data.  Also,  with  large 
motors  some  unbalance  in  the  rotor  may  create  vibration  problems.  However,  for  large 
models  with  fins  which  have  high  damping  characteristics,  the  variable  frequency  motors 
are  often  the  best  system  and  frequently  simplify  the  design. 

A  typical  Magnus  force  test  installation92,  using  an  electric  drive  motor,  is  shown 
in  Figure  96.  The  electric  motor  and  gear  box  is  an  integral  part  of  the  model  and 
is  the  connecting  link  between  the  stationary  balance  and  the  rotating  model  through 
a  flexible  coupling  at  the  upstream  end.  The  inner  shell  mounted  to  the  balance 
supports  the  two  ball  bearings  on  which  the  model  rotates. 

High  model  spin  rates,  not  generally  available  with  electric  drive  motors,  have 
been  obtained  by  NOL,  BRL,  and  WADD  through  the  use  of  an  air  turbine  housed  within 
the  model.  The  air  turbines  can  be  made  smaller  than  commercially  available  electric 
motors,  and  weights  are  a  fraction  of  those  for  equivalent  electric  motors,  thus 
eliminating  low  resonant  frequencies  which  can  interfere  with  strain  gage  measurements. 
Also,  they  are  more  versatile  for  designing  in  that  there  is  no  direct  linkage  required 
as  with  an  electric  motor  where  there  must  be  a  direct  drive  between  model  and  motor. 
Procedures  for  design  of  small  air  turbines  for  Magnus  force  balances  are  given  in 
Reference  93. 

Magnus  force  test  systems  utilizing  air  turbine  drive  units  are  shown  in  Figures  97 
to  100.  For  the  systems  shown  in  Figures  97  to  99,  the  model  is  mounted  on  ball  bearings 
which  transmit  loading  to  a  strain  gage  balance,  whereas  in  the  system  shown  in  Figure 
100  a  gas  bearing  is  used  for  model  support.  All  of  these  have  the  air  turbine  near 
the  base  of  the  model,  with  the  air  exhausting  out  of  the  model  base  into  the  wake 
region.  Some  designs  have  been  developed  where  the  turbine  is  forward  of  the  balance 
and  a  return  loop  is  used  for  the  air  exhaust.  The  return  loop  system  has  the  advantage 
that  tests  can  be  made  while  increasing  speed  with  the  turbine  without  fear  of  an  effect 
of  the  air  exhaust  on  the  Magnus  measurements.  However,  the  turbine,  when  enclosed 
for  the  return  loop,  is  out  of  necessity  smaller  for  a  given  size  model,  and  therefore 
will  produce  less  torque.  In  addition,  problems  may  arise  when  air  is  passed  through 
the  center  of  the  balance  because  of  the  changing  balance  temperature,  thus  reducing 
the  accuracy  of  the  data. 
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The  system  shown  in  Figure  98  is  a  design  used  by  WADD  and  is  similar  to  one  deve¬ 
loped  by  BRL.  Rotational  speeds  up  to  90,000  r.p.m.  have  been  obtained  with  this  unit. 
The  rotating  assembly,  consisting  of  the  model  shell,  turbine  assembly,  and  tachometer 
permanent  magnets,  mounts  on  two  ball  bearings  which  support  the  model  on  a  strain 
gage  balance.  More  detailed  information  on  this  system  may  be  found  in  Reference  94. 

The  system  shown  in  Figure  99  is  used  in  the  AEDC-VKF  supersonic  and  hypersonic 
tunnels.  The  model  mounts  on  ball  bearings  on  an  inner  shell  which  fits  over  the 
balance.  The  turbine  is  at  the  rear  of  the  balance  and,  after  it  has  powered  the 
model  to  the  required  spin  rate,  it  can  be  disconnected  from  the  model  by  a  sliding 
clutch.  The  balance  is  a  four-component  moment  type  incorporated  in  the  system  as  a 
separate  unit.  Thus,  different  capacity  balances  can  easily  be  used  without  changing 
the  other  parts  of  the  system. 

An  air  bearing  has  been  used  by  NOL  in  V  .ystem  shown  in  Figure  100  to  provide 
either  a  Magnus  force  or  roll -damping  test  capability  with  the  same  balance  equipment 
without  modifications.  The  model  :s  spun  up  with  an  air  turbine  drive,  and  when  the 
required  spin  rate  has  been  achieved,  the  turbine  drive  is  disconnected  from  the  model 
by  an  air-operated  clutch.  Magnus  measurements  are  made  with  a  four-component  balance 
which  is  interchangeable  to  provide  different  load  capabilities.  This  balance  is 
more  complex  than  one  normally  used  with  a  ball  bearing  support,  since  supply  air  for 
the  gas  bearing  and  balance  coolant  must  be  routed  through  the  balance. 

3.4.3. 1  Bearings 

A  critical  factor  in  the  rotating  model  Magnus  force  measuring  system  is  the  bearings, 
since  they  must  in  some  cases  be  capable  of  withstanding  rotational  speeds  of  100,000 
r.p.m.  or  more.  It  is  a  general  practice  to  mount  the  bearings  so  that  they  are  sub¬ 
jected  to  a  thrust  load  at  all  times.  This  preload  confines  the  ball  rotation  to  one 
axis  and  prevents  brinelling  of  the  races  that  could  occur  from  dynamic  axial  movement. 
Temperature  rises  in  bearings  result  from  operation  at  high  speeds  and  these  should 
be  monitored  to  given  an  indication  of  the  bearing  performances. 

Bearings  used  in  the  WADD  Magnus  force  balance  system  are  run  at  10,000  and  20,000 
r.p.m.  and  temperatures  are  continuously  monitored  until  the  temperature  levels  off. 

A  plateau  in  temperature  at  the  forward  end  of  the  support  strut  is  an  indication  that 
the  bearings  are  in  good  condition  and  that  the  balls  have  worn  a  track  cycle. 

Temperatures  are  continually  monitored  during  tests,  and  excessive  rises  above  about 
150°F  indicate  that  the  bearings  should  be  disassembled,  rinsed  in  solvent,  and  re¬ 
greased.  Experience  indicates  that  bearing  temperatures  much  above  150°F  signal  the 
onset  of  freezing  of  the  bearing  assembly. 

As  noted  earlier,  the  Magnus  force  system  shown  in  Figure  98  was  designed  by  WADD 
using  much  of  the  experience  from  BRL.  A  series  of  tests  was  made  by  WADD  to  extend 
the  use  of  the  BRL  developed  bearings  to  rotational  speeds  considerably  above  the 
40.000  to  50,000  r.p.m.  specified  by  the  manufacturer  as  the  upper  limit. 

Bearing  failures  were  generally  found  to  occur  when  the  inner  race  temperatures 
reached  200  to  250°F.  Because  of  the  heat  generated,  carbon  particles  from  the  lubricant 
and  phenolic  retainer  were  produced,  and  under  some  conditions  the  carbon  particles 
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scored  the  races,  produced  additional  heat,  and  ultimately  caused  failure.  This  problem 
was  overcome  by  a  conditioning  procedure  devised  by  WADD  (Ref. 94)  which  consisted  of 
bearing  run  in  to  seat  the  balls  and  also  form  a  glazed  surface  on  the  phenolic  retainer 
at  critical  contact  points. 

Lubrication  of  bearings  is  dependent  lor  the  most  part  on  the  type  of  operation  they 
will  be  subjected  to.  For  intermittent  service,  where  the  bearings  are  allowed  to 
cotfl,  lubrication  consists  of  light  applications  of  Aircraft  and  Instrument  Grease, 
MIL-G- 25760  and  a  few  drops  of  Instrument  Lubricating  Oil,  MIL-L-6085A.  No  grease 
seals  are  used,  since  a  portion  of  the  grease  remains  on  the  bearings  for  intermittent 
operations.  Conditioning  of  the  bearing  in  this  manner  has  resulted  in  a  life  span 
of  10  to  50  runs  with  cleaning  between  every  10  to  20  spinups. 

Continuous  operation  at  high  rotational  speeds  of  70,000  to  90,000  r.p.m.  requires 
continuous  bearing  cooling  to  prevent  excessive  heating  and  bearing  freeze  up.  The 
WADD  system  shown  in  Figure  98  was  designed  with  a  bleed  line  from  the  main  air  supply 
to  provide  cooling  air  to  the  bearings.  Reduction  of  the  front  bearing  temperature 
was  effectively  achieved;  however,  cooling  of  the  rear  bearing  was  ineffective.  The 
bearing  temperatures  are  shown  in  Figure  101  for  various  throttle  value  settings  on 
the  air  cooling  air  supply.  Although  this  method  of  cooling  is  effective,  it  reduced 
the  number  of  runs  between  lubrications  because  the  lubricant  was  carried  away  with 
the  cooling  air.  A  final  development  in  the  system  included  mist  lubrication  by  an 
alemite  system  to  maintain  lubrication  while  cooling  the  bearings.  Details  of  the 
lubrication  system  are  shown  in  Figure  102. 

Bearings  with  phenolic  retainers  were  used  exclusively  in  the  WADD  Magnus  force 
system98,  since  these  were  found  to  possess  the  best  properties  for  rapid  accelerations 
and  high  speeds.  Radial-type  ball  bearings  provided  the  best  combination  with  the 
inner  race  spring-loaded.  A  precision  bearing  with  controlled  minimum  play  was  used 
for  the  front  bearing  and  a  class  3  bearing  for  the  rear  bearing.  Stainless  steel 
bearings  were  used  at  the  rear  station  because  the  rear  bearing  operated  at  higher 
temperatures  and  greater  dissipation  of  heat  was  achieved  with  stainless  steel  bearings. 

J.4.J.2  Model  Construction 

Dynamic  conditions  can  exist  at  high  rotational  speeds  which  render  Magnus  force 
data  unreliable,  increase  bearing  wear,  cause  fouling  in  the  air  turbine  drive  system, 
and  also  reduce  turbine  performance.  These  adverse  conditions  may  be  eliminated  or 
minimized  by  dynamic  balancing  of  the  model.  The  balancing  should  be  done  carefully 
and  to  a  high  degree  of  accuracy,  since  only  a  very  small  unbalance  can  significantly 
influence  the  accuracy  of  the  data. 

The  ideal  model  is  one  with  low  weight  and  high  strength  and  designed  to  maintain 
good  fits  regardless  of  temperature  gradients,  normal  loading,  and  spin  rates.  A 
material  combination  used  frequently  is  75ST  aluminum  and  heat-  dated  stainless 
steel.  Ibis  combination  has  low  weight  and  high  strength,  and  good  fits  may  be 
maintained.  It  is  particularly  important  that  a  minimum  of  parts  be  used  and  that  all 
mating  parts  be  made  to  close  tolerances,  so  that  dynamic  balancing  will  not  be  re¬ 
quired  after  each  assembly. 
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3. 4. 3. 3  Strain  Gage  Balance 

Four- component  strain  gage  balances  are  required  to  record  forces  and  moments 
generated  in  the  pitch  and  yaw  planes  by  model  spin.  Since  Magnus  forces  are  at  the 
maximum  1/10  to  1/20  of  the  normal  force  and  can  be  much  smaller,  the  Magnus  balance 
must  be  stiff  enough  to  withstand  large  normal-force  loads  while  at  the  same  time 
maintain  the  sensitivity  required  to  accurately  measure  the  small  side  loads.  NOL 
developed  a  balance  with  eccentric  side  beams  so  that  when  the  balance  is  subject  to 
yaw  loads  a  secondary  bending  moment  is  induced  in  the  eccentric  column  which  acts  as 
a  mechanical  amplifier79.  An  AEDC-VKF  balance  of  this  type  is  shown  in  Figure  103(b). 
The  AEDC-VKF  balance  has  a  normal  force  to  side  force  capacity  ratio  of  10.  Balances 
have  also  been  developed  which  have  extremely  small  outrigger  side  beams  that  act  as 
pure  tension  members.  Normal  force  to  side  force  capacities  up  to  25  have  been 
obtained  with  this  design  (Fig. 103(a)). 

Satisfactory  results  can  be  obtained  with  the  forementioned  methods  utilizing  con¬ 
ventional  wire  strain  gages.  However,  an  attractive  method  for  improvement  of  Magnus 
force  measurements  is  to  use  semiconductor  strain  gages  as  substitutes  for  wire-wound 
gages.  Wire  gages  have  a  gage  factor  of  2,  whereas  semiconductor  gage  factors  are  on 
the  order  of  60  to  200.  Utilization  of  semiconductor  strain  gages  would  produce 
voltage  outputs  30  to  100  times  greater  than  the  conventional  wire  gages,  since  the 
voltage  output  for  a  given  stress  is  directly  proportional  to  the  gage  factor. 

The  feasibility  of  using  semiconductor  strain  gages  for  Magnus  force  measurements 
has  been  investigated  by  Platou  of  BRL96.  Results  from  these  investigations  show  that 
the  semiconductor  strain  gage  is  very  sensitive  to  temperature.  These  adverse  effects 
can  be  diminished  by  temperature  compensation  provided  factory  matched  gages  are  used 
in  each  bridge.  In  addition,  the  usual  precautions  in  mounting,  bonding,  and  soldering 
connections  should  be  strictly  adhered  to  for  good  performance.  An  example  of  the 
increased  sensitivity  obtained  using  semiconductor  strain  gages  is  shown  in  Figure  104. 


3  5  PRESSURE  MEASUREMENTS 

Most  of  the  interest  in  the  last  20  years  or  more  has  been  in  the  development  of 
instrumentation  for  measuring  total  forces  and  moments  on  a  model  rather  than  local 
values.  However,  some  efforts  have  been  directed  toward  measuring  forces  and  moments 
on  segmented  wings  at  supersonic  speeds97, 98  to  obtain  section  static  aerodynamic 
characteristics.  Distributions  of  forces  and  moments,  as  well  as  the  total  loading, 
are  of  considerable  value,  especially  to  validate  theories.  Detailed  study  of  aero¬ 
dynamic  damping  may  be  enhanced  by  knowledge  of  local  loading  especially  for  the  cases 
of  models  having  flow  separations,  boundary  layer  transition  which  varies  with  angle 
of  attack,  and  models  constructed  with  low  temperature  eolation  simulators  releasing 
gas  to  the  boundary  layer. 

Measurements  of  detailed  pressure  distributions  on  an  oscillating  model  and  integra¬ 
tion  of  the  distributions  provide  a  means  to  obtain  local  loading  relative  to  the 
motion  of  the  model  and  in-phase  and  quadrature  lift  and  moment. 

With  the  exception  of  some  unreported  tests  at  AEDC-VKF  on  an  8-degree  cone  at  Mach 

4  and  measurements  of  pressure  distributions  on  an  oscillating  wedge  at  Mach  1.4  and 
1.8  (Ref.  99),  most  of  the  experimental  work100'102  has  been  conducted  at  low  subsonic 
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and  transonic  speeds.  The  authors  reporting  on  these  tests  have  described  the  trans¬ 
ducer  installations,  test  techniques,  and  data  reduction  in  detail,  and  therefore  their 
work  has  been  used  liberally  in  preparing  the  following  sections. 

3.5.1  Transducer  Installation 

The  most  important  considerations  in  developing  a  system  for  measuring  pressures  on 
an  oscillating  model  are  the  transducer  characteristics  and  installation. 

Some  of  the  desirable  characteristics  for  a  transducer  to  measure  time-varying 
pressures  on  oscillating  models,  in  addition  to  resolving  power  are103, 104: 

(a)  Transducer  outputs  should  be  insensitive  to  the  inertia  loads  imposed  under 
oscillatory  conditions. 

(b)  Transducer  outputs  should  be  linear  with  applied  pressure. 

(c)  The  transducer  should  be  snail  to  allow  installation  of  several  units  spaced 
for  adequate  coverage  and  installation  close  to  the  measuring  point  to  avoid 
a  lag  in  response. 

(d)  The  resonant  frequencies  of  the  transducer  components  should  be  well  above  the 
frequency  range  for  pressure  measurements. 

(e)  The  diameter  of  the  pressure  receiver  should  be  smaller  than  the  half  wave 
length  of  the  pressure  wave. 

(f)  The  transducer  should  be  insensitive  to  environmental  temperature  change. 

In  practice,  there  may  not  be  one  transducer  with  characteristics  which  embody  all 
of  these  requirements  for  a  particular  application,  and  therefore  a  practical  desitn 
must  be  developed  with  the  sacrifice  of  some  of  the  desirable  features.  In  all  cases, 
however,  the  transducer  should  be  selected  giving  major  weight  to  both  high  resolving 
power  and  high  frequency  response. 

The  next  important  consideration  following  selection  of  a  transducer  using  these 
guidelines  is  the  installation.  The  method  of  installation  is  governed  to  a  large 
extent  by  the  size  and  geometry  of  the  model  and  the  transducer  response  characteristics. 
The  transducer  should  be  located  near  enough  to  the  surface  to  accurately  measure 
rapidly  varying  pressures.  Ideally,  it  is  most  desirable  to  have  the  transducer  dia¬ 
phragm  right  at  the  measuring  point  to  have  zero  tube  length  and  transducer  volume. 

This,  however,  is  not  very  feasible  for  most  installations,  and,  therefore,  the  trans¬ 
ducer  must  be  connected  to  the  pressure  orifice  by  a  connecting  tube.  Transducers 
used  are  usually  the  differential-pressure  type  which  require  a  reference  pressure 
system  on  one  side  of  the  transducer  diaphragm. 

Connecting  tubing  between  the  orifice  and  transducer  must  be  selected  with  proper 
regard  to  the  important  variables  -  namely  frequency  response  and  phase  shift. 

3. 5.  2  Pressure  System  Response 

In  the  range  of  typical  tubing  lengths  ranging  from  0. !  to  6  in. ,  tubing  diameters 
ranging  from  0.015  to  0.125  in.  and  transducer  volumes  from  0.01  to  0.1  cubic  in., 
Davis104  notes  that  the  viscous  and  thermal  effects  predominate,  with  a  dependency  on 
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the  frequency  of  oscillation.  An  example  of  the  magnitude  pressure  response  and  phase 
angle  variation  with  frequency  from  Reference  104  is  presented  in  Figure  105  to  illus¬ 
trate  these  effects. 

A  method  is  presented  by  Davis104  for  calculating  the  dynamic  response  of  a  system 
and  phase  angle  for  a  sinusoidal  pressure  input.  The  assumptions  used  in  the  formu¬ 
lation,  based  on  the  work  of  References  105,  106  and  107,  include 

(a)  constant  pressure  over  the  tube  cross  section, 

(b)  constant  pressure  over  the  gage  volume, 

(c)  laminar  flow  throughout  the  system, 

(d)  small  pressure  perturbations. 

The  expression  for  the  pressure  at  the  transducer  or  gage  is  given  by 


£■ 
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This  equation  should  only  be  used  as  a  rough  guide,  since  the  pressure  system  has 
been  represented  by  a  differential  equation  with  frequency-dependent  coefficients. 

Exact  characteristics  of  a  pressure  system  should  be  determined  by  calibration. 

3.5.3  Calibration  Technique 

Calibration  of  a  transducer  tubing  system  which  is  to  be  used  for  unsteady  pressure 
measurements  is  necessary  because  of  the  numerous  uncertainties  which  exist  in  the 
analytical  treatment. 

The  dynamic  response  calibration  of  a  pressure  transducer  and  tubing  system  is 
characterized  by  frequency  response  curves  such  as  those  in  Figure  105  showing  magnitude 
and  phase  angle.  Specialized  calibration  equipment  is  necessary  for  calibrations  which 
must  be  made  over  a  wide  range  of  frequencies.  Techniques  for  calibrating  microphones 
are  directly  applicable  to  calibration  of  pressure  transducers  and  are  used  here  to 
outline  general  specifications  for  pressure  transducer  calibration  equipment104,108. 

(1)  The  pressure  oscillation  amplitude  should  be  constant  over  long  periods  and 
at  different  frequencies  and  easily  variable  over  the  range  of  interest. 

(11)  The  frequency  of  oscillation  should  be  continuously  variable  over  the  frequency 
range  of  interest. 

(ill)  Local  conditions  such  as  mean  pressure  and  temperature  should  be  variable 
over  the  range  of  Interest  and  known. 

(iv)  The  wave  form  of  the  pressure  oscillation  should  be  relatively  free  of  harmonics. 

(v)  The  calibration  equipment  should  have  a  standard  pickup  near  the  point  of 
measurement. 
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For  dynamic  calibrations  at  frequencies  of  25  c/s  and  lower,  a  system  of  the  type 
shown  in  Figure  106  will  provide  a  sinusoidal,  oscillating  pressure.  Oscillation  of 
the  pump  plunger  by  motion  of  the  motor  drive  cam  pressurizes  the  diaphragm  and  pro¬ 
duces  an  oscillating  pressure  input  to  the  transducer.  The  phase  lag  between  the 
pressure  application  and  transducer  response  is  determined  from  components  of  the 
transducer  output  in  phase  and  in  quadrature  with  the  motion  of  the  plunger.  The  mean 
pressure  in  the  calibration  may  be  adjusted  for  the  pressure  range  of  interest. 

For  frequencies  between  10  and  100  c/s,  a  piston- type  calibrator  may  be  used; 
however,  at  higher  frequencies  up  to  5000,  Davis104  recommends  a  cam- type  pulsator  of 
the  type  shown  in  Figure  107.  The  pulsator  chamber  is  supplied  with  high  pressure  air 
which  exits  through  an  opening  controlled  by  the  rotating  cam.  The  space  between  the 
cam  and  the  exit  is  adjusted  until  the  wave  form  is  sinusoidal.  Reference  pressure 
gages  and  the  transducer  and  tube  system  to  be  calibrated  are  mounted  to  sense  the 
pressure  within  the  chamber  of  the  pulsator. 
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Orientation  of  coordinate  systems 


Multipiece  crossed  flexures 


Cross-flexure 


Pig.  5  Cross- flexure,  crucifom  and  torque  tube  pivots 


Cantilever  Support 


All  dimensions  are  in  inches. 


Double  End  Support  Cantilever  Support 


Pic.  6  SB  ell  scale  instrument  flexure  pivot 
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Pig.  9  Ball  bearing  pivot 


Model  Mounting  Bulkhead 


Pi*. 10  High  amplitude  forced  oscillation  balance 
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Outer  Ring 


Shield 
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End  Plate 


Mounting  Pads 


All  dimensions  are  in  inches. 


Source  AEDC-VKF,  Ref.  19 


Pig.  12  Gas  bearing  load  characteristics 


ft  -  lb/sec\  w  in4  /iLJWsecX 
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Pig. 13  Damping  characteristics  of  AEDC-VXF  gas  bearing 


Source  HSO,  Ref.  23 
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Fig.  18  Half-models  and  reflection  planes 


(a)  Mgj  =  2. 5  ,  d8/d  =  0.4  .  a  =  ±2  deg 


Pig. 19  Support  interference  effects 


Source  AEDC-VKF,  Ref.  28 
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(b)  M*  =  2.5  .  d8/d  =  0.8  .  a  =  ±1.5  deg 
Fig. 19  Support  Interference  effects 
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Pig. 20  Support  interference  effects 
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Pig. 21  Support  interference  effects 
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Pig. 23  Model  balance  and  support  assembly  for  Hotshot  (AEDC-VKF  tunnel  F)  tunnel 


Fig. 24  Free  oscillation  balance  system 
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Rod  Diameter  or  Thickness,  In. 
Flared  Afterbody  Configuration 


Rod  Diameter  or  Thickness.,  in. 

Straight  Afterbody  Configuration  (with  Nose  Roughness) 


Pig. 27  Damping  derivatives  for  the  flared  and  straight  afterbody  configurations 
measured  in  conjunction  with  a  series  of  simulated  transverse  rod 
configurations,  forced  oscillation  tests,  amplitude  1.5  deg,  Mach  5 
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15.5° 


Mach  Number 

Pig. 30  Contribution  of  aerodynamic,  still-air,  and  structural  damping  to 

total  damping 


Fig. 31  Sketch  of  the  vacuum  system  used  to  evaluate  still-air  damping 
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Fig. 32  Variation  of  structural  damping  with  static  pressure 


Mach  Number 


Fig. 33  Variation  of  the  ratio  of  model  still-air  damping  to  model  aerodynamic 

damping  with  Mach  number 
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Pig. 34  Variation  of  the  angular  viscous- damping-moment  parameter  with  angular 

restoring-moment  parameter 


Pig. 35  Tricyclic  pitching  and  yawing 
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Source  Ref.  41 

Fig.  36  Typical  free  flight  epicyclic  pitching  and  yawing  motions  of  a  nonrolling 

statically  stable  symmetrical  missile 


Source  Ref.  41 

Fig. 37  typical  free  flight  epicyclic  pitching  and  yawing  motions  of  a  rolling 

statically  stable  symmetrical  missile 
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Forced  Oscillation  Off  of  Resonance 


Forced  Vibration  at  Resonance 


Pig.  42  Phase  diagraa 


Fig. 43  High  frequency  free  oscillation  balance  for  Hotshot  tunnel  tests 


Fig. 43  Continued 
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FiK. 44  Model  displacement  system 
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free  oscillation  balance 


Fig. 47  Free-oscillation  transverse  rod  balance 
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Pig. 48  Photographs  of  the  free-oscillation  (±35  deg)  transverse  rod  balance 
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Pig. 49  Free-oscillation  balance 
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Fig.  50  Transverse  rod  balance  for  low  density  hypersonic  flow  tunnel 


Fig. 51 


Photograph  of  the  assembled  balance 
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Fig. 52  Small  scale  free-oscillation  balance 
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Source  LTV,  Ref.  51 
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Fig. 53  Basic  mass  injection  balance  and  porous  model 


Fig. 54  Dynamic  stability  balance  suitable  for  gas  injection 
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Fig. 55  Oltronix  dampometer 


Source  Oltronix,  Ref.  53 
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Typical  calibrated  slotted 


disks  and  principle  of  dampometer  operation 
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Source  NAE,  Ref.  54 


Pig. 57  Free-oscillation  data  reduction 
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Fig. 58  Sting  balance  assembly  of  the  forced  oscillation  balance  system 
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Linkage 
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Fig. 59  Snail  amplitude  (±3  dog),  forced  oscillation  balance 
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on  dynamic  stability  testing 
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Pig. 62  Control  panel  for  forced-oscillation  system 


g.  63  Continued 
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Pig. 64  System  block  diagram  for  forced-oscillation  balance 


ASSEMBLED  OSCILLA 
BALANCE  MECHANIS 
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Pig. 65  Photograph  of  the  forward  portion  of  the  forced-oscillation  balance  nechani! 


Fig. 66  Schematic  view  of  model  and  driving-system  components 


Pig. 67  Details  of  oscillating  mechanism 


Pig. 68  Block  diagram  of  electronic  circuits  used  to  measure  model  displacement 

applied  moment 


Schematic  arrangement  for  yaw-tidealip  oacillationa. 


Source  RAE,  Ref.  61 
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Forced-oscillation  dynastic  stability  balance 
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Pig. 70  Yaw-roll  balance  and  accelerometers 


Fig. 71  Forced-oscillation  balance  with  electromagnetic  excitation 


Pivot  Axis 


Fig. 72  Cutaway  view  of  torsional  flexure  pivot  with  strain  gages  for 
forced-oscillation  balance  with  electromagnetic  excitation 
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Pig.  73  If&gnetic  ■odel  suspension  sjstei 


Pig. 74  Cone  model  for  magnetic  suspension 


Pig. 75  Basic  L  configuration  magnetic  suspension  system 
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Pig. 77  Magnetic  suspension  working  section.  7  in.  x  7  in.  hypersonic  wind  tunnel 


Pig.  78  Schematic  dlagraa  of  the  suspension  signet  array  of  a  six-component 
■agnetlc  wind  tunnel  balance  and  suspension  systea 
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VERTICAL  SUSPENSION  SPRING 
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Pig. 81  Apparatus  for  elevator  oscillation 
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Model  A 


Pig. 83  Sketch  of  control  surface  plucking  Bechaniam 
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Fig.  85  Free  oscllUtlon  »lth  autooatlcally  recycled  feedback  excitation 


Dynamic  Steady-Roll  Balance,  Internal  Details 
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Pig.  90  Dynaalc  steady-roll  balance 
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Pig. 91  Forced-oscillation  balance 
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pie- 92  Pictorial  drawing,  forced-oscillation  balance 
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Boundary  layer  Is  shifted  to 
asymmetric  shape  by  spin. 


<- — Cy  H 


(a)  Magnus  force  at  low  angles  of  attack 


Separated  region  Is  shifted 
to  asymmetric  shape  by  spin. 


«•—  Cy  H 


(b)  Magnus  force  at  large  angles  of  attack 


Pig. 93  Illustrations  of  the  various  Magnus  effects 


Effect  of  spin  parameter  on  the  Magnus  characteristics 
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Pig.  95 


Effect  of  boundary  layer  on 


the  Magnus  force 


Pit- 97  *Mnaa  force  »m«.  turbine  dri 
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Pig. 98  Magnus  force  system,  turbine  drive  within  the  model 
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Pig. 100  Magnus- roll  damping  balance  with  air  bearings 


Bearing  Teap 


Pig.  101  Variation  of  bearing  te*>erature  with  spin  rate  for  bypass  air  cooling 


lubricating  system  and  other  instrumentation 


Semiconductor  Gage 


Frequency ,  f ,  cps 


Fig.  105  Experimental  response  magnitude  and  phase-angle  calibration  for  1.38  in.  of 
0.030  in.  nominal-diameter  tubing  connected  to  a  gage  volume  of  0.02  cu. in. 
for  air  at  sea- level  conditions 
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Source  Ref.  101 


Pig. 106  Puap  for  sinusoidal  oscillating  pressure 
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